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EXTREMELY  LOW  FREQUENCY  (ELF)  PROPAGATION  FORMULAS 
FOR  DIPOLE  SOURCES  RADIATING  IN  A 
SPHERICAL  EARTH-IONOSPHERE  WAVEGUIDE 


1.  INTRODUCTION 


Extremely  low  frequency  (ELF)  propagation  formulas  for  dipole  sources  radiating  in  a 
spherical  earth-ionosphere  waveguide  (figure  1-1)  have  been  derived  by  various  authors. 
Developments  of  such  formulas  are  provided  in  the  texts  written  by  J.  R.  Wait  (reference  1)  and 
J.  Galejs  (reference  2).  These  spherical  waveguide  formulas  have  been  derived  for  a  uniform  earth 
and  ionosphere  and,  depending  on  the  field  component,  have  a  range  dependence  characterized  by 
either  the  Legendre  function  of  the  first  kind  of  complex  degree  and  order  zero  or  one  of  its  first 
two  derivatives.  In  the  ELF  band  (defined  here  as  30  to  300  Hz),  because  the  ionospheric 
reflection  height  is  less  than  one-half  of  a  free-space  wavelength,  the  only  propagating  mode  in  the 
earth-ionosphere  waveguide  is  the  zeroth-order  transverse  magnetic  (TM)  mode,  which  is 
commonly  referred  to  as  the  quasi-transverse  electromagnetic  (quasi-TEM)  mode. 


Figure  1-1.  Spherical  Coordinate  System  Description  of  the  Earth-Ionosphere  Waveguide 


To  predict  the  fields  more  accurately  at  ranges  closer  to  the  source  or  to  account  for  more 
complicated  boundary  conditions,  various  investigators  have  incorporated  the  earth-flattening 
approximation  into  their  ELF  propagation  formulas.  In  this  approximation,  the  Legendre 
function  range  dependence  is  approximated  by  the  product  of  a  Hankel  function  and  a  curvature 
correction  term  (spherical  earth  spreading  factor).  As  a  result,  the  fields  are  derived  from  a 
planar  earth-ionosphere  waveguide  model  and  then  multiplied  by  the  curvature  correction  term. 
For  example,  Bannister  (reference  3)  has  derived  ELF  propagation  formulas  based  on  the  earth¬ 
flattening  approximation  that  extend  the  results  of  Wait  and  Galejs  to  closer  ranges  from  the 
source.  Whereas  the  spherical  waveguide  formulas  given  by  Wait  and  Galejs  are  valid  for  ranges 
greater  than  approximately  three  ionospheric  reflection  heights  from  the  source,  Bannister’s 
formulas  are  valid  in  the  quasi-nearfield  range,  which  is  defined  as  the  range  where  the 
measurement  distance  is  greater  than  an  earth  wavelength,  but  much  less  than  a  free-space 
wavelength.  However,  Bannister’s  approximate  formulas  are  not  valid  at  field  points  close  to  the 
antipode  where  the  simple  spherical  focusing  factor  fails. 

To  account  for  the  anisotropic  surface  impedance  in  the  vicinity  of  a  horizontal  electric 
dipole  (HED)  at  ELF,  Wolkoff  and  Kraimer  (references  4  and  5)  have  listed  propagation 
formulas  that  are  modifications  of  Bannister’s  HED  formulas.  Wolkoff  and  Kraimer’s  formulas 
account  for  the  anisotropic  ground  through  the  use  of  two  complex-valued  antenna  pattern 
factors.  These  antenna  pattern  factors  are  unique  for  a  given  HED  and  must  be  determined  from 
near-field  measurements  of  the  antenna.  Wolkoff  and  Kraimer  have  determined  the  antenna 
pattern  factors  for  each  of  the  U.S.  Navy’s  ELF  transmitting  antennas  (reference  4).  Wolkoff 
and  Kraimer’s  propagation  formulas  have  been  formally  derived  from  Bannister’s  formulas  by 
Casey  (reference  6)  through  use  of  the  reciprocity  theorem. 

For  prediction  of  the  ELF  fields  from  dipole  sources  at  antipodal  ranges,  propagation 
formulas  that  are  valid  out  to  approximately  20  Mm  from  the  source  must  be  applied.  In  a  recent 
report  (reference  7),  approximate  formulas  for  a  HED  source  that  contain  the  proper  range 
dependence  at  antipodal  ranges  in  a  spherical  earth-ionosphere  waveguide,  referred  to  as 
“antipode-centered  propagation  formulas,”  were  derived.  These  HED  formulas  are  based  on 
Burrow’s  simple  parallel-plate  waveguide  approximation  of  the  earth-ionosphere  waveguide 
(reference  8)  and  include  a  curvature  correction  factor.  In  reference  7,  the  antipode-centered 
propagation  formulas  were  compared  with  Bannister’s  HED  formulas  (direct  and  indirect  great- 
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circle  path  fields  were  combined)  under  various  propagation  conditions,  where  both  the  source 
and  field  points  are  located  on  the  earth’s  surface.  The  results  showed  that  Bannister’s  vertical 
electric  field  and  radial  magnetic  field  formulas  (magnitude  only)  agree  to  within  1  dB  of  the 
corresponding  antipode-centered  formulas  for  ranges  greater  than  0.97  Mm  to  1 . 1 3  Mm  from  the 
antipode,  depending  on  the  propagation  conditions.  In  addition,  Bannister’s  azimuthal  magnetic 
field  formula  agrees  to  within  1  dB  of  the  corresponding  antipode-centered  formula  for  ranges 
greater  than  3.17  Mm  to  3.72  Mm  from  the  antipode,  depending  on  the  propagation  conditions. 

In  this  report,  ELF  propagation  formulas  for  dipole  sources  radiating  in  a  spherical  earth- 
ionosphere  waveguide  are  derived  from  first  principles.  These  derivations  are  presented  because 
the  developments  given  by  previous  authors  were  found  to  be  difficult  to  follow.  The  formulas 
derived  here  are  based  on  the  assumptions  of  a  homogeneous,  isotropic  earth  and  a 
homogeneous,  isotropic  ionosphere  of  constant  reflection  height.  As  a  result,  the  earth  and 
ionosphere  boundaries  are  modeled  as  scalar  surface  impedances.  The  spherical  waveguide 
formulas  are  expressed  in  terms  of  series  expansions  of  TM  and  transverse  electric  (TE)  modes. 
However,  at  ELF,  the  only  mode  of  practical  interest  is  the  quasi-TEM  mode.  The  computed 
results  for  the  spherical  waveguide  formulas  presented  in  this  report  are  based  on  an  exact 
hypergeometric  series  representation  of  the  Legendre  function  of  the  first  kind.  The  propagation 
formulas  for  a  horizontal  magnetic  dipole  (HMD)  can  be  derived  from  the  HED  formulas 
through  application  of  the  duality  principle  (reference  2). 

In  addition,  in  this  report,  through  appropriate  approximations  of  the  range  dependence,  it 
will  be  shown  how  the  spherical  waveguide  formulas  (for  both  vertical  electric  dipole  (VED)  and 
HED  sources)  can  reduce  to  Bannister’s  formulas  or  to  the  antipode-centered  formulas.  It  will 
also  be  shown  how  the  quasi-TEM  spherical  waveguide  formulas  derived  for  a  HED  located  at 
the  surface  of  the  earth  can  be  modified  to  account  for  the  anisotropic  surface  impedance  in  the 
vicinity  of  the  antenna.  These  modified  spherical  waveguide  formulas  will  include  Wolkoff  and 
Kraimer’s  antenna  pattern  factors  and  will  be  useful  for  prediction  of  the  electromagnetic  (EM) 
fields  radiated  by  the  U.  S.  Navy’s  four  transmitting  antennas  at  antipodal  ranges.  Detailed 
derivations  of  approximate  formulas  for  the  Legendre  function  range  dependence  are  given  in  the 
appendices.  Comparisons  of  the  spherical  waveguide  propagation  formulas  with  both 
Bannister’s  direct  great-circle  path  and  total  field  (direct  plus  indirect  great-circle  paths) 
formulas  are  presented  for  the  surface  magnetic  field  and  vertical  electric  field  components  under 
various  propagation  conditions.  The  field  comparisons  are  presented  at  several  frequencies  for 
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both  VED  and  HED  sources  at  ranges  that  extend  from  1  Mm  from  the  source  to  the  antipode, 
where  both  the  source  and  field  points  are  located  at  the  earth’s  surface.  These  comparisons  will 
help  to  more  accurately  establish  the  maximum  ranges  over  which  Bannister’s  approximate 
formulas  are  valid. 
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2.  APPLICATION  OF  PARALLEL-PLATE  WAVEGUIDE  MODEL  FOR 
DETERMINATION  OF  APPROXIMATE  MODE  CUTOFF  FREQUENCIES 


To  obtain  a  basic  understanding  of  EM  propagation  in  a  spherical  earth-ionosphere  waveguide 
in  the  ELF  band,  it  is  helpful  to  formulate  the  problem  under  the  assumption  of  planar  boundaries. 
Such  an  approach  has  been  applied  by  Burrows  (reference  8)  and  Budden  (reference  9).  In  the 
ELF  band,  both  the  earth  and  ionosphere  appear  as  nearly  perfect  conductors  (references  2  and  8). 
Therefore,  in  this  section,  the  TEM,  TM,  and  TE  modes  are  derived  for  a  parallel-plate  waveguide 
with  perfectly  conducting  walls. 

The  analysis  of  this  section  will  show  that,  because  the  ionospheric  reflection  height  is  much 
less  than  a  free-space  wavelength  at  ELF,  the  TEM  mode  is  the  only  propagating  mode.  The 
TEM  mode  is  characterized  by  electric  and  magnetic  fields  that  are  oriented  perpendicular  to  the 
direction  of  propagation.  If  the  finite  conductivity  of  the  earth  and  ionosphere  are  considered, 
the  TEM  mode  configuration  is  perturbed  to  the  quasi- TEM  mode.  The  quasi-TEM  mode 
includes  surface  electric-  and  magnetic -field  components  that  lie  along  the  direction  of 
propagation.  Therefore,  the  quasi-TEM  mode  is  the  only  mode  of  practical  concern  at  ELF. 

The  cutoff  frequencies  for  the  dominant  TE  and  TM  modes  are  given  for  several  representative 
ionospheric  reflection  heights.  The  propagation  formulas  for  dipole  sources  radiating  in  a 
parallel-plate  waveguide  with  finitely  conducting  boundaries  have  been  derived  by  Burrows 
(reference  8)  and  will  not  be  presented  in  this  report. 

2.1  TEM  AND  TM  MODES 

Consider  the  parallel-plate  waveguide  model  of  the  earth- ionosphere  waveguide,  as  shown 
in  figure  2- 1 .  In  this  model,  the  z-direction  is  normal  to  the  planar  boundaries  and  the  x-direction 
denotes  the  direction  of  propagation.  Because  the  parallel-plate  waveguide  is  a  two-dimensional 
model,  it  is  assumed  that  there  is  no  variation  along  the  ^-direction,  i.e.,  ?  I?y  =  0  for  all  field 
components.  For  simplicity,  the  earth  and  ionosphere  are  both  assumed  to  be  perfectly 
conducting. 

Appendix  A  provides  a  derivation  of  the  EM  fields  in  terms  of  potentials.  These  formulas 
are  based  on  an  assumed  time-harmonic  dependence  of  eim,  where  a>  =  2nfis  the  angular 
frequency  (rad/s),/is  the  frequency  in  Hertz  (Hz),  and  j  =  1  .  A  time-harmonic  dependence 

will  be  assumed  in  each  of  the  field  expressions  presented  in  this  report. 
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HIGHLY  CONDUCTING 

Figure  2-1.  Parallel-Plate  Model  of  the  Earth-Ionosphere  Waveguide 


In  this  section,  the  TEM  and  TM  modes  (with  respect  to  the  z-direction)  are  derived  in  terms 
of  the  z  component  of  the  electric  Hertz  potential  n  e. ;  and  the  TE  modes  (with  respect  to  the  z- 
direction)  are  derived  in  terms  of  the  z  component  of  the  magnetic  Hertz  potential  .  The  z 
component  of  each  potential  is  chosen  in  order  to  generate  modes  that  will  correspond  with  the 
spherical  waveguide  modes  that  are  derived  in  appendix  C. 

To  derive  the  TEM  and  TM  modes  in  the  parallel-plate  waveguide,  consider  the  electric 
Hertz  potential  ite  that  is  given  as 

ne  =  z  3i*(x,z)  ,  (2-1) 

where  z  denotes  the  unit  vector  along  the  z-direction.  The  electric  Hertz  vector  is  defined  in 
terms  of  the  magnetic  vector  potential  A  in  appendix  A,  expression  (A-26a).  The  substitution  of 
this  definition  into  the  inhomogeneous  vector  Helmholtz  equation  for  the  magnetic  vector 
potential  (A-6)  yields 

VxVx^-ity  =  -Vr+-L  J  ,  (2-2) 

J(0eo 

where  k0  =  (o  fp~e~0  is  the  wave  number  in  free  space,  p0  and  e0  are  the  permeability  and 
permittivity,  respectively,  of  free  space,  &e  is  the  electric  scalar  potential,  and  J  is  the  electric 
current  density.  Note  that  expression  (2-2)  applies  to  a  point  lying  in  free  space. 
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To  eliminate  &e  from  the  above  equation,  the  Lorentz  gauge  condition  (reference  10)  is 
applied,  i.e.. 


V  •  jre  +  0e  =  0 


(2-3) 


The  substitution  of  expression  (2-3)  into  (2-2)  yields 


W2jte  +  k2jte  =  ■ 


1 


J<°£o 


J. 


(2-4) 


To  derive  the  TEM  and  TM  modes,  it  is  assumed  that  there  are  no  sources  so  that  /=  0  in 
formula  (2-4).  The  following  is  obtained  if  expression  (2-1)  is  substituted  into  equation  (2-4): 


d2ne„  d2Jt ?  . 

ljr+ljr  +  Mf=0- 


(2-5) 


Because  the  modes  are  assumed  to  propagate  along  the  x-direction  in  the  waveguide,  jie  can  be 


expressed  as 

ne.(x,z)  =  f(z)  e~Jkx 


(2-6) 


where  k  is  the  wave  number  (along  the  direction  of  propagation)  and f(z)  is  a  function  to  be 
determined.  The  substitution  of  expression  (2-6)  into  (2-5)  yields 

f"(z)  +  k2cf(z)  =  0,  ( 


where 

k20  =  k2o-k2  .  (2-8) 

In  the  above  formula,  kc  is  commonly  referred  to  as  the  cutoff  wave  number  (reference  10). 


Expressions  for  the  EM  fields  in  terms  of  the  electric-  and  magnetic -Hertz  vectors  are  given 
in  formulas  (A-27)  and  (A-28).  For  the  TEM  and  TM  modes,  the  magnetic-Hertz  vector  JCh  is 
zero.  Therefore,  the  substitution  of  equation  (2-1)  into  the  EM  field  formulas  (A-27)  and  (A-28) 
yields  the  following: 


^  K? 

1 

II 

(2-9a) 

Ez =  k2jll  , 

(2-9b) 
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and 


Hy  =  -  a>£0kjiez  •  (2-9c) 

Note  that  the  x  variation  defined  in  expression  (2-6)  has  been  accounted  for  in  the  above  field 

formulas.  The  field  components  listed  above  are  the  only  nonzero  components.  Also,  note  that 
because  the  TM  modes  are  “transverse  magnetic”  with  respect  to  the  z-direction,  Hz  =  0. 

The  solution  of  the  differential  equation  (2-7)  is  given  as 

f{z)  =  c,  cos  kcz  +  c2  sin  kcz  ,  (2- 10) 

where  c  j  and  c2  are  arbitrary  constants.  To  determine  c  j  and  c2,  the  boundary  conditions  at  the 

surfaces  of  the  earth  and  ionosphere  must  be  applied.  Because  the  earth  and  ionosphere  are 
assumed  to  be  perfectly  conducting,  the  following  boundary  conditions  apply: 

Ex(x,0)  =  0,  z  =  0,h  ■  (2-11) 

The  substitutions  of  expressions  (2-9a)  and  (2-6)  into  the  boundary  conditions  above  yields 
/'(0)  =  f\h)  =  0  •  (2-12) 

From  formula  (2-10),/'(z)  is  given  as 

f{z)  =  -  kccl  sin  kcz  +  kcc2  cos  kcz  ■  (2-13) 


The  application  of  the  boundary  condition  for/'  at  z  =  0  to  formula  (2-13)  gives 

kcc2  =  0  •  (2-14) 

Therefore,  the  solutions  of  the  above  equation  are  kc  =  0  or  c  2  =  0.  If  c2  =  0,  then  the  boundary 
condition  f'(h)  =  0  results  in 

*c»  =  X  ’  W  =  (U’2’- ■  (2‘15) 

The  index  n  in  the  above  formula  denotes  the  mode  index. 

The  substitution  of  the  above  result  for  /  (z)  into  expression  (2-6)  for  ne_  yields 
Jtez„(x,z)  =  Cj  cos(^f)  e~jkx ,  u  =  0,  1,  2, ... ,  (2-16) 
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where  k  is  given  from  formulas  (2-8)  and  (2-15)  as 

.  (2-17) 

A  mode  is  said  to  be  at  cutoff  when  k  =  0  or  equivalently,  when  kc  =  k0.  The  cutoff  wave 
numbers  for  the  TEM  («  —  0)  and  TM  modes  (n  =  1,  2, ...)  are  given  in  formula  (2-15). 

Therefore,  the  cutoff  frequencies  for  these  modes  are  given  by 

fen  —  2ft  ^ cn  ~  2h  ’  n  —  *  (2-18) 

where  c  denotes  the  speed  of  light  in  free  space.  Table  2-1  lists  the  cutoff  frequencies  for  the 
dominant  modes  in  the  parallel-plate  waveguide  as  computed  from  expression  (2-18)  for  several 
representative  ionospheric  reflection  heights.  The  table  indicates  that  the  TEM  mode  propagates  at 
all  frequencies  while  the  TMj  mode  does  not  propagate  for  frequencies  below  1.67  kHz.  Because 
the  frequency  range  of  interest  here  lies  below  300  Hz,  the  TM  modes  are  of  no  concern. 


Table  2-1.  Cutoff  Frequencies  for  TEM,  TMlt  and  TM2  Modes  for  a 
Parallel-Plate  Waveguide  with  Perfect  Electrically  Conducting 
Boundaries  for  Several  Waveguide  Heights 


/co(Hz) 

/cl  (kHz) 

/c2(kHz) 

50 

0 

3.00 

— 

75 

0 

2.00 

MEM 

90 

0 

1.67 

3.33 

The  TEM  and  TM  fields  are  determined  through  the  substitution  of  expression  (2-16)  into 
formulas  (2-9)  to  give  the  following: 


and 


£„  =  ^s in(^f)e-^,n=l,2,3,..., 
Etn  =  k 2  cos(^f)  e~jkx ,  n  =  0, 1, 2, ... , 


(2- 19a) 
(2- 19b) 


Hyn  =  -  we0  e  jkx ,  n  =  0, 1, 2, ... 


(2- 19c) 
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In  formulas  (2-19),  note  that  the  arbitrary  constant  has  been  suppressed.  It  should  be  noted 
that  Ex  vanishes  for  the  TEM  mode  («  =  0).  This  is  consistent  with  the  fact  that  the  TEM  mode 
has  no  field  component  along  the  direction  of  propagation  (reference  10).  However,  for  a 
finitely  conducting  boundary,  Ex  does  not  vanish  but  is  related  to  the  surface  magnetic  field  Hy 
through  the  surface  impedance  of  the  boundary  surface.  For  example,  for  a  finitely  conducting 
earth  with  surface  impedance  r\e,  Ex  is  related  to  Hy  at  z  =  0  as  Ex  =  -J]e  Hy. 

2.2  TE  MODES 

The  derivation  of  the  TE  modes  in  the  parallel-plate  waveguide  are  very  similar  to  the 
derivation  of  the  TM  modes.  Consider  the  magnetic  Hertz  potential  Kh  that  is  given  as 

jc h-zjr/’(x,z)  •  (2-20) 

The  magnetic  Hertz  vector  is  defined  in  terms  of  the  electric  vector  potential  F  in  expression 
(A-26b).  The  substitution  of  this  definition  into  the  inhomogeneous  vector  Helmholtz  equation 
for  the  electric  vector  potential  (A- 1 1 )  yields 

VxVxjt*- k2nh  =  -  V(ph  +  — - — M,  (2-21) 

where  <Ph  is  the  magnetic  scalar  potential  and  M  is  the  magnetic  current  density.  It  should  be 
noted  that  expression  (2-21)  applies  to  a  point  in  free  space. 

To  eliminate  <Ph  from  the  above  equation,  the  Lorentz  gauge  condition  (reference  10)  is 
applied,  i.e., 

V  •  nh  +  <Ph  =  0  •  (2-22) 

The  substitution  of  expression  (2-22)  into  (2-21)  yields 

V2ah  +  k2nh  = - —  M  .  (2-23) 

To  derive  the  TE  modes,  it  is  assumed  that  there  are  no  sources  so  that  M  =  0  in  the  above 
formula.  If  expression  (2-20)  is  substituted  into  equation  (2-23),  the  following  expression  is 
obtained: 
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(2-24) 


d2jr£  d2Jth 
dx2  +  dz2 


+  k2Ji^  =  0  . 


Because  the  modes  are  assumed  to  propagate  along  the  x-direction  in  the  waveguide,  n1}  can  be 
expressed  as 

nh-i  jrh(x,z )  =  z  g(z)  e~jkx ,  (2-25) 


where  g(z)  is  a  function  to  be  determined.  The  substitution  of  expression  (2-25)  into  (2-24) 
yields 

g"(z)  +  k2cg(z)  =  0  ,  (2-26) 

where  the  cutoff  wave  number  kc  is  defined  in  formula  (2-8). 


The  TE  mode  EM  fields  can  be  obtained  through  the  substitution  of  expression  (2-20)  into 
formulas  (A-27)  and  (A-28)  with  the  electric  Hertz  vector  ne  set  to  zero  to  yield  the  following: 


Ey  =  omjaih  , 

(2-27a) 

(2-27b) 

H*  =  *2*?  ■ 

(2-27c) 

Note  that  the  x  variation  defined  in  expression  (2-25)  has  been  accounted  for  in  the  above  field 
formulas.  The  field  components  listed  above  are  the  only  the  nonzero  components.  Also,  note 
that  because  the  TE  modes  are  transverse  electric  with  respect  to  the  z-direction,  we  have  Ez  =  0. 

The  solution  of  the  differential  equation  (2-26)  is  given  as 

g(z)  =  c3  cos  kcz  +  c4  sin  kcz  ,  (2-28) 

where  c3  and  c4  are  arbitrary  constants.  To  determine  c3  and  c4,  the  boundary  conditions  at  the 

surfaces  of  the  earth  and  ionosphere  must  be  applied.  Because  the  earth  and  ionosphere  are 
assumed  to  be  perfectly  conducting,  the  following  boundary  conditions  exist: 

£y(x, 0)  =  0,  z  =  0,h  .  (2-29) 
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The  substitutions  of  expressions  (2 -27a)  and  (2-25)  into  the  boundary  conditions  above  yields 


5(0)  =  g(h)  =  0 


(2-30) 


The  substitution  of  the  boundary  conditions  (2-30)  into  expression  (2-28)  results  in  the  following 
solution: 


gm(z)  =  c4 sin(^p)  ,m  =  1,2, 3, ... . 


(2-31) 


A  comparison  of  formulas  (2-28)  and  (2-31)  shows  that  the  cutoff  wave  number  kcm  for  the  mth 
TE  mode  is  given  as 


k  -  HML  ni  —  i  93 

Kcm  fo  5  ni  ~~  1?  w **• 


(2-32) 


The  substitution  of  the  above  result  for  gm(z)  into  formula  (2-25)  for  n1}  followed  by  another 
substitution  into  the  TE  field  expressions  (2-27)  yields 


Eym  =  (Ofxji  sin(^p)  e  Jkx ,  m=  1, 2, 3, ... , 

[mjr)e~Jkx’ m= 1’2’3’  — 


_  Jkmjr  „-jkx 


Hxm  =  ~  C°S 


and 


Hvn  =  k  sin 


{mr]e~jkx ’  m~  i’2’3—  ’ 


(2-33a) 

(2-33b) 

(2-33c) 


where  the  wave  number  k  is  given  from  formulas  (2-8)  and  (2-32)  as 

The  cutoff  frequencies  for  the  TE  modes  are  given  by 


(2-34) 


fa 


2n 


me 
2  h 


,  m  -  1,  2,  3, ... . 


(2-35) 


A  comparison  of  formula  (2-18)  with  (2-35)  indicates  that  the  cutoff  frequencies  of  the  TMn  and 
TEm  modes  are  the  same  provided  n  =  m.  Therefore,  the  cutoff  frequencies  for  the  TEj  and  TE2 
are  the  same  as  those  of  the  TMj  and  TM2  modes,  respectively,  given  in  table  2-1.  Thus, 
because  the  frequency  range  of  interest  here  lies  below  300  Hz,  the  TE  modes  are  also  of  no 
concern  in  this  report. 
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3.  ELF  PROPAGATION  FORMULAS  BASED  ON  THE 
EARTH-FLATTENING  APPROXIMATION 


To  reduce  the  mathematical  complexities  associated  with  the  problem  of  the  propagation  of 
radio  waves  over  a  spherical  earth,  Pryce  (reference  11)  introduced  the  earth-flattening 
approximation  in  which  the  problem  reduces  to  the  propagation  over  a  planar  earth  with  an 
atmosphere  having  a  modified  refractive  index.  In  his  investigation  of  the  accuracy  of  the  earth¬ 
flattening  approximation,  Pekeris  (reference  12)  assumed  that  the  range-dependence  function  can 
be  expanded  in  an  asymptotic  series  involving  increasing  negative  powers  of  the  earth’s  radius. 

In  Pekeris’  asymptotic  series  (for  the  range  dependence),  the  first  term  corresponds  to  the  model 
of  a  flat  earth  and  the  succeeding  terms  are  corrections  for  curvature.  Koo  and  Katzin  (reference 
13)  extended  the  work  of  Pryce  and  Pekeris  to  obtain  exact  differential  equations  for  the 
spherical  geometry  in  terms  of  equations  of  planar  earth  type,  resulting  in  solutions  that  are 
applicable  for  arbitrary  ranges  and  heights.  Wait  (reference  1)  has  applied  the  earth-flattening 
approximation  to  obtain  an  approximate  range  dependence  in  the  ELF  band. 

In  this  section,  the  earth-flattening  approximation  as  applied  to  the  range  dependence  of  the 
fields  is  presented.  Bannister’s  ELF  propagation  formulas,  which  incorporate  the  earth¬ 
flattening  approximation,  are  presented  for  VED  and  HED  antennas  located  on  the  surface  of  the 
earth.  In  addition,  Wolkoff  and  Kraimer’s  HED  propagation  formulas,  which  are  modifications 
of  Bannister’s  formulas  to  account  for  the  anisotropic  surface  impedance  in  the  vicinity  of  the 
source,  are  also  presented.  For  each  set  of  formulas  presented  in  this  section,  both  the  direct  and 
indirect  great-circle  path  fields  are  given  because  both  fields  must  be  combined  in  order  to 
predict  the  propagation  at  antipodal  ranges.  In  section  5  of  this  report,  Bannister’s  formulas  are 
compared  with  the  spherical  waveguide  formulas  to  determine  the  maximum  ranges  of  validity 
of  Bannister’s  formulas. 

3.1  EARTH-FLATTENING  APPROXIMATION  TO  THE 
LEGENDRE  FUNCTION  OF  THE  FIRST  KIND 

As  shown  in  appendices  F,  G,  and  H,  the  range  dependencies  of  the  EM  fields  in  the  earth- 
ionosphere  waveguide  are  expressed  in  terms  of  the  Legendre  function  of  the  first  kind, 

Pv(-  cos  6) ,  or  one  of  its  first  two  derivatives.  From  appendix  J,  the  complex  degree  v  of  the 
Legendre  function  is  related  to  the  wave  number  k  in  the  waveguide  as 

ksf}-ja=v+U2,  (3-1) 
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where  /?  =  2n IX  and  a  are  the  phase  and  attenuation  constants,  respectively,  and  X  is  the 
wavelength.  It  should  be  mentioned  that  k  corresponds  to  the  quasi-TEM  mode,  the  only 
propagating  mode  in  the  spherical  waveguide  at  ELF. 


A  mathematical  development  of  the  earth  flattening  approximation  is  given  in  appendix  K. 
From  expression  (K-33),  the  Legendre  function  is  approximated  as 


Pv(~  cos  0) 
j  sin  vji 


\[^r„ 


Ivl  »  1,  \kp\  »  1  , 


(3-2) 


where  p  =  ad  is  the  great-circle  path  distance  along  the  earth  from  the  source  point  to  the  field 
point  and  denotes  the  Hankel  function  of  the  second  kind  with  order  0.  In  the  above 
approximation,  H^\kp)  corresponds  to  the  range  dependence  for  a  flat  earth  and  the  square- 
root  term  is  the  correction  for  curvature.  Expression  (3-2)  is  sometimes  referred  to  as  the  earth¬ 
flattening  approximation  with  curvature  correction.  It  should  be  noted  that  expression  (3-2) 
differs  by  a  factor  of  j  with  the  one  given  by  Wait  (reference  1).  The  propagation  formulas 
presented  in  the  remainder  of  this  section  each  incorporate  approximation  (3-2). 


3.2  BANNISTER’S  FORMULAS 

As  mentioned  in  section  1,  Bannister  (reference  3)  has  derived  ELF  propagation  formulas  for 
dipole  sources  that  incorporate  the  earth-flattening  approximation.  These  formulas  are  valid  in 
the  quasi-near  field  range,  which  corresponds  to  measurement  distances  that  are  greater  than  an 
earth  wavelength.  Although  only  Bannister’s  formulas  for  VED  and  HED  sources  are  presented 
here,  the  field  expressions  for  an  HMD  source  can  be  readily  obtained  from  the  HED  formulas 
(reference  14)  through  replacement  of  the  electric  dipole  moment  p  with  jk/n,  where  m  is  the 
current  moment  (current-area  product)  of  the  HMD,  and  ke  is  the  wave  number  in  the  earth.  It 

should  be  added  that  the  equivalent  HMD  is  oriented  perpendicular  to  the  HED,  i.e.,  ay-directed 
HMD  is  equivalent  to  an  jc-directed  HED  located  on  the  surface  of  the  earth,  where  the  z 
direction  points  radially  outward  from  the  earth’s  surface. 

It  should  be  mentioned  that  Bannister’s  formulas  presented  here  have  been  slightly  modified 
in  order  to  compute  the  phase  as  well  as  the  magnitude  of  the  field.  The  same  formula 
modifications  were  used  by  Wolkoff  and  Kraimer  (references  4  and  5)  in  their  HED  formulas. 
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3.2.1  VED 


Consider  a  VED  located  on  the  earth’s  surface  and  oriented  along  the  positive  z  direction  in 
the  coordinate  system  defined  in  figure  3-1.  Under  the  assumption  of  a  scalar  earth  conductivity 
cre,  Bannister  (reference  3)  has  derived  propagation  formulas  that  are  valid  at  the  earth’s  surface 
(z  -  0)  for  ranges  greater  than  an  earth’s  wavelength  (i.e.,  p  >  Ae).  These  formulas  are  given  as 


Fd=  MoP 
z  2jc kap3 


Vh(t)  e-aP  +  j%  Gh(u )  (kp) 2  H$Xkp) 


pi  a 


1/2 


sin  (p  /  a ) 


(3-3  a) 


pl  a 


11/2 


and 


sin  (p  /  a) 


pl  a 


(3-3b) 


sin  (p  /  a) 


1/2 


where 


Gh(u )  =  ^  coth(w)  +  (l  -  ^]«2csch2(«)  * 
VA(r)  =  t3  coth(f)  csch2(r) » 


(3-3c) 

(3-4) 

(3-5) 


and 


p  =  Idl  (electric  dipole  moment). 


k0-T~lT  (wave  number  in  free  space), 


k=®-ja.  (wave  number  in  the  earth-ionosphere  waveguide), 


S -  'jr  -  v~  J  ~b~  (normalized  wave  number), 

Ko  Ko 


Ve  = 


JWc 


cre  +  jo)ee 


(intrinsic  impedance  of  the  earth), 


u  = 


np 

2h 


and 
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The  superscript  d  in  the  field  components  refers  to  the  direct  great-circle  path  contribution.  In 
the  above  formulas,  /  denotes  the  antenna  current,  dl  is  the  effective  length  of  the  antenna,  a  is 
the  radius  of  the  earth,  X0  is  the  wavelength  in  free  space,  £e  is  the  permittivity  of  the  earth,  v  is 

the  speed  of  propagation  in  the  earth-ionosphere  waveguide,  a  is  the  attenuation  constant  in 
Np/m,  h  is  the  ionospheric  reflection  height,  and  p  =  a# is  the  direct  great-circle  path  distance 
from  the  source  to  the  field  point.  Note  that  a  time-harmonic  dependence  of  e*m  is  assumed  in 
expressions  (3-3),  where  (0=  2rfis  the  angular  frequency  (rad/s),  and/is  the  frequency  in  Hertz. 

The  functions  Gh  and  Vh  were  included  by  Bannister  to  extend  the  range  of  the  previous  ELF 
formulas  down  to  the  quasi-nearfield  range  (i.e.,  p  >  Xe  and  p  «  A0).  Therefore,  for  an  earth 
conductivity  ae  =  2\  10'4  S/m  and  frequency/=  76  Hz,  Bannister’s  formulas  are  valid  for  ranges 
greater  than  30  km  from  the  source.  Bannister  has  noted  that  for  w  <  0.5,  Gh{u)  s  1 ;  for  t  <  0.5, 
Gh(t)  s  Vh(t)  s  1 ;  for  u  >  2.5,  Gh(u)  s  =  £ ;  for  t  >  4.5,  Vh(t)  =  0 .  Note  that  the  spherical 
earth  spreading  (curvature)  factor  [(/>/«)  /  sin  (p/a)  J  appears  in  each  VED  field  component 
expression. 


A 


Figure  3-1.  Coordinate  System  Used  by  Bannister  (Reference  3 ) 

(The  VED  and  HED  sources  are  each  located  at  the  origin  with  the  VED  oriented 
along  the  positive  z-direction  and  the  HED  along  the  positive  x-direction.) 


The  VED  field  expressions  given  above  are  associated  with  the  fields  that  propagate  along 
the  direct  great-circle  path  that  connects  the  source  and  field  points  (figure  3-2).  At  antipodal 
ranges,  the  contribution  from  the  indirect  great -circle  path  field  must  be  included  in  order  to 
account  for  the  interference  produced  by  the  two  paths.  The  indirect  great-circle  path  field 
components  for  the  VED  can  be  obtained  from  the  direct  great-circle  path  fields  in  equations 
(3-3)  through  replacement  of  the  range  p  with  the  indirect  great-circle  path  range  p,  =  2jta  -  p 
(figure  3-2)  and  with  the  inclusion  of  the  appropriate  phase  shift  as  obtained  from  table  J-l  in 
appendix  J.  As  explained  in  appendix  J,  this  added  phase  shift  of  ±  j  depends  on  the  primary 
range  dependence  of  the  field  and  therefore  varies  with  each  field  component.  Therefore,  the 
indirect  great-circle  path  VED  fields  are  given  as  follows: 


pi  =  VoP 
z  231 k0pf 


W  ap-  +  J  f  w2  HjfKkp) 


p./a 


sin  (p  /  a) 


1/2 


(3-6a) 


Pi  I  a 


and 


sin  (p  /  a) 


Pi/a 


1/2 


(3-6b) 


sin  (p  /  a) 


1/2 


(3-6c) 


Tip  JlP  1 

where  u{  =  and  t{  =  ■  The  superscript  i  in  the  field  components  refers  to  the  indirect 

great-circle  path  contribution.  For  each  component,  the  total  field  is  obtained  through  addition  of 
the  direct  and  indirect  great-circle  path  contributions,  i.e., 


Ez  =  Ef  +  E‘z, 

(3-7a) 

hi 

II 

hi 

+ 

hi 

7D 

(3-7b) 

(3-7c) 
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2a 


z 


P  +  Pj  =  2na 
P  +  P8  ~  !ta 


Figure  3-2.  The  Earth  and  the  Two  Great-Circle  Paths  to  the  Field  Point 


3.2.2  Horizontal  Electric  Dipole 


Consider  a  HED  located  on  the  earth’s  surface  and  oriented  along  the  positive  x  direction  of 
the  coordinate  system  defined  in  figure  3- 1 .  Under  the  same  assumptions  given  for  the  VED,  the 
HED  propagation  formulas  are  given  as  follows: 


p  I  a 


sin  (p  /  a) 


1/2 


cos  cp 


(3-8a) 


£d=  r^pG^t) 
p  4a)p0p 


Hf{kp)-^Hf\kp) 


p  /  a 


sin  (p  /  a) 


1/2 


cos  q)  5 


(3-8b) 


VjkpHjP) 

4(op0p2 


H^\kp) 


p  /  a 

sin  (p  /  a) 


1 3/2 

sin  cp  > 


VckPHh({) 
4  iop0p2 


HfXkp ) 


p  /  <7 

sin  (p  /  a) 


3/2 

sin  > 


(3-8c) 


(3-8d) 
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and 


K 


4mficp 


H$Xkp)-±HfXkp) 


p  /  a 


sin  (p  /  a) 


where 


1/2 


COS  (p  f 


(3-8e) 


=  Gh(t)  +  v*«) . 


(3-9) 


Bannister  has  noted  that  for  t  >  2.5,  Gh(t)  s  yjf  =  j-  ( I  and  for  t  >  4.5,  Vh(t  )  =  0  and 
2 t_p(k  A2  1  ^ 


1/2 


Hhd)  -  GhW -It-J  \~§  j  •  Note  that  the  spherical  earth  spreading  factor  [(p/a)  /  sin  (p/a)J 
appears  in  each  HED  field  component  expression  along  with  an  additional  factor  of 
[(p/a)  /  sin  (p/a)  j  appearing  in  the  and  formulas.  This  factor  of  [(p/a)  /  sin  (p/a)  j  also 
appears  in  the  corresponding  HED  spherical  waveguide  formulas  derived  in  the  appendix  H.  Note 
that  the  surface  electric  field  components  and  are  related  to  the  surface  magnetic  field 
components  and  ,  respectively,  through  the  impedance  of  the  earth  i}e. 


As  previously  mentioned,  at  antipodal  ranges,  the  contribution  from  the  indirect  great-circle 
path  field  must  be  included  in  order  to  predict  the  interference  produced  by  the  two  paths.  The 
indirect  great-circle  path  field  components  for  the  HED  can  be  obtained  from  the  direct  great- 
circle  path  fields  in  equations  (3-8)  through  replacement  of  the  range  p  with  the  indirect  great- 
circle  path  range  p,  and  with  the  inclusion  of  the  appropriate  phase  shift  as  obtained  from  table 
J-l  in  appendix  J.  As  mentioned  earlier  in  this  section,  the  added  phase  shift  of  ±  j  depends  on 
the  primary  range  dependence  of  the  field  and  therefore  varies  with  each  field  component.  Thus, 
the  indirect  great -circle  path  HED  fields  are  given  as  follows: 


_  -  VekpGh(u{) 
4p  • 


Pi /a 


sin  (p  /  a) 


1/2 


cos  cp  , 


(3- 10a) 


EP  = 


4a)p0Pi 


Pi  I  a 


sin  (p  /  a) 


1/2 


COS  (p, 


(3- 10b) 


=  jp2ekpH  h(t  j) 
*  4amopf 


Pij_ a 

sin  (p  /  a) 


3/2 

sin  <p , 


(3- 10c) 
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Pi  I  a 


sin  (p  /  a) 


3/2 


sm  (p 


(3-10d) 


and 


m<fi2p  i 


Hv  4o>p0  Pi 


«*('*> 


Pi  I  a 


sin  (p  /  g) 


1/2 


COS  (p 


(3-10e) 


where  p(,  w„  and  tj  were  previously  defined  for  the  VED.  For  each  field  component,  the  total 
field  is  obtained  through  addition  of  the  direct  and  indirect  great-circle  path  contributions  as  is 
shown  for  the  VED  fields  in  expressions  (3-7). 


3.3  FORMULAS  FOR  A  HED  ABOVE  AN  ANISOTROPIC  GROUND 

To  account  for  the  anisotropic  earth  conductivity  beneath  a  HED  antenna,  Wolkoff  and 
Kraimer  (references  4  and  5)  have  extended  Bannister’s  formulas  (3-8).  Wolkoff  and  Kraimer’s 
formulas  have  been  formally  derived  by  Casey  (reference  6).  The  coordinate  system  used  by 

Wolkoff  and  Kraimer  to  describe  the  EM  field  radiated  by  a  HED  at  a  given  point  on  the  surface 
of  the  earth  is  described  in  figure  3-3.  In  this  illustration,  Et  and  Nt  refer  to  the  geometric  (or 

true)  east  and  north  directions,  respectively,  at  the  transmitter  location,  while  Er  and  Nr  refer  to 
the  geometric  east  and  north  directions,  respectively,  at  the  receiver  location  (or  field  point). 

Note  that  the  geometric  east  and  north  directions  vary  as  functions  of  location  along  the  earth’s 
surface  and  lie  on  the  tangent  plane  to  the  spherical  earth.  In  addition,  E  and  N  correspond  to  the 
x  and  y  directions,  respectively,  in  the  local  rectangular  coordinate  system  at  each  point  on  the 
spherical  earth,  and  the  vertical  direction  z  points  radially  outward  from  the  center  of  the  earth  at 
each  point  as  illustrated  in  figure  3-4. 

In  figure  3-3,  the  range  p  is  the  direct  great-circle  path  distance  from  the  HED  center  to  the 
field  point,  (p  is  the  azimuthal  angle  measured  counterclockwise  from  true  east  at  the  transmitter, 
and  %  is  the  azimuthal  angle  measured  clockwise  from  true  north  at  the  transmitter.  The  unit 
vectors  (p,  q>,  z )  form  a  right-handed  cylindrical  coordinate  system  in  which  z  is  directed 
radially  outward  from  the  center  of  the  earth.  This  cylindrical  coordinate  system  differs  from  the 
convention  used  by  Bannister  and  others  (compare  with  figure  3-1),  who  define  the  azimuthal 
angle  (p  to  be  measured  counterclockwise  with  respect  to  the  electrical  axis  of  the  HED.  Wolkoff 
and  Kraimer  (references  4  and  5)  use  the  revised  coordinate  system  described  here  because  their 
antenna  pattern  factors  are  determined  from  magnetic  field  measurements  that  are  referenced 
with  respect  to  the  local  north  and  east  directions. 
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Figure  3-3.  Coordinate  Systems  at  the  Transmitter  and  Receiver  Locations 
Used  by  Wolkoff  and  Kraimer  (References  4  and  5) 


MQRTH  POLE 


EQUATOR 


Figure  3-4.  Local  Rectangular  Coordinate  System  (E,  N,  z)  at  a 
Given  Point  on  a  Spherical  Earth 


Wolkoff  and  Kraimer’s  propagation  formulas  for  the  direct  wave  fields  from  a  HED  are 
given  as  follows: 


Ed  JkIGhM  Hm{kp) 

‘  4 p  1 


p/ a 

1/2 

r  "j 

sin  (p  /  a) 

Tesin  cpn-Tn  cos  (f>n 

(3-1  la) 


.  klHAt)  ... 

K  =  H?\kP) 


p  4  a)fxy 


p  /  a 

3/2 

r  i 

sin  (p  /  a) 

Te  cos  <pfI  +  Tn  sin  <pH 

(3-1  lb) 


and 


Hd  = 


k2IGh(t) 


v  4  <°P0P 


H£\kp)-±HfKkp) 


p  1  a 

1/2 

sin  (p  /  a) 

- 

Te  sin  (p„-Tncos(pn 


.  (3-1 lc) 


In  the  above  formulas,  Te  and  Tn  denote  the  antenna  pattern  factors  in  the  east  and  north 

directions,  respectively.  These  complex-valued  quantities  account  for  the  antenna  length  as  well 
as  the  anisotropic  ground  conductivity  beneath  the  antenna.  For  a  given  HED,  Te  and  Tn  are 

determined  through  measurement  of  the  surface  magnetic  field  at  a  location  lying  in  the  near 
field  of  the  antenna  (reference  4).  It  should  be  noted  that  in  the  above  formulas,  Te  and  Tn  have 

been  interchanged  from  their  original  definitions  given  in  reference  6.  A  detailed  discussion  of 
the  antenna  pattern  factors  is  given  in  a  technical  report  that  is  still  in  preparation  (reference  15). 


Wolkoff  and  Kraimer’s  propagation  formulas  for  the  indirect -wave  fields  from  a  HED  are 
obtained  from  the  direct  great -circle  path  fields  in  equations  (3-11)  through  replacement  of  the 
range  p  with  the  indirect  great-circle  path  range  p,  and  with  the  inclusion  of  the  appropriate 
phase  shift  of  ±  j  for  each  field  component  as  obtained  from  table  J-l  in  appendix  J.  Thus, 

Wolkoff  and  Kraimer’s  HED  indirect  wave  formulas  are  given  as  follows: 
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where  p,  is  defined  in  figure  3-2  and  w,  and  tt  are  defined  in  section  3.2.1.  For  each  field 
component,  the  total  field  is  obtained  through  addition  of  the  direct  and  indirect  great-circle  path 
contributions  as  is  shown  for  the  VED  fields  in  expressions  (3-7). 

A  comparison  of  Wolkoff  and  Kraimer’s  HED  formulas  (3-11)  with  Bannister’s  formulas 
(3-8)  yields  the  conversions  in  table  3-1. 

Table  3-1.  Conversion  Table  Relating  Bannister ’s  HED  Propagation 
Formulas  with  Those  of  Wolkoff  and  Kraimer 


Field  Component 

Bannister 

Wolkoff  and  Kraimer 

e9h9 

-  T]e  dl  cos  (p 

Te  sin  (pn  -  Tn  cos  (p„ 

Hp 

-  rje  dl  sin  (p 

Te  cos  (p„  +  Tn  sin  <pn 

To  convert  from  one  of  Bannister’s  HED  formulas  to  the  corresponding  Wolkoff  and  Kraimer 
formula,  the  quantity  in  the  second  column  of  table  3-1  is  replaced  by  the  corresponding  quantity 
given  in  the  third  column.  Table  3-1  is  a  revision  of  a  similar  table  originally  presented  by  Casey 
(reference  7).  A  comparison  of  columns  two  and  three  above  shows  that  the  antenna  length  dl  in 
Bannister’s  formulas  is  absorbed  into  the  antenna  pattern  factors. 

Numerical  results  based  on  Wolkoff  and  Kraimer’s  HED  formulas  will  not  be  presented  in 
this  report.  For  results  based  on  these  formulas,  the  reader  is  referred  to  a  forthcoming  report 
(reference  15). 
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4.  SPHERICAL  WAVEGUIDE  PROPAGATION  FORMULAS 


ELF  propagation  formulas  for  VED  and  HED  sources  that  radiate  in  a  spherical  earth- 
ionosphere  waveguide  have  been  derived  in  appendices  F  and  H,  respectively.  These  spherical 
waveguide  formulas  are  based  on  the  assumptions  of  a  uniform  isotropic  earth  and  a  uniform 
isotropic  ionosphere  of  constant  height,  where  the  waveguide  boundaries  are  represented  by 
scalar  surface  impedances.  Through  use  of  the  thin -shell  approximation  derived  in  appendix  I,  it 
is  shown  that  only  the  quasi-TEM  mode  can  propagate  in  the  30-Hz  to  300-Hz  frequency  band. 

In  this  section,  the  quasi-TEM  spherical  waveguide  formulas  for  VED  and  HED  sources  are 
presented,  where  the  source  and  observation  points  are  located  on  the  surface  of  the  earth.  From 
the  series  approximation  of  the  Legendre  function  of  the  first  kind  (see  appendix  J,  section  J.2),  it 
is  shown  that  the  HED  spherical  waveguide  formulas  reduce  to  expressions  that  are  closely 
related  to  the  antipode-centered  formulas  that  were  previously  derived  (reference  7)  for  antipodal 
ranges.  In  addition,  through  use  of  the  earth-flattening  approximation  derived  in  appendix  K,  it 
is  shown  that  the  spherical  waveguide  formulas  reduce  to  Bannister’s  formulas  given  in  section 
3.2.  Finally,  through  use  of  table  3-1,  spherical  waveguide  formulas  for  a  HED  above  an 
anisotropic  ground  are  presented. 


4.1  PROPAGATION  PARAMETERS  FOR  TM  AND  TE  MODES  AT  ELF 


The  radial  wave  numbers  for  the  quasi-TEM,  TM,  and  TE  modes  in  the  spherical  earth- 
ionosphere  waveguide  at  ELF  are  derived  in  appendix  I.  These  formulas  are  based  on  the  thin- 

shell  approximation  that  has  been  shown  (reference  2)  to  be  suitable  in  the  ELF  band.  The  radial 
wave  number  krv  has  been  defined  in  terms  of  the  ffee-space  wave  number  k0  and  the  wave 

number  k  in  the  direction  of  propagation  as 

(4-i) 

From  appendix  I,  approximate  formulas  for  the  radial  wave  number  krv  for  the  quasi-TEM 

o 

mode  and  the  radial  wave  numbers  krv  for  the  TM  modes  are  given  as 
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(4-3) 


k  clEL  + 
K"'n~  ll 


rw  ,  AAg+Al*o 


mi 


,  AJ<Ji  «!;//=  1,2,3,.... 


Similarly,  approximate  formulas  for  the  radial  wave  numbers  k  for  the  TE  modes  are  given  by 


nm 
h 


k  s 
^  -v 

m 


1  - 


7(4  +  4) 


kji 


;  w  =  1,2,3, .... 


(4-4) 


In  the  above  approximations,  it  is  assumed  that  the  spherical  waveguide  boundaries  appear  as 
nearly  perfect  conductors.  Therefore,  in  formulas  (4-2)  and  (4-3),  \Ag  |«  1  and  \Ae  |«  1,  where  Ag 
denotes  the  normalized  surface  impedance  of  the  earth  and  Ae  denotes  the  normalized  surface 
impedance  of  the  ionosphere  for  the  quasi- TEM  and  TM  modes.  Similarly,  in  formula  (4-4), 

\Ag  | «  1  and  \Ah  |  «  1,  where  A/,  denotes  the  normalized  surface  impedance  of  the  ionosphere  for 
the  TE  modes. 


For  a  given  mode,  the  wave  number  k  in  the  direction  of  propagation  is  given  from  formula 
(4-1)  as 

k  =  P~ja  ~{k20- &rV) 1/2  ,  (4-5) 

where  ft  =  271/ A  and  a  are  the  phase  and  attenuation  constants,  respectively,  and  A  is  the 
wavelength.  In  the  above  formula,  note  that  the  radial  wave  number  kn  corresponds  to  a 
particular  mode  and  is  given  by  one  of  the  approximations  (4-2),  (4-3),  or  (4-4).  A  subscript  or 
superscript  has  not  been  included  with  £  because  only  quasi-TEM  mode  propagation  formulas 
are  referred  to  in  the  remainder  of  section  4.  In  the  ELF  propagation  literature,  the  mode 
propagation  parameters  are  given  by  c/v,  a,  and  h,  where  c  is  the  speed  of  light  in  free  space, 
v  is  the  phase  velocity  of  the  mode,  and  h  is  the  ionospheric  reflection  height.  The  mode 
parameter  c/v  is  expressed  in  terms  of  the  wave  number  k  as 


c  _  R e{k}  _  l 

v“  K  "V 


(4-6) 


Tables  4-la  and  4-lb  provide  a  listing  of  the  phase  velocity  ratio  (c/v)  and  attenuation  a , 
respectively,  for  equally  spaced  frequencies  across  the  30-Hz  to  300-Hz  frequency  band  under 
daytime  propagation  conditions  with  h  =  50  km,  earth  conductivity  og  =  10"3  S/m,  and 
ionospheric  conductivity  cr,  =  10'5  S/m.  These  tables  are  based  on  the  thin-shell  approximation 
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and  are  computed  from  formulas  (4-2)  to  (4-6).  Note  that  the  attenuation,  expressed  in  decibels 
per  megameter  (dB/Mm),  can  be  obtained  from  a(Np/m)  through  the  relation  1  Np  =  20  log10(e) 
dB  =  8.686  dB. 


Table  4-la.  Phase  Velocity  Ratio  c/v  for  the  Dominant  Modes  in  a  Spherical 
Earth-Ionosphere  Waveguide  at  ELF  Frequencies  Under  Daytime 
Conditions  with  h=50  km,  ag  - 10~3  S/m,  and  Gi  - 10'5  S/m 


Table  4- lb  shows  that  the  quasi- TEM  mode  is  the  only  propagating  mode  in  the  ELF  band 
with  an  attenuation  rate  of  approximately  0.75  dB/Mm  at  30  Hz  that  increases  with  frequency  to 
2.63  dB/Mm  at  300  Hz.  Table  4- lb  also  shows  that  the  attenuation  rates  of  the  TM  and  TE 
modes  are  too  high  to  be  of  any  practical  concern  across  the  30-Hz  to  300-Hz  frequency  band. 
Note  that  the  attenuation  increases  with  the  order  of  the  mode.  In  addition,  the  attenuation  rates 
of  the  TE  j  and  TE2  modes  are  less  than  those  of  the  TM  j  and  TM2  modes,  respectively.  Table 
4- la  indicates  that  the  phase  velocity  of  the  quasi-TEM  mode  monotonically  increases  with 
frequency  over  the  30-Hz  to  300-Hz  band. 
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Table  4-lb.  Attenuation  a  in  dB/Mm  for  the  Dominant  Modes  in  a  Spherical 
Earth-Ionosphere  Waveguide  at  ELF  Frequencies  Under  Daytime 
Conditions  with  h=  50  km,  og  =  10~3  S/m,  and  <T;  =  10~s  S/m 


Frequency 

(Hz) 

Quasi-TEM 

TMi 

tm2 

TE, 

te2 

30 

0.7543 

545.7 

1091 

390.6 

781.3 

60 

1.110 

545.6 

1091 

430.4 

861.0 

90 

1.385 

545.4 

1091 

449.5 

899.5 

120 

1.617 

545.2 

1091 

461.3 

923.3 

150 

1.822 

544.9 

1091 

469.4 

939.9 

180 

2.007 

544.5 

1091 

475.3 

952.3 

210 

2.178 

544.1 

1091 

479.8 

961.9 

240 

2.337 

543.6 

1090 

483.4 

969.7 

270 

2.486 

543.1 

1090 

486.2 

976.1 

300 

2.628 

542.5 

1090 

488.5 

981.5 

Tables  4-2a  and  4-2b  list  the  phase  velocity  ratio  and  attenuation  rate,  respectively,  for 
equally  spaced  frequencies  across  the  30-Hz  to  300-Hz  ELF  band  under  nighttime  propagation 
conditions  with  h  =  15  km,  Gg  =  10"^  S/m,  and  G[  =  10'^  S/m.  The  numbers  in  these  tables  are 

also  based  on  the  thin-shell  approximation.  Table  4-2b  shows  that  the  quasi-TEM  mode  is  the 
only  propagating  mode  in  the  ELF  band  at  nighttime  with  an  attenuation  rate  of  approximately 
0.53  dB/Mm  at  30  Hz  that  monotonically  increases  with  frequency  to  1.78  dB/Mm  at  300  Hz  . 
Note  that  the  attenuation  rates  of  the  TM  and  TE  modes  are  again  too  high  to  be  of  any  practical 
concern.  Table  4-2a  indicates  that  the  phase  velocity  of  the  quasi-TEM  mode  monotonically 
increases  with  frequency  over  the  tabulated  frequency  band. 
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A  comparison  of  tables  4-1  with  tables  4-2  shows  that  the  mode  attenuation  rates  are  less 
at  nighttime.  In  addition,  the  phase  velocity  of  the  quasi-TEM  mode  is  greater  at  nighttime. 
Because  the  quasi-TEM  mode  is  the  only  propagating  mode  in  the  30-Hz  to  300-Hz  band,  the 
remainder  of  this  report  will  focus  on  the  quasi-TEM  mode. 


Table  4-2a.  Phase  Velocity  Ratio  c/v  for  the  Dominant  Modes  in  a  Spherical 
Earth-Ionosphere  Waveguide  at  ELF  Frequencies  Under  Nighttime 
Conditions  with  h  =  75  km,  <jg  =  10'3  S/m,  and  <rf  =  10'5  S/m 

Frequency  guasi.TEM  TM,  TM,  TE,  TE, 


30 


1.106 

3.199  xlO’3 

1.599  xlO’3 

9.360 

1.075 

4.525  x  lO'3 

2.262x1 0'3 

3.729 

1.061 

5.546  xlO'3 

2.771  x  lO’3 

2.143 

1.053 

6.409  xlO'3 

3.200  xlO'3 

1.439 

150 

1.048 

7.172  xlO'3 

3.579  xlO'3 

1.054 

180 

1.043 

7.867  x  lO'3 

3.921  x  lO’3 

0.8163 

210 

1.040 

8.510  xlO'3 

4.237  xlO'3 

0.6572 

240 

1.038 

9.113  xlO'3 

4.532  xlO'3 

0.5445 

270 

1.035 

9.685  xlO'3 

4.809x1  O’3 

0.4612 

2.103 


1.627 


1.023  xlO’2  5.072  x  10’3 


Table  4-2b.  Attenuation  a  in  dB/Mm  for  the  Dominant  Modes  in  a  Spherical 
Earth-Ionosphere  Waveguide  at  ELF  Frequencies  Under  Nighttime 
Conditions  with  h  =  75  km,  cg  =  Id3  S/m,  and  <T;  =  10~5  S/m 


Frequency 

(Hz) 

Quasi-TEM 

TMj 

tm2 

TEj 

te2 

30 

0.5263 

363.8 

727.6 

290.9 

581.9 

60 

0.7655 

363.6 

727.6 

310.7 

621.6 

90 

0.9496 

363.4 

727.4 

120 

1.105 

363.0 

727.3 

325.1 

651.2 

150 

1.242 

362.6 

727.1 

328.6 

658.8 

180 

1.366 

362.1 

726.8 

331.0 

664.4 

1.480 

726.5 

332.7 

668.6 

1.586 

360.8 

726.2 

333.9 

672.0 

270 

1.686 

360.0 

725.8 

334.6 

674.7 

300 

1.780 

359.2 

725.3 

335.1 

676.8 

4.2  QUASI-TEM  FIELDS 
4.2.1  VED 


The  spherical  waveguide  propagation  formulas  for  a  VED  located  at  an  arbitrary  height 
above  the  surface  of  the  earth  are  derived  in  appendix  F.  These  formulas  simplify  considerably 
for  the  case  when  both  the  source  and  observation  points  are  located  on  the  surface  of  the  earth. 
Under  this  assumption,  from  appendix  F,  section  F.2,  the  quasi-TEM  fields  for  a  VED  are  given 
as  follows: 


-.ve 


-,ve 
'■0  : 


and 


h  :  = 


XT  1  (JL 

4 kjia2 

U  U  '  p 

sin  van  1 

'ioP\  1  ;i 

Aha 

sin  v0n  36 

.  P 

At _  *„ 

1  0  p  < 

Pv  (—  cos  6)  , 


(4- 7a) 


(4- 7b) 


(4-7c) 
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The  reader  is  referred  to  the  spherical  coordinate  system  ( r ,  0,  (p)  described  in  figure  1-1  with 
unit  vectors  ( r,d,q> )  illustrated  in  figure  4-1.  In  these  formulas,  the  VED  is  located  at  (r,  0)  =  (a, 
0)  and  the  observation  point  is  at  (r,  0,  (p)  =  (a,  0,  (p).  Note  that  because  of  azimuthal  symmetry, 
the  above  formulas  are  independent  of  (p.  In  addition,  the  excitation  factor  Ae0  has  been  replaced 
by  0.5,  a  suitable  approximation  at  ELF,  as  is  shown  in  appendix  I,  section  1.2.  In  the  above 
formulas,  Ag  denotes  the  normalized  surface  impedance  of  the  earth  and  Pv°  is  the  Legendre 
function  of  the  first  kind  of  degree  v0  and  order  zero,  where  v0  is  the  n  =  0  solution  of  the  TM- 
mode  characteristic  equation. 


z 


From  appendix  J,  section  J.l,  v0  is  related  to  the  quasi-TEM  wave  number  k  as 


K~  a 


(4-8) 


As  discussed  in  appendix  K,  if  |  v0  \ »  1,  then 
v0+H2  s  [v0(vo+  1)]1/2  . 


(4-9) 


The  above  result  is  generally  valid  in  the  ELF  band.  Therefore,  under  this  condition,  the  wave 
number  k  can  be  approximated  as 
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(4-10) 


k  s 


[v0(v0+1)] 


1/2 


a 


For  a  homogeneous  and  isotropic  earth,  the  normalized  surface  impedance  A„  is  given  by 


4 

6  ~  Vo  ’ 


(4-11) 


where  the  intrinsic  impedance  of  the  earth  rje  is  given  in  section  3.2.1. 


4.2.2  Horizontal  Electric  Dipole 

The  spherical  waveguide  propagation  formulas  for  a  HED  located  at  an  arbitrary  height 
above  the  surface  of  the  earth  are  derived  in  appendix  H.  As  with  the  VED  formulas,  the  HED 
formulas  simplify  considerably  for  the  case  when  both  the  source  and  observation  points  are 
located  on  the  surface  of  the  earth.  Under  this  assumption,  from  appendix  H,  section  H.2,  the 
quasi-TEM  fields  for  a  HED  are  given  as  follows: 


vhe  ^o^gP  i  d 

E.  = - jjr - ^ —  -th  Pv  (-  cos  6 )  cos  w  , 

T  Aha  sinv^jr  86  V  '  T  ’ 

(4- 12a) 

rhe  PoVcAlP  1  1  d2  n  , 

E»=-  4 h  vJy<>  +  1)  sin  K(-  cos  0)  COS  <f . 

(4- 12b) 

E%=Ja”\f!P  ,  1  »  ‘  *P>(-  cos0)sin«>, 

v  4 h  vo(v0  +  1)  sin  vQn  sin  6  86  vc 

(4- 12c) 

zjhe  JkoAgP  1  1  1  d  r>  4  m 

=  a  i*  ,  .  — Tn  (-  cos  0)  sin  w , 

6  4^  v0(v0  +  1)  sin  vcJt  sin  6  86  v<, 

(4-1 2d) 

and 

H  =  . , 5  *  .  *  2  F  (—  cos  0)  cos  cp  . 

v  4^  v0(v0  +  1)  sin  v0Jt  302  ° 

(4-12e) 

The  reader  is  again  referred  to  the  spherical  coordinate  system  (r,  6 ,  ^)  described  in  figure  1-1. 
In  these  formulas,  the  HED  is  located  at  (r,  0)  =  (a,0)  and  the  observation  point  is  at  (r,  6,  <p)  = 
(a,  9,  <p).  Note  that  the  above  formulas  have  the  same  azimuthal  dependencies  as  Bannister’s 
HED  formulas  that  are  listed  in  section  3.2.2.  As  in  the  VED  case,  the  excitation  factor  Ae0  has 
been  approximated  by  0.5.  The  complex  degree  v0  in  the  above  HED  formulas  corresponds  to 

the  n  =  0  solution  of  the  TM-mode  characteristic  equation. 
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4.3  ANTIPODE-CENTERED  FORMULAS 


In  order  to  predict  the  EM  fields  produced  by  a  HED  at  antipodal  ranges,  the  author 
(reference  7)  derived  some  simple  approximate  formulas  that  are  based  on  an  earth-flattening 
type  of  approximation  that  has  the  proper  range  dependence  in  the  vicinity  of  the  antipode  and 
accounts  for  spherical  curvature  at  further  distances  from  the  antipode.  The  resulting  “antipode- 
centered”  formulas  have  range  dependencies  that  are  similar  in  form  to  the  corresponding  ones 

that  incorporate  the  earth-flattening  approximation  except  that  the  Hankel  function  is  replaced  by 
a  Bessel  function  of  the  first  kind  and  the  range  argument  is  replaced  by  pa,  where  pa  =  a(n  -  6) 

is  the  great-circle  path  distance  from  the  antipode  to  the  field  point  (figure  3-2).  The  range 
dependencies  in  the  antipode-centered  formulas  are  derived  from  an  approximate  series 
representation  of  Pv(-  cos  0)  that  is  suitable  in  the  vicinity  of  the  antipode  (0  =  n).  From  the 
primary  terms  in  this  series,  approximate  formulas  for  Pv(-  cos  6)  and  its  first  two  derivatives 
with  respect  to  6  are  derived  in  appendix  J,  section  J.2.  These  formulas  are  given  as  follows: 


Pv(-  cos  6)  =  J0(kpa ) 


Pa/a 
sin  pja 


1/2 

9 


(4- 13a) 


d 


~  Pv(-  cos  0)  £  ka  J JkpJ) 


Pa'a 


sin  pja 


1/2 


and 


(4- 13b) 


dj_ 

dd2 


Pv(—  cos  0)  s  -  (ka)2 


J°(kpa)  ~  Jl(kpa) 


PcJa 
sin  pja 


1/2 


(4- 13c) 


Through  use  of  the  above  approximations,  antipode-centered  formulas  for  both  VED  and  HED 
sources  are  derived  from  the  spherical  waveguide  formulas  presented  in  section  4.2.  The  resulting 
HED  formulas  are  then  compared  with  the  corresponding  formulas  derived  in  reference  7. 

4.3.1  VED 


The  application  of  the  approximate  formulas  (4-13)  and  (1-7)  to  the  spherical  waveguide 
formulas  for  a  VED  source  given  in  expressions  (4-7)  yields  the  following  field  approximations: 


ive  jVokZP 
r  "  W0h 


1 

sin  v0n 


J()(kPa) 


pJA 

sin  pja 


1/2 

9 


(4- 14a) 
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1/2 


(4- 14b) 
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ve 
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VokPAg 

4  h 


1 

sin  vQn 
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sin  pa/« 
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V<? 

V 


¥  1 

4/2  sin 


J  l(kP  a) 


PjH_ 

sin  pa/nr 


1/2 

5 


(4- 14c) 


where  pa  =  a(n-  6)  is  the  great-circle  path  distance  from  the  antipode  to  the  field  point  (figure 
3-2).  In  the  above  formulas,  note  that  sin  6  =  sin  (p/a )  =  sin  (pa  /a)  and  p  +  pa  =  7ta. 

4.3.2  Horizontal  Electric  Dipole 


The  application  of  the  approximate  formulas  (4-13)  and  (1-7)  to  the  spherical  waveguide 
formulas  for  a  HED  source  given  in  expressions  (4-12)  results  in  the  following  field 
approximations: 


vMgP _ 

4 h  sin  v0k 


^(kpa) 


Pa !a 
sin  pja 


1/2 

COS  (p  , 


(4- 15  a) 


rphe 
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(4- 15b) 


(f 


4 khpa  sin  voji 


Pj° 


sin  pja 


3/2 


sin  (p 


4 khpa  sin  v0ji 


Pja 


sin  pja 


3/2 


sin  q> 
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Hke  K\P  I 

v  4 h  sin  v  it 


J°(kpa)  -  ij  Jl(kpa) 


Pala 


sin  pja 


1/2 


cos  q? 


(4- 15c) 


(4-15d) 


(4-15e) 


A  comparison  of  expressions  (4-15a),  (4-15d),  and  (4-15e)  with  the  corresponding  formulas 
given  in  reference  7  shows  that  the  above  formulas  differ  by  a  factor  C  that  is  defined  as 
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(4-16) 


eJkm  =  j  eJvQ* 

2  sin  v0n  2  sin  vja 

Note  that  the  final  expression  in  formula  (4-16)  was  obtained  through  use  of  equation  (4-8). 

Table  4-3  lists  the  magnitude  and  phase  of  C  for  several  frequencies  in  the  ELF  band.  The 
magnitude  of  C  is  tabulated  to  three  significant  figures  while  the  phase  of  C  is  given  to  the 
nearest  thousandth  of  a  degree.  Because  the  antipode-centered  propagation  formulas  in  reference 
7  were  derived  via  Bannister’s  HED  formulas  and  because  there  is  a  sign  difference  between 
Bannister’s  HED  formulas  and  the  spherical  waveguide  formulas  (see  section  4.4),  C  has  been 
multiplied  by  -1  in  table  4-3.  The  results  show  that  C  is  very  small  in  magnitude  and  phase  over 
the  ELF  band,  indicating  that  the  HED  antipode-centered  propagation  formulas  of  reference  7  are 
sufficiently  accurate  to  be  used  in  place  of  expressions  (4- 15a),  (4-15d),  and  (4-15e). 

Table  4-3.  Magnitude  and  Phase  of  Difference  Factor  C  in  HED  Antipode-Centered  ELF 
Propagation  Formulas  at  Several  Frequencies  under  Daytime  Conditions 


v0 

|C|(dB) 

Phase{-C}  (deg) 

1 

-  2.56  x  10"2 

3.599 

76 

1 

- 1.52  x  10-2 

0.204 

100 

15.9-7  L32 

1.92  x  10-3 

0.007 

300 

46.4-  /  3.67 

-  5.73  x  10"10 

0.000 

The  suitability  of  the  antipode-centered  ELF  propagation  listed  in  formulas  (4-14)  and  (4-15) 
depends  entirely  on  the  accuracy  of  the  approximations  for  Pv(-  cos  8)  and  its  first  two 

derivatives  with  respect  to  6  as  given  in  expressions  (4-13).  Comparison  plots  of  the  exact  series 
formulas  for  Pv(-  cos  8)  (and  its  first  two  derivatives  with  respect  to  8)  and  the  approximate 

formulas  (4-13)  are  given  in  appendix  J,  section  J.3,  for  frequencies  of  30  Hz  and  76  Hz.  The 
plots  show  that  the  approximate  formulas  are  in  close  agreement  with  the  corresponding  exact 
formulas.  This  comparison  improves  with  increasing  frequency  and  degrades  with  increasing 
derivative  order.  (For  more  details  regarding  this  comparison,  refer  to  appendix  J,  section  J.3.) 

In  summary,  the  antipode-centered  propagation  formulas  (4-14)  and  (4-15)  are  suitable 
substitutes  for  the  spherical  waveguide  formulas  presented  in  section  4.2. 
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4.4  REDUCTION  TO  BANNISTER’S  FORMULAS  VIA 
THE  EARTH-FLATTENING  APPROXIMATION 

The  earth-flattening  approximation  to  the  Legendre  function  of  the  first  kind  is  derived  in 
appendix  K.  Earlier  in  this  report,  it  was  mentioned  that  this  approximation  has  been  incorporated 
into  Bannister’s  VED  and  HED  formulas  given  in  section  3.2.  Through  use  of  the  earth-flattening 
approximation,  it  is  shown  that  the  VED  and  HED  spherical  waveguide  formulas  can  reduce  to 
Bannister’s  formulas. 


4.4.1  VED 


From  expression  (K-33)  in  appendix  K,  the  Legendre  function  of  the  first  kind  and  degree  v 
is  approximated  as 


Pv{-  cos  6)  s  j  sin  vk  H^\kp) 


pi  a 

sin  (p  /  a) 


1/2 


(4-17) 


where  p  =  a6  is  the  great -circle  path  distance  along  the  earth  from  the  source  point  to  the  field 
point  (figure  3-1)  and  H{2)  denotes  the  Hankel  function  of  the  second  kind  with  order  0.  As 
discussed  in  section  3.1,  H^\kp)  corresponds  to  the  range  dependence  for  a  flat  earth  and  the 
square-root  term  is  the  correction  for  curvature.  The  spherical  waveguide  propagation  formulas 
for  a  VED  are  given  in  expressions  (4-7).  In  these  formulas,  the  range  dependencies  are 
expressed  in  terms  of  either  Pv(—  cos  6)  or  its  first  derivative.  From  the  approximation  (4-17), 
the  derivative  of  Pv(-  cos  6)  with  respect  to  0  is  given  by 
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where  the  ^-derivative  of  H(2\kp)  is  given  by 
Jq  HgKkp)  =  ka  H^'(kp)  =  -ka  H\2\kp)  . 


(4-19) 


In  the  approximation  (4- 1 8),  note  that  the  6  derivative  was  applied  only  to  the  Hankel  function 
term  because  it  varies  more  rapidly  than  the  curvature  correction  term. 


The  substitution  of  the  approximations  (4-10),  (4-17)  and  (4-18)  along  with  expression 
(4-11)  into  the  VED  formulas  (4-7)  yields 
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In  the  above  formulas,  note  that  the  spherical  components  r  and  flhave  been  replaced  by  the 
equivalent  cylindrical  components  z  and  p,  respectively.  The  above  formulas  are  generally  valid 
for  ranges  greater  than  three  ionospheric  reflection  heights,  i.e.,  p  >  3 h.  In  order  to  show  that  the 
above  VED  formulas  are  equivalent  to  those  of  Bannister  given  in  section  3.2.1,  the  functions  Gh 
and  Vh  must  be  replaced  by  their  approximate  forms.  For  p  >  3h,  these  functions  can  be 
approximated  as 

G*(«>  -  #  =  f  ,  (4-21) 

and 


Vh(t)m  0. 


(4-22) 


The  substitutions  of  the  above  approximations  into  Bannister’s  VED  formulas  (3-3)  yields 
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A  comparison  of  formulas  (4-20)  with  (4-23)  shows  that  the  VED  spherical  waveguide  formulas 
reduce  to  Bannister’s  formulas  for  p>3h. 
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4.4.2  Horizontal  Electric  Dipole 


The  spherical  waveguide  propagation  formulas  for  a  HED  are  given  in  expressions  (4-12).  In 
these  formulas,  the  range  dependencies  are  expressed  in  terms  of  either  the  first  or  second 
derivatives  of  Pv(-  cos  0).  An  approximation  for  the  first  derivative  of  Pv(-  cos  0)  with  respect  to 
6  is  given  in  formula  (4- 1 8).  From  this  approximation,  the  second  derivative  of  Pv(-  cos  6)  is 
given  as 
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(4-24) 


In  the  above  approximation,  note  that  the  0  derivative  was  only  applied  to  the  Hankel  function 
term  because  it  varies  more  rapidly  than  the  curvature  correction  term. 


The  substitution  of  the  approximations  (4-10),  (4-17),  (4-18),  and  (4-24),  along  with 
expression  (4-1 1),  into  the  HED  formulas  (4-12)  yields 
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In  the  above  formulas,  note  that  the  spherical  components  r  and  0have  again  been  replaced  by 
the  equivalent  cylindrical  components  z  and  p,  respectively.  The  above  formulas  are  generally 

valid  for  p  >  3/z.  In  order  to  show  that  the  above  HED  formulas  are  equivalent  to  those  of 
Bannister  given  in  section  3.2.2,  the  functions  Gh,  Vh,  and  Hh  must  be  replaced  by  their 

approximate  forms.  Approximations  for  Gh(u)  and  Vh(t)  are  given  in  expressions  (4-21)  and  (4- 
22),  respectively.  For  p  >  3 h,  Gh(t)  and  Hh(t)  can  be  approximated  as  follows: 
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The  substitutions  of  the  approximations  (4-21),  (4-26),  and  (4-27)  into  Bannister’s  HED 
formulas  (3-8)  yields 
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A  comparison  of  formulas  (4-25)  with  (4-28)  indicates  that  each  of  the  HED  spherical  waveguide 
formulas  reduces  to  within  a  factor  of -1  of  Bannister’s  formulas  for  p  >  3 h.  At  this  time,  there 
is  no  explanation  for  the  sign  difference. 
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4.5  MODIFICATION  OF  HED  SPHERICAL  WAVEGUIDE  FORMULAS 
TO  ACCOUNT  FOR  AN  ANISOTROPIC  GROUND 


The  HED  spherical  waveguide  formulas  (4-12)  are  based  on  the  assumption  of  a 
homogeneous,  isotropic  earth.  These  formulas  can  be  modified  to  account  for  an  anisotropic- 
ground  conductivity  in  the  vicinity  of  the  HED  through  use  of  the  reciprocity  theorem.  The 
derivations  of  the  modified  HED  spherical  waveguide  formulas  follow  that  given  by  the  author 

(reference  6)  for  Bannister’s  HED  formulas  and  will  not  be  presented  here.  In  addition,  the 
conversion  (table  3-1)  may  be  applied  to  the  spherical  waveguide  formulas  if  Ez  and  Hp  are 

replaced  by  Er  and  Hq,  respectively.  Therefore,  the  modified  HED  spherical  waveguide 

formulas  are 
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In  the  above  formulas,  Te  and  Tn  are  the  antenna  pattern  factors  in  the  east  and  north  directions, 
respectively,  and  (pn  is  defined  in  figure  3-3.  As  with  Wolkoff  and  Kraimer’s  formulas  given  in 
section  3.3,  only  the  field  components  of  practical  interest  are  listed  above. 
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5.  COMPARISONS  OF  PROPAGATION  FORMULAS 


In  section  3,  Bannister’s  ELF  propagation  formulas  were  presented  for  both  VED  and  HED 
sources  radiating  in  a  spherical  earth-ionosphere  waveguide,  where  both  the  source  and  field 
points  are  located  on  the  earth’s  surface.  In  section  4,  spherical  waveguide  formulas  for  VED 
and  HED  sources  were  given.  The  spherical  waveguide  formulas  are  expected  to  be  more 
accurate  at  antipodal  ranges  because  they  have  a  range  dependence  that  is  more  accurate  than 
Bannister’s  formulas,  which  are  based  on  the  earth-flattening  approximation  with  curvature 
correction. 

In  this  section,  plots  of  the  spherical  waveguide  and  Bannister’s  ELF  propagation  formulas 
are  presented  as  functions  of  range  for  both  daytime  and  nighttime  propagation  conditions  at 
several  frequencies  in  the  ELF  band.  Results  are  presented  for  VED  and  HED  sources,  where 
both  the  source  and  field  points  are  located  on  the  surface  of  the  earth.  In  the  plots  of  Bannister’s 
formulas,  separate  graphs  of  the  direct  great-circle  path  field  and  the  summation  of  the  direct  and 
indirect  great-circle  path  fields  are  given.  The  comparisons  presented  here  will  help  to  more 
accurately  establish  the  maximum  ranges  of  validity  of  Bannister’s  formulas.  Plots  of  only  the 
dominant  field  components  are  presented,  i.e.,  Evze ,  ,  Eze ,  Hpe ,  and  H^e . 

For  a  HED  source,  the  plots  presented  in  this  report  are  at  the  azimuth  angle  (p  in  which  the 
magnitude  of  the  field  component  is  a  maximum.  In  particular,  the  plots  of  E1?6  and  are 
given  for  (p  =  0°  and  the  plots  of  Hpe  are  given  for  <p  =  90°.  In  addition,  for  the  spherical 
waveguide  formulas,  the  spherical  coordinates  r  and  6  are  replaced  by  the  cylindrical 
coordinates  z  and  p,  respectively,  where  z  =  r-a,p  =  ad,  and  a  is  the  earth’s  mean  radius.  The 
unit  vectors  in  these  coordinate  systems  are  related  as  follows:  z  =  r  and  p  =  6 . 

Table  5-1  is  a  list  of  the  propagation  parameters  at  several  ELF  frequencies  under  typical 
daytime  and  nighttime  conditions.  The  propagation  parameters  c/v,  a,  and  h  were  obtained  from 
Bannister  (reference  16)  and  are  based  on  an  exponential  conductivity  profile  of  the  ionosphere 
(reference  17).  The  degrees  of  the  Legendre  function  v0  in  the  spherical  waveguide  formulas 
that  are  listed  in  this  table  are  computed  from  formula  (4-8).  It  should  be  mentioned  that  the  76- 
Hz  parameters  listed  in  this  table  are  different  from  those  obtained  from  Wolkoff  (reference  18) 
for  use  in  an  earlier  investigation  (reference  7).  In  that  study,  Wolkoff  determined  the  76-Hz 
propagation  parameters  from  the  measured  data  listed  by  Bannister  (reference  17).  The  76-Hz 
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frequency  is  of  particular  interest  in  applications  because  it  is  the  center  frequency  of  the  U.S. 
Navy  Submarine  ELF  Communications  System. 


Table  5-1.  Propagation  Parameters  for  Typical  Daytime  and  Nighttime  Conditions  at 
Several  Frequencies  as  Obtained  from  Bannister  (Reference  16) 


Frequency 

(Hz) 

Propagation 

Condition 

c/v 

a  (dB/Mm) 

h  (km) 

V0 

30 

Daytime 

1.31 

0.60 

51 

4.75  -7  0.440 

30 

Nighttime 

1.14 

0.70 

73 

4.07-/  0.513 

76 

Daytime 

1.25 

1.4 

53.5 

12.2-/  1.03 

76 

Nighttime 

1.12 

0.90 

77 

10.9-7  0.660 

100 

Daytime 

1.23 

1.8 

55 

15.9-/  1.32 

100 

Nighttime 

1.12 

1.15 

77 

14.5-7  0.844 

300 

Daytime 

1.17 

5.0 

59 

46.4  -  /  3.67 

300 

Nighttime 

1.10 

2.7 

81 

43.6-/  1.98 

5.1  RESULTS  AT  76  HZ 

Plots  of  the  computed  magnitude  and  phase  of  the  vertical  electric  field  produced  by  a  VED 
under  daytime  propagation  conditions  at  76  Hz  are  given  in  figures  5-la  and  5-lb,  respectively.  In 
addition,  plots  of  the  magnitude  and  phase  of  the  azimuthal  magnetic  field  from  a  VED  source  are 
presented  for  the  same  propagation  conditions  in  figures  5-2a  and  5-2b,  respectively.  Each  of 
these  graphs  are  given  as  a  function  of  the  distance  measured  from  the  antipode  pa,  which  is 
defined  in  terms  of  the  direct  great-circle  path  distance  p  measured  from  the  source  as  pa  =  na  -  p 
and  illustrated  in  figure  3-2.  Figures  5- la  and  5 -2a  show  an  interference  pattern  that  is  greatest  in 
the  vicinity  of  the  antipode  and  slowly  diminishes  with  increasing  distance  from  the  antipode. 

The  interference  is  caused  by  the  combination  of  the  direct  and  indirect  great -circle  path  fields. 

In  figures  5-la  and  5-2a,  note  that  Bannister’s  direct  great-circle  path  field  components  do  not 
show  an  interference  pattern  because  they  do  not  include  the  indirect  great-circle  path 
contribution. 
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In  figure  5-la,  Bannister’s  total  field  (direct  plus  indirect  great-circle  path  formulas)  result 
for  Evz  shows  good  agreement  with  the  spherical  waveguide  formula  for  all  ranges  except  those 
that  are  very  close  to  the  antipode.  In  particular,  Bannister’s  total  field  result  for  Evze  agrees  to 
within  1  dB  of  the  spherical  waveguide  formula  for  pa  -  1.15  Mm  or,  equivalently,  p  =  18.87 
Mm.  In  comparison,  Bannister’s  direct  great-circle  path  formula  for  E\e  agrees  to  within  1  dB 
of  the  spherical  waveguide  formula  for  pa  =  6.09  Mm  or,  equivalently,  p  =  13.93  Mm.  Clearly, 
under  these  propagation  conditions,  Bannister’s  total  field  results  are  much  more  accurate  than 
Bannister’s  direct  great-circle  path  results  at  antipodal  ranges. 

A  comparison  of  figures  5-la  and  5-2a  indicates  that  Bannister’s  total  field  result  for 
does  not  agree  as  closely  with  the  corresponding  spherical  waveguide  formula  as  Evze .  In  figure 
5-2a,  Bannister’s  total  field  result  for  Hv*  agrees  to  within  1  dB  of  the  spherical  waveguide 
formula  for  pa  =  1.89  Mm  or  equivalently,  p  =  18.13  Mm.  In  addition,  Bannister’s  direct  great- 
circle  path  formula  for  Hv*  agrees  to  within  1  dB  of  the  spherical  waveguide  formula  for 
pa  =  6.78  Mm  or  equivalently,  p  =  13.24  Mm.  The  poorer  agreement  in  Bannister’s  total  field 
result  for  H™  is  because  its  range  dependence  is  proportional  to  the  first  derivative  of 
Pv(-  cos  6)  as  compared  with  Ev’e ,  which  has  a  range  dependence  that  is  directly  proportional  to 
Pv(-  cos  6)  •  This  observation  is  attributed  to  the  fact  that  in  the  earth-flattening  approximation 
for  p  (-  cos  6) ,  each  successive  derivative  introduces  additional  error  into  the  approximation. 

The  phase  plots  for  Evze  and  Hv‘ ,  given  in  figures  5-lb  and  5-2b,  respectively,  show  a 
similar  comparison  of  Bannister’s  results  with  the  spherical  waveguide  formulas  as  was  observed 
for  the  field  magnitudes.  In  figure  5- lb,  Bannister’s  total  field  result  for  E'’e  agrees  to  within  5° 
in  phase  of  the  spherical  waveguide  formula  result  for  pa  =  1.35  Mm  or,  equivalently,  p  =  18.67 
Mm.  In  figure  5-2b,  Bannister’s  total  field  result  for  H™  agrees  to  within  5°  in  phase  of  the 
spherical  waveguide  formula  result  for  pa  =  3.67  Mm  or,  equivalently,  p  =  16.35  Mm. 

Bannister’s  direct  great-circle  path  formulas  for  Evze  and  H™  agree  to  within  5°  in  phase  of  the 
corresponding  spherical  waveguide  formula  results  for  pa  =  7.29  Mm  and  pa  =  7.30  Mm, 
respectively.  As  expected,  the  phase  results  computed  from  Bannister’s  total  field  formulas 
agree  more  closely  with  the  spherical  waveguide  formulas  than  do  Bannister’s  direct  great-circle 
path  formulas. 
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Plots  of  the  computed  magnitudes  and  phases  of  the  HED  field  components  E^e ,  H^e , 
and  H'y  under  daytime  propagation  conditions  at  76  Hz  are  given  in  figures  5-3,  5-4,  and  5-5, 
respectively.  The  magnitude  plots  show  that  Bannister’s  total  field  results  for  EE ,  Hhe ,  and 

^  r 

H ^  agree  to  within  1  dB  of  the  corresponding  spherical  waveguide  formula  results  for  pa=  1.89 
Mm,  pa-  2.12  Mm,  and  pa  =  2.67  Mm,  respectively.  Note  that  the  ranges  of  agreement  for  H^e 
and  FJ'e  are  the  same.  This  observation  is  not  surprising  because  the  range  dependencies  of 
these  field  components  are  each  proportional  to  the  first  derivative  of  Pv(-  cos  6)  •  The  poorer 
agreement  for  H^e  is  because  its  range  dependence  is  proportional  to  the  second  derivative  of 
j P  (-  cos  6) ,  resulting  in  additional  error  in  the  earth-flattening  approximation. 

In  figures  5-3a,  5-4a,  and  5-5a,  Bannister’s  direct  great -circle  path  results  for  EE ,  H^e ,  and 
H^e  agree  to  within  1  dB  of  the  corresponding  spherical  waveguide  formula  results  for  pa  =  6.78 
Mm,  pa  =  6.77  Mm,  and  pa  =  6.08  Mm,  respectively.  Note  that  the  ranges  of  agreement  for  Hvqe , 
EE ,  and  Hhne  are  nearly  identical  because  the  formulas  have  the  same  range  dependence. 
Surprisingly,  H^e  shows  a  better  agreement  with  the  spherical  waveguide  formulas  than  the 
other  direct  great-circle  path  field  components. 

The  phase  plots  for  EE ,  HE ,  and  ,  given  in  figures  5-3b,  5-4b,  and  5-5b,  respectively, 
show  a  similar  comparison  of  Bannister’s  results  with  the  spherical  waveguide  formulas  as  was 
observed  for  the  field  magnitudes.  In  these  plots,  Bannister’s  total  field  results  for  EJ}C ,  ,  and 

Hy  agree  to  within  5°  in  phase  of  the  corresponding  spherical  waveguide  formula  results  for 
pa  =  3.67  Mm,  pa  =  3.67  Mm,  and  pa  =  6. 14  Mm,  respectively.  Based  on  these  comparisons,  the 
phase  agreement  for  each  of  these  field  components  is  worse  than  observed  in  the  corresponding 
magnitude  comparisons.  Bannister’s  direct  great-circle  path  formulas  for  EE ,  H1™ ,  and  H^e 
agree  to  within  5°  in  phase  of  the  corresponding  spherical  waveguide  formula  results  for  pa  = 

7.30  Mm,  pa  =  7.3 1  Mm,  and  pa  =  9.56  Mm,  respectively.  As  was  observed  for  the  VED  fields, 
the  phase  results  computed  from  Bannister’s  total  field  HED  formulas  agree  more  closely  with 
the  spherical  waveguide  formulas  than  do  Bannister’s  direct  great-circle  path  formulas. 


-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


(The  propagation  parameters  are  c/v  =  1.25,  a  =  1.4  dB/Mm,  h  =  53.5  km,  and  the  dipole  moment  is  p  =  1  Am.) 


Figure  5-lb.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Phase  of  the  Vertical 
Electric  Field  Produced  by  a  VED  Under  Typical  Daytime  Propagation  Conditions  at  76  Hz 

(The  propagation  parameters  are  c/v  =  1.25,  a=  1.4  dB/Mm,  h  =  53.5  km,  and  the  dipole  moment  isp  =  1  Am.) 
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Figure  5-2a.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the 
Azimuthal  Magnetic  Field  Produced  by  a  VED  Under  Typical  Daytime  Propagation  Conditions  at  76  Hz 
(The  propagation  parameters  are  civ  =  1.25,  a-  1.4  dB/Mm,  h  =  53.5  km,  and  the  dipole  moment  is  p  =  1  Am.) 
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Figure  5-2b.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Phase  of  the  Azimuthal 
Magnetic  Field  Produced  by  a  VED  Under  Typical  Daytime  Propagation  Conditions  at  76  Hz 

(The  propagation  parameters  are  c/v  =  1.25,  a=  1.4  dB/Mm,  h  =  53.5  km,  and  the  dipole  moment  is  p  =  1  Am.) 


Electric  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  76  Hz 

(The  propagation  parameters  are  c/v  =  1.25,  a=  1.4  dB/Mm,  h  =  53.5  km;  the  azimuthal  angle  is  <p=  0°;  and  the  dipole  moment  isp  =  1  Am.) 


Figure  5-3b.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Phase  of  the  Vertical 
Electric  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  76  Hz 

(The  propagation  parameters  are  c/v  =  1.25,  a=  1.4  dB/Mm,  h  =  53.5  km;  the  azimuthal  angle  is  <p=  0°;  and  the  dipole  moment  is  p  =  1  Am.) 
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Magnetic  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  76  Hz 
(The  propagation  parameters  are  c/v  =  1.25,  a=  1.4  dB/Mm,  h  =  53.5  km;  the  azimuthal  angle  is  <p=  90°;  and  the  dipole  moment  is  p  -  1  Am.) 
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Figure  5-4b.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Phase  of  the  Radial 
Magnetic  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  76  Hz 

(The  propagation  parameters  are  c/v  =  1.25,  a=  1.4  dB/Mm,  h  =  53.5  km;  the  azimuthal  angle  is  <p=  90°;  and  the  dipole  moment  is  p  =  1  Am.) 


Azimuthal  Magnetic  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  76  Hz 

(The  propagation  parameters  are  civ  =  1.25,  a=  1.4  dB/Mm,  h  =  53.5  km;  the  azimuthal  angle  is  <p=  0°;  and  the  dipole  moment  is  p  =  1  Am.) 


Figure  5-5b.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Phase  of  the  Azimuthal 
Magnetic  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  76  Hz 

(The  propagation  parameters  are  c/v  =  1.25,  a-  1.4  dB/Mm,  h  =  53.5  km;  the  azimuthal  angle  is  <p= 0°;  and  the  dipole  moment  is  p  =  1  Am.) 


Plots  of  the  computed  magnitudes  of  the  VED  field  components  Eve  and  Hvne  under  typical 

4  ly 

nighttime  propagation  conditions  at  76  Hz  are  given  in  figures  5-6a  and  5-6b,  respectively. 
Similar  plots  of  the  computed  magnitudes  of  the  HED  field  components  E*e ,  H1™ ,  and  H*E 

**  r  t' 

under  typical  nighttime  propagation  conditions  at  76  Hz  are  given  in  figures  5-7 a,  5-7b,  and  5-7c, 
respectively.  In  these  plots,  Bannister’s  total  field  results  for  Eve ,  HVJ ,  E,e ,  H^e ,  and 
agree  to  within  1  dB  of  the  corresponding  spherical  waveguide  formula  results  for  pa  >  1.49  Mm, 
pa  >  2.50  Mm,  pa  >  2.50  Mm,  pa  >  2.5 1  Mm,  and  pa  >  3.47  Mm,  respectively.  Although  these 
agreements  are  not  quite  as  good  as  the  corresponding  daytime  propagation  results,  they  follow 
the  same  trend.  In  comparison,  Bannister’s  direct  great-circle  path  field  results  for  Evze ,  H™ , 
E!E ,  Hhne ,  and  H’E  agree  to  within  1  dB  of  the  corresponding  spherical  waveguide  formula 
results  for  pa  >  10.20  Mm,  pa  >  10.24  Mm,  pa  >  10.24  Mm,  pa  >  10.25  Mm,  and  pa  >  10.21  Mm, 
respectively.  These  ranges  of  agreement  are  considerably  worse  than  what  was  observed  for  the 
direct  great-circle  path  field  results  under  daytime  conditions.  This  observation  is  attributed  to 
the  lower  attenuation  at  nighttime  which  results  in  an  interference  pattern  from  the  two  great- 
circle  path  fields  that  extends  to  greater  distances  from  the  antipode. 

A  summary  listing  of  the  ranges  over  which  Bannister’s  direct  great-circle  path  and  total 
field  formulas  agree  to  within  1  dB  in  magnitude  of  the  corresponding  spherical  waveguide 
formulas  are  given  in  tables  5-2a  and  5-2b,  respectively,  for  both  typical  daytime  and  nighttime 
propagation  conditions  at  76  Hz.  In  these  tables,  the  normalized  range  fipa  is  included  in 
parentheses,  where  ft  =  Re{fc}  =  2nlA  is  the  phase  constant  of  the  quasi-TEM  wave  and  /t  is  the 
corresponding  wavelength.  (Normalized  ranges  will  be  discussed  later  in  this  section.) 

A  comparison  of  these  tables  clearly  shows  that  Bannister’s  total  field  formulas  have  a  much 
greater  range  of  agreement  with  the  spherical  waveguide  formulas  than  do  Bannister’s  direct 
great-circle  path  formulas.  From  table  5-2a,  the  ranges  of  satisfactory  agreement  for  Bannister’s 
direct  great-circle  path  formulas  extend  from  pa  >  6.08  Mm  to  pa  >  10.25  Mm,  depending  on  the 
field  component  and  propagation  condition.  In  comparison,  from  table  5-2b,  the  ranges  of 
satisfactory  agreement  for  Bannister’s  total  field  formulas  extend  from  pa>  1.15  Mm  to  pa  > 

3.47  Mm.  As  previously  mentioned,  the  poorer  agreement  observed  for  Bannister’s  direct  great- 
circle  path  formulas  is  attributed  to  their  inability  to  predict  the  interference  produced  by  the 
superposition  of  the  direct  and  indirect  great-circle  path  fields. 
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Table  5-2a.  Ranges  Over  Which  Bannister’s  Direct  Great-Circle  Path  Field  Formulas  Agree 
to  Within  1  dB  in  Magnitude  of  the  Spherical  Waveguide  Formulas  at  76  Hz 


Field 

Component 

Propagation 

Condition 

Range 

Ele 

4 

Daytime 

pa > 6.09 Mm  (J3pa>  12.13) 

Nighttime 

pa  >  10.20  Mm  0 Qpa  >  18.20) 

H;e 

Daytime 

pa  >  6.78  Mm  (J3pa  >  13.50) 

Nighttime 

pa>  10.24 Mm  (J3pa>  18.27) 

E t 

Daytime 

pa  >  6.78  Mm  (J3pa  >  13.50) 

Nighttime 

pa  >  10.24  Mm  (fipa  >  18.27) 

h* 

Daytime 

Pa  ^  6.77  Mm  (#>a>  13.48) 

Nighttime 

/?a  >  10.25  Mm  0 &Pa  >  18-29) 

jet  he 

“  (p 

Daytime 

Pa  >  6.08  Mm  (J3pa>  12.11) 

Nighttime 

>  10.21  Mm  >  18.21) 

Table  5-2b.  Ranges  Over  Which  Bannister’s  Total  Field  Formulas  Agree  to  Within 
1  dB  in  Magnitude  of  the  Spherical  Waveguide  Formulas  at  76  Hz 


Field 

Component 

Propagation 

Condition 

Range 

Daytime 

pa  >  1.15  Mm  (J3pa>  2.29) 

Nighttime 

pa  >  1.49  Mm  (J3pa>  2.66) 

H;e 

Daytime 

/?a >  1.89 Mm  (J3pa>2.16) 

Nighttime 

pa>  2.50  Mm  (#?a>4.46) 

E*e 

4 

Daytime 

/?a  >  1.89  Mm  {fipa  >  3.76) 

Nighttime 

>  2.50  Mm  (/?/?a  >  4.46) 

jet  he 

Daytime 

pa  >  2.12  Mm  (/?/?a  >  4.22) 

Nighttime 

pa  >  2.51  Mm  (/?pa  >  4.48) 

jet  he 
**  cp 

Daytime 

pa > 2.67  Mm  (J3pa>532) 

Nighttime 

/?a>3.47Mm  (/?/9a  >  6.19) 
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-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Electric  Field  Produced  by  a  VED  Under  Typical  Nighttime  Propagation  Conditions  at  76  Hz 
(The  propagation  parameters  are  cfv  =  1.12,  a =  0.9  dB/Mm,  h-11  km,  and  the  dipole  moment  is  p  =  1  Am.) 


Figure  5-6b.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the 
Azimuthal  Magnetic  Field  Produced  by  a  VED  Under  Typical  Nighttime  Propagation  Conditions  at  76  Hz 
(The  propagation  parameters  are  civ  =1.12,  a=  0.9  dB/Mm,  h  =  77  km,  and  the  dipole  moment  isp  =  1  Am.) 
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Electric  Field  Produced  by  a  HED  Under  Typical  Nighttime  Propagation  Conditions  at  76  Hz 
(The  propagation  parameters  are  cfv  =  1.12,  a=  0.9  dB/Mm,  h-11  km;  the  azimuthal  angle  is  <p=  0°;  and  the  dipole  moment  is  p  =  1  Am.) 
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Figure  5-7b.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the  Radial 
Magnetic  Field  Produced  by  a  HED  Under  Typical  Nighttime  Propagation  Conditions  at  76  Hz 

(The  propagation  parameters  are  civ  =  1.12,  a=  0.9  dB/Mm,  h  =  77  km;  the  azimuthal  angle  is  (p=  90°;  and  the  dipole  moment  is  p  =  1  Am.) 
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Azimuthal  Magnetic  Field  Produced  by  a  HED  Under  Typical  Nighttime  Propagation  Conditions  at  76  Hz 

(The  propagation  parameters  are  c/v  =  1.12,  a=  0.9  dB/Mm,  h  =  77  km;  the  azimuthal  angle  is  <p~  0°;  and  the  dipole  moment  is  p  =  1  Am.) 


5.2  RESULTS  AT  OTHER  FREQUENCIES 


In  order  to  obtain  a  better  assessment  of  the  maximum  ranges  of  validity  of  Bannister’s 
propagation  formulas,  it  is  necessary  to  compare  them  with  the  spherical  waveguide  formulas  at 
additional  frequencies  across  the  ELF  band.  Plots  of  the  computed  magnitudes  and  phases  of  the 
VED  field  components  Ev_e  and  H™  under  typical  daytime  propagation  conditions  at  30  Hz  are 
given  in  figures  5-8  and  5-9,  respectively.  Similar  plots  of  the  HED  field  components  E^e ,  Hpe , 
and  Hy  under  typical  daytime  propagation  conditions  at  30  Hz  are  given  in  figures  5-10,  5-11, 
and  5-12,  respectively.  A  comparison  of  these  plots  with  the  corresponding  76-Hz  graphs 
indicates  that  the  30-Hz  field  values  are  generally  larger  in  magnitude,  and  the  interference 
produced  by  the  two  great-circle  path  fields  at  30  Hz  extends  to  much  greater  distances  from  the 
antipode.  These  observations  are  the  result  of  the  lower  attenuation  at  30  Hz,  which  allows  the 
indirect  great-circle  path  fields  to  propagate  over  greater  ranges  from  the  source.  At  30  Hz, 
Bannister’s  total  field  results  for  Ev_e ,  Hv‘ ,  E^e ,  Hhpe ,  and  H^e  agree  to  within  1  dB  of  the 
corresponding  spherical  waveguide  formula  results  for  pa  >  2.71  Mm,  pa  >  4.55  Mm,  pa  >  4.55 
Mm,  pa  >  4.52  Mm,  and  pa  >  6.44  Mm,  respectively.  These  results  indicate  that  Bannister’s  total 
field  formulas  show  a  poorer  agreement  with  the  spherical  waveguide  formulas  than  the  76-Hz 
results.  Similarly,  Bannister’s  direct  great-circle  path  formulas  do  not  agree  as  well  with  the 
spherical  waveguide  formulas  as  at  76  Hz. 

Magnitude  plots  of  the  HED  field  components  E^e ,  Hpe ,  and  under  typical  nighttime 
propagation  conditions  at  30  Hz  are  given  in  figures  5-13a,  5-13b,  and  5-13c,  respectively.  In 
these  plots,  Bannister’s  total  field  formulas  for  E*}e ,  Hhne ,  and  H^e  agree  to  within  1  dB  of  the 
corresponding  spherical  waveguide  formula  results  for  pa  >  2.76  Mm,  pa  >  2.83  Mm,  and  pa  > 
4.69  Mm,  respectively.  Unlike  the  76-Hz  results,  the  30-Hz  data  show  that  Bannister’s  total  field 
formulas  produce  a  better  agreement  with  the  spherical  waveguide  formulas  at  nighttime  than  in 
the  daytime.  The  one  exception  to  this  observation  is  Ev,e ,  which  shows  a  worse  agreement  at 
nighttime.  The  general  improvement  in  the  nighttime  results  at  30  Hz  is  probably  attributed  to 
the  larger  nighttime  attenuation  as  indicated  in  table  5-1. 

Table  5-3  is  a  summary  list  of  the  ranges  over  which  Bannister’s  total  field  formulas  agree 
to  within  1  dB  in  magnitude  of  the  corresponding  spherical  waveguide  formulas  for  both  daytime 
and  nighttime  propagation  conditions  at  30  Hz.  A  comparison  of  tables  5-2b  and  5-3  clearly 
shows  that  Bannister’s  total  field  results  improve  at  the  higher  frequency. 
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Table  5-3.  Ranges  Over  Which  Bannister’s  Total  Field  Formulas  Agree  to  Within 
1  dB  in  Magnitude  of  the  Spherical  Waveguide  Formulas  at  30  Hz 


Field 

Component 

Propagation 

Condition 

Range 

E? 

Daytime 

pa  >  2.71  Mm  {fipa  >  2.23) 

Nighttime 

pa  >  3.10  Mm  (fipa  >  2.22) 

Daytime 

pa  >  4.55  Mm  {fipa  >3.75) 

Nighttime 

pa  > 2.76  Mm  {fipa>  1.98) 

£he 

Daytime 

pa  >  4.55  Mm  0 6pa  >  3.75) 

Nighttime 

pa  ^  2.76  Mm  (J3pa  >  1.98) 

Daytime 

pa  >  4.52  Mm  (J3pa  >3.72) 

Nighttime 

pa  >  2.83  Mm  (J3pa  >  2.03) 

rjhe 

ti  y 

Daytime 

pa  >  6.44  Mm  (fipa>  5.30) 

Nighttime 

pa  >  4.69  Mm  (J3pa  >  3.36) 

Plots  of  the  computed  magnitudes  of  the  field  components  Eve ,  HVJ ,  Ehe ,  Hhe ,  and  Hhe 
under  typical  daytime  propagation  conditions  at  300  Hz  are  given  in  figures  5-14a,  5-14b,  5-15a, 
5-15b,  and  5-15c,  respectively.  The  plots  show  that  the  interference  pattern  produced  by  the 
superposition  of  the  two  great-circle  path  fields  attenuates  rapidly  with  distance  from  the 
antipode.  This  observation  is  the  result  of  the  large  value  of  attenuation  at  this  higher  frequency. 
The  plots  also  show  that  Bannister’s  formulas  produce  a  better  agreement  with  the  corresponding 
spherical  waveguide  formulas  than  at  the  lower  frequencies.  In  particular,  Bannister’s  total  field 
results  for  Evze ,  H™ ,  Eze ,  Hpe ,  and  H^e  agree  to  within  1  dB  of  the  corresponding  spherical 
waveguide  formula  results  for  pa  >  0.32  Mm,  pa  >  0.59  Mm,  pa  >  0.59  Mm,  pa  >  0.59  Mm,  and 
pa  >  0.82  Mm,  respectively.  Bannister’s  direct  great-circle  path  formulas  also  produce  a  better 
agreement  with  the  spherical  waveguide  formulas  at  300  Hz  than  at  30  Hz  and  76  Hz  because  the 
interference  pattern  does  not  extend  to  as  great  of  a  distance  from  the  antipode  as  at  the  lower 
frequencies. 
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-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Figure  5-8a.  Comparison  of  Spherical  Waveguide  and  Bannister's  ELF  Propagation  Formulas  for  the  Magnitude  of  the  Vertical 
Electric  Field  Produced  by  a  VED  Under  Typical  Daytime  Propagation  Conditions  at  30  Hz 

(The  propagation  parameters  are  civ  =  1.31,  a=  0.6  dB/Mm,  h  =  51  km,  and  the  dipole  moment  is  p  =  1  Am.) 


Spherical  waveguide 
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Bannister  (dir.+ind.) 


Electric  Field  Produced  by  a  VED  Under  Typical  Daytime  Propagation  Conditions  at  30  Hz 
(The  propagation  parameters  are  civ  -  1.31,  a-  0.6  dB/Mm,  h-  51  km,  and  the  dipole  moment  is  p  =  1  Am.) 


Figure  5-9a.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the 
Azimuthal  Magnetic  Field  Produced  by  a  VED  Under  Typical  Daytime  Propagation  Conditions  at  30  Hz 

(The  propagation  parameters  are  c/v  =  1.31,  a=  0.6  dB/Mm,  h  =  51  km,  and  the  dipole  moment  is  p  -  1  Am.) 
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Magnetic  Field  Produced  by  a  VED  Under  Typical  Daytime  Propagation  Conditions  at  30  Hz 
(The  propagation  parameters  are  civ  =  1.31,  a=  0.6  dB/Mm,  h  =  51  km,  and  the  dipole  moment  is  p  =  1  Am.) 
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-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Figure  5-  10a.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the  Vertical 
Electric  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  30  Hz 

(The  propagation  parameters  are  c/v  =  1.31,  a=  0.6  dB/Mm,  h  =  51  km;  the  azimuthal  angle  is  (p=  0°;  and  the  dipole  moment  is  p  =  1  Am.) 
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Electric  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  30  Hz 
(The  propagation  parameters  are  c/v  =  1.31,  a=  0.6  dB/Mm,  h  =  51  km;,  the  azimuthal  angle  is  <p=  0°;  and  the  dipole  moment  is/?  =  1  Am.) 


Spherical  waveguide 
Bannister  (dir.) 


Figure  5-1  la.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the  Radial 
Magnetic  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  30  Hz 

(The  propagation  parameters  are  c/v  =  1 .3 1 ,  a  =  0.6  dB/Mm,  h  =  5 1  km;  the  azimuthal  angle  is  (p  =  90°;  and  the  dipole  moment  is  p  =  1  Am.) 
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Magnetic  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  30  Hz 

(The  propagation  parameters  are  civ  =  1.31,  a=  0.6  dB/Mm,  h  =  51  km;  the  azimuthal  angle  is  (p=  90°;  and  the  dipole  moment  is  p  =  1  Am.) 
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Magnetic  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  30  Hz 
(The  propagation  parameters  are  c/v  =  1.31,  a- 0.6  dB/Mm,  h  =  51  km;  the  azimuthal  angle  is  (p- 0°;  and  the  dipole  moment  isp  =  1  Am.) 


-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Figure  5-13a.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the  Vertical 
Electric  Field  Produced  by  a  HED  Under  Typical  Nighttime  Propagation  Conditions  at  30  Hz 

(The  propagation  parameters  are  c/v  =  1.14,  a=  0.7  dB/Mm,  h  =  73  km;  the  azimuthal  angle  is  <p=  0°;  and  the  dipole  moment  is  p  =  1  Am.) 
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-  Bannister  (dir.+ind.) 


Magnetic  Field  Produced  by  a  HED  Under  Typical  Nighttime  Propagation  Conditions  at  30  Hz 
(The  propagation  parameters  are  c/v  =  1.14,  a=  0.7  dB/Mm,  h  =  73  km;  the  azimuthal  angle  is  <p=  90°;  and  the  dipole  moment  isp  =  1  Am.) 


-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Figure  5- 13  c.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the 
Azimuthal  Magnetic  Field  Produced  by  a  HED  Under  Typical  Nighttime  Propagation  Conditions  at  30  Hz 

(The  propagation  parameters  are  c/v  =  1.14,  oc=  0.7  dB/Mm,  h  =  73  km;  the  azimuthal  angle  is  <p=  0°;  and  the  dipole  moment  is  p  =  1  Am.) 


-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Electric  Field  Produced  by  a  VED  Under  Typical  Daytime  Propagation  Conditions  at  300  Hz 
(The  propagation  parameters  are  c/v  =1.17,  a=  5.0  dB/Mm,  h  ~  59  km,  and  the  dipole  moment  is  /;  =  1  Am.) 


Figure  5-14b.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the 
Azimuthal  Magnetic  Field  Produced  by  a  VED  Under  Typical  Daytime  Propagation  Conditions  at  300  Hz 
(The  propagation  parameters  are  c/v  =1.17,  a=  5.0  dB/Mm,  h  =  59  km,  and  the  dipole  moment  is  p  =  1  Am.) 
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Figure  5- 15a.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the  Vertical 
Electric  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  300  Hz 
(The  propagation  parameters  are  c/v  =  1.17,  a=  5.0  dB/Mm,  h  =  59  km;  the  azimuthal  angle  is  (p- 0°;  and  the  dipole  moment  is  p  =  1  Am.) 


-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Figure  5-15b.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the  Radial 
Magnetic  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  300  Hz 

(The  propagation  parameters  are  c/v  =  1.17,  a=  5.0  dB/Mm,  h  =  59  km;  the  azimuthal  angle  is  <p=  90°;  and  the  dipole  moment  is  p  =  1  Am.) 
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Figure  5-15c.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the 
Azimuthal  Magnetic  Field  Produced  by  a  HED  Under  Typical  Daytime  Propagation  Conditions  at  300  Hz 
(The  propagation  parameters  are  civ  =  1.17,  a=  5.0  dB/Mm,  h  =  59  km;,  the  azimuthal  angle  is  <p= 0°;  and  the  dipole  moment  is  /;  =  1  Am.) 


Magnitude  plots  of  the  HED  field  components  Ef ,  Hf ,  and  H’f  under  typical  nighttime 
propagation  conditions  at  300  Hz  are  given  in  figures  5-16a,  5-16b,  and  5-16c,  respectively. 
Because  of  the  lower  attenuation  at  nighttime,  the  interference  pattern  produced  by  the  two  great- 
circle  path  fields  extends  to  greater  distances  from  the  antipode.  As  a  result,  at  nighttime, 
Bannister’s  total  field  results  at  300  Hz  do  not  compare  as  well  with  the  spherical  waveguide 
formulas  as  under  daytime  conditions.  A  similar  observation  was  made  at  76  Hz.  In  the  300-Hz 
plots,  Bannister’s  total  field  formulas  for  E^e ,  Hhpe ,  and  agree  to  within  1  dB  of  the 
corresponding  spherical  waveguide  formula  results  for  pa  >  0.65  Mm,  pa  >  0.65  Mm,  and  pa  > 
1.22  Mm,  respectively. 

Table  5-4  is  a  list  of  the  ranges  over  which  Bannister’s  total  field  formulas  agree  to  within 
1  dB  in  magnitude  of  the  corresponding  spherical  waveguide  formulas  for  both  typical  daytime 
and  nighttime  propagation  conditions  at  300  Hz.  A  comparison  of  tables  5-2b,  5-3,  and  5-4 
indicates  that  Bannister’s  total  field  results  improve  with  increasing  frequency. 

Table  5-4.  Ranges  Over  Which  Bannister's  Total  Field  Formulas  Agree  to  Within 
1  dB  in  Magnitude  of  the  Spherical  Waveguide  Formulas  at  300  Hz 


Field 

Component 

Propagation 

Condition 

Range 

Elf 

Daytime 

pa>  0.32  Mm  (J3pa  >  2.35) 

Nighttime 

pa  >  0.40  Mm  (J3pa  >  2.77) 

Daytime 

pa  >  0.59  Mm  (J3pa>4.34) 

Nighttime 

pa > 0.65  Mm  (J3pa> 4.50) 

E * 

Daytime 

pa  >  0.59  Mm  (J3pa  >  4.34) 

Nighttime 

pa  >  0.65  Mm  (J3pa  >  4.50) 

Daytime 

pa  >  0.59  Mm  (J3pa>  4.34) 

nighttime 

pa  >  0.65  Mm  (J3pa>  4.50) 

rrhe 
*1  cp 

Daytime 

pa  > 0.82  Mm  (J3pa> 6.03) 

Nighttime 

pa  >  1.22  Mm  (J3pa  >  8.44) 

Table  5-5  is  a  list  of  the  average  normalized  ranges  (expressed  in  terms  of  J3pa)  over  which 
Bannister’s  total  field  formulas  agree  to  within  1  dB  in  magnitude  of  the  corresponding  spherical 
waveguide  formulas.  The  ranges  have  been  normalized  in  order  to  scale  them  with  respect  to  a 
wavelength  in  the  waveguide.  The  ranges  are  averaged  over  the  three  frequencies  discussed  in 
this  report,  namely,  30  Hz,  76  Hz,  and  300  Hz.  The  third  column  in  the  table  lists  the  average 
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ranges  for  each  field  component  under  typical  daytime  and  nighttime  propagation  conditions. 

The  fourth  row  in  the  table  lists  the  ranges  averaged  over  daytime  and  nighttime  propagation 
conditions  for  each  field  component. 

The  third  column  of  table  5-5  indicates,  on  the  average,  that  Bannister’s  total  field  formulas 
for  Evze  and  H^e  produce  better  agreement  with  the  spherical  waveguide  formulas  at  daytime 
than  at  nighttime.  The  opposite  result  is  true  for  the  other  field  components,  HVJ ,  E^e ,  Hhne . 

t  r 

Note  that  the  better  agreement  in  the  field  components  at  nighttime  was  only  observed  at  30  Hz, 
where  the  nighttime  attenuation  is  greater  than  in  the  daytime.  The  results  in  the  fourth  column 
of  table  5-5  show  that  Bannister’s  total  field  formula  for  Elf  provides  the  best  agreement  with 
the  spherical  waveguide  formula  and  produces  the  worst  agreement.  As  mentioned,  this 
observation  is  attributed  to  the  fact  that  in  the  earth-flattening  approximation  for  Pv(-  cos  6) , 
each  successive  derivative  introduces  additional  error  into  the  approximation.  Also  note  in  the 
fourth  column  of  table  5-5  that  HVJ ,  ,  and  H^f  each  produce  a  similar  average  agreement 

with  the  corresponding  spherical  waveguide  formulas.  This  result  is  attributed  to  the  fact  that  the 
range  dependence  of  each  of  these  components  is  proportional  to  the  first  derivative  of 
Pv(-  cos  6) .  In  summary,  table  5-5  clearly  shows  that  the  spherical  waveguide  propagation 
formulas  are  necessary  for  prediction  of  the  ELF  fields  for  ranges  that  are  close  to  the  antipode. 


Table  5-5.  Average  Normalized  Ranges  Over  Which  Bannister’s  Total  Field  Formulas 
Agree  to  Within  1  dB  in  Magnitude  of  the  Spherical  Waveguide  Formulas 
(Averages  Are  Taken  Over  30,  76,  and  300  Hz) 


Field 

Component 

Propagation 

Condition 

Average 

Normalized  Range 

Avg  Norm.  Range 
(Day  and  Night) 

E\e 

Daytime 

Ppa>  2-29 

Ppa>  2-42 

Nighttime 

Ppa>  2-55 

Daytime 

Ppa>  3-95 

ppa>  3.80 

Nighttime 

Ppa  >  3.65 

Eze 

Daytime 

0Pa>  3-95 

ppa>  3.80 

Nighttime 

Ppa>  3-65 

H* 

Daytime 

Pp*>  4-09 

fipa  >  3.88 

Nighttime 

Ppa>  3-67 

t rhe 
tty 

! 

Daytime 

ppa>  5.55 

0Pa>  5-77 

Nighttime 

fra>  600 
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-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Figure  5-16a.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the  Vertical 
Electric  Field  Produced  by  a  HED  Under  Typical  Nighttime  Propagation  Conditions  at  300  Hz 

(The  propagation  parameters  are  c/v  =  1 . 1 ,  a-  2.7  dB/Mm,  h  =  8 1  km;  the  azimuthal  angle  is  (p=  0°;  and  the  dipole  moment  is  p  =  1  Am.) 


-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Magnetic  Field  Produced  by  a  HED  Under  Typical  Nighttime  Propagation  Conditions  at  300  Hz 

(The  propagation  parameters  are  civ  =  1.1,  a-  2.7  dB/Mm,  h  -  81  km;  the  azimuthal  angle  is  (p~  90°;  and  the  dipole  moment  is  p  =  1  Am.) 


-  Spherical  waveguide 

-  Bannister  (dir.) 

-  Bannister  (dir.+ind.) 


Figure  5-16c.  Comparison  of  Spherical  Waveguide  and  Bannister’s  ELF  Propagation  Formulas  for  the  Magnitude  of  the 
Azimuthal  Magnetic  Field  Produced  by  a  HED  Under  Typical  Nighttime  Propagation  Conditions  at  300  Hz 

(The  propagation  parameters  are  c/v  =  1.1,  a=  2.7  dB/Mm,  h  =  81  km;  the  azimuthal  angle  is  (p=  0°;  and  the  dipole  moment  is  p  =  1  Am.) 


6.  SUMMARY  AND  CONCLUSIONS 


ELF  propagation  formulas  have  been  derived  for  dipole  sources  radiating  in  a  spherical 
earth-ionosphere  waveguide,  where  the  waveguide  boundaries  are  approximated  as  scalar  surface 
impedances.  The  range  dependencies  in  these  formulas  involve  the  Legendre  function  of  the 
first  kind  of  complex  degree  and  order  zero  or  one  of  its  first  two  derivatives.  Several 
approximations  to  the  Legendre  function  were  derived.  Through  use  of  the  earth-flattening 
approximation  to  the  Legendre  function,  it  was  shown  how  the  spherical  waveguide  formulas 
reduce  to  Bannister’s  simplified  propagation  formulas.  In  addition,  through  use  of  another 
approximation  to  the  Legendre  function  that  is  suitable  for  antipodal  ranges,  it  was  also  shown 
how  the  spherical  waveguide  formulas  reduce  to  formulas  that  are  similar  to  the  antipode- 
centered  formulas  that  were  previously  derived  by  the  author.  Numerical  results  focused  on  the 
quasi-TEM  mode,  the  only  propagating  mode  in  the  30-Hz  to  300-Hz  band. 

Comparisons  of  Bannister’s  formulas  with  the  spherical  waveguide  formulas  were  made  for 
VED  and  HED  sources  for  the  case  where  both  the  source  and  field  points  are  located  along  the 
surface  of  the  earth.  The  300-Hz  results  have  shown  that  Bannister’s  total  field  formulas 
produce  a  good  agreement  with  the  spherical  waveguide  formulas  to  distances  within  less  than 
1  Mm  from  the  antipode.  This  agreement  degrades  with  decreasing  frequency  because  of  lower 
attenuation.  It  was  also  found  that  the  agreement  between  the  formulas  is  worse  for  the  field 
components  proportional  to  the  derivatives  of  the  Legendre  function.  Bannister’s  direct  great- 
circle  path  formulas  do  not  agree  as  closely  with  the  spherical  waveguide  formulas  because  they 
do  not  account  for  the  interference  produced  by  the  indirect  great-circle  path  contribution  in  the 
vicinity  of  the  antipode. 

This  report  also  presented  derivations  of  the  quasi-TEM  spherical  waveguide  formulas  for  a 
HED  located  above  an  anisotropic  surface  impedance.  These  modified  spherical  waveguide 
formulas  are  extensions  of  Wolkoff  and  Kraimer’s  formulas  (references  4  and  5)  and  include  the 
antenna  pattern  factors  that  were  previously  derived  in  reference  6.  These  formulas  will  be  useful 
for  the  prediction  of  the  EM  fields  radiated  by  the  U.  S.  Navy’s  four  transmitting  antennas  at 
antipodal  ranges. 

At  antipodal  ranges,  the  spherical  model  of  the  earth-ionosphere  waveguide  may  not  be 
suitable  because  of  the  variation  in  the  ionospheric  reflection  height  along  the  great-circle  paths. 
As  a  result,  the  spherical  waveguide  formulas  need  to  be  modified  in  order  to  account  for  the 
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variations  in  the  propagation  parameters  along  the  direct  and  indirect  great-circle  paths  as 
functions  of  solar  elevation.  Wolkoff  and  Casey  (reference  15)  present  formulas  for  the  effective 
propagation  parameters  along  a  nonuniform  earth -ionosphere  waveguide.  A  development  of 
ELF  propagation  formulas  for  a  nonspherical  earth -ionosphere  waveguide  will  be  the  subject  of  a 
future  investigation. 
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APPENDIX  A 

ELECTROMAGNETIC  FIELDS  IN  TERMS  OF  POTENTIALS 
IN  SPHERICAL  COORDINATES 

A.1  ELECTRIC  AND  MAGNETIC  VECTOR  POTENTIALS 

Consider  a  homogeneous,  isotropic  medium  of  permittivity  e,  permeability  p,  and 
conductivity  <7.  In  this  medium,  there  exist  electric  current  and  charge  distributions  J  and  p, 
respectively,  and  fictitious  magnetic  current  and  charge  distributions  M  and  pm,  respectively.  For 

a  time-harmonic  dependence  e-7®*,  Maxwell’s  equations  are  written  as  follows  (references  19, 20, 
and  21): 


V  x  E  =  -  jwpH  -  M , 

(A- la) 

V  x  H  =  jmeE  +  J , 

(A- lb) 

V  •  E  =  — 

V  L,  £  , 

(A-lc) 

0.1=5- 

II 

tu 

> 

(A-ld) 

In  the  above  expressions,  E  and  H  denote  the  electric  and  magnetic  fields,  respectively,  (0  =  2nf 
is  the  angular  frequency,  and  e  =  e  -  jo/co  is  the  effective  permeability  of  the  medium.  The 
electric  and  magnetic  fields  can  each  be  expressed  in  terms  of  a  pair  of  vector  potentials.  The 
pair  of  vector  potentials  that  are  most  commonly  used  are  the  magnetic  vector  potential  A  and 
the  electric  vector  potential  F.  The  magnetic  vector  potential  is  associated  with  electric  sources 
(M= 0,  pm  =  0)  and  the  electric  vector  potential  is  associated  with  magnetic  sources  (J=0, 

P  =  0). 


Consider  the  case  when  only  electric  sources  exist  ( M  =  0,  pm  =  0).  From  equation  (A-ld), 
the  magnetic  field  is  solenoidal  and  can  be  represented  as  the  curl  of  another  vector  .4, 

H  =  jiWxA.  (A-2) 

The  multiplicative  factor  1/p  has  been  included  to  conform  with  the  conventional  definition.  The 
substitution  of  equation  (A-2)  into  (A- la)  yields 


A-l 


Vx|fi  +  jwA )  =  0. 


(A-3) 


In  expression  (A-3),  the  term  within  the  parentheses  is  irrotational  and  is,  therefore,  equal  to  the 
gradient  of  a  scalar  potential  &e.  Thus, 


E  =  - ju)A  -  V<Pe  . 


(A-4) 


A  negative  sign  in  front  of  the  potential  has  been  included  to  conform  with  the  convention.  An 
alternate  expression  for  the  electric  field  can  be  obtained  from  expressions  (A- la)  and  (A-2)  as 


E  =  — i—  (v 
jaofie  v 


x  V  x  A  - 


(A-5) 


The  substitution  of  equations  (A-2)  and  (A-4)  into  (A- lb)  yields  the  general  equation  for  the 
magnetic  vector  potential,  i.e., 

V  x  V  x  A  -  k2A  =  -  j(o^eV<Pe  +  fij  ,  (A-6) 

where  k  =  co^fjHF  is  the  wave  number  of  the  medium.  From  a  theorem  in  vector  calculus 
(reference  22),  a  vector  is  uniquely  specified  by  giving  its  divergence  and  curl  within  a  given 
region  and  its  normal  component  over  the  boundary.  Therefore,  in  the  definition  of  the  magnetic 
vector  potential  in  (A-2),  A  is  underdetermined  because  its  divergence  has  not  been  specified. 
The  specification  of  the  divergence  of  A  is  referred  to  as  the  gauge  condition  (reference  23).  The 
gauge  condition  chosen  usually  removes  the  scalar  potential  term  from  equation  (A-6),  thereby 
simplifying  the  equation.  For  the  problems  of  interest  in  this  report,  spherical  coordinates  (r,  6, 
(p)  are  used  and  only  radially -directed  current  sources  are  of  interest. 

Consider  the  case  when  only  magnetic  sources  exist  (J  =  0,  p  =  0).  From  equation  (A-lc), 
the  electric  field  is  solenoidal  and  can  be  represented  as  the  curl  of  another  vector  F,  i.e., 

E  =  -  \  V  x  F  .  (A-7) 

t 

The  multiplicative  factor  -Me  has  been  included  to  conform  with  the  conventional  expression. 
The  substitution  of  expression  (A-7)  into  (A- lb)  yields 

Vx(ff  +  j(')F)  =  0  .  (A-8) 


A-2 


In  equation  (A-8),  the  term  within  the  parentheses  is  irrotational  and  is  therefore  equal  to  the 
gradient  of  a  scalar  potential  <Ph.  Thus, 

H  =  - jwF  -  V  <Ph  .  (A-9) 

A  negative  sign  in  front  of  the  scalar  potential  has  been  included  to  conform  to  the  convention.  An 
alternate  expression  for  the  magnetic  field  can  be  obtained  from  expressions  (A- la)  and  (A-7)  as 

H  =  —  (v  x  V  x  F  -  eM]  .  (A-10) 

jwfis  V  1 

The  substitution  of  equations  (A-7)  and  (A-9)  into  (A- la)  yields  the  general  equation  for  the 
electric  vector  potential,  thus, 

V  x  V  x  F  -  k2F  =  -  ja)[xsV<Ph  +  eM .  (A-ll) 

In  a  homogeneous  source-free  region  (i.e.,  J=M=  0),  the  electromagnetic  (EM)  fields  can 
be  expressed  in  terms  of  both  vector  potentials.  From  expressions  (A-5)  and  (A-7),  the  electric 
field  is  given  by 

E  =  — - — V  xVxA-iVxF.  (A- 12) 

JCOfiE  £ 

From  formulas  (A-2)  and  (A-10),  the  magnetic  field  is  given  by 

H  =  }j  VxA  +  -L  VxVxF.  (A- 13) 

"  ycojue 

Equations  (A- 12)  and  (A- 13)  are  the  general  formulas  for  the  EM  fields  in  terms  of  the  vector 
potentials  in  a  source-free  region. 

In  the  representation  of  the  EM  fields  in  terms  of  potentials,  suitable  expressions  can  be 
obtained  from  only  one  scalar  component  of  both  A  andF.  In  spherical  coordinates,  suitable 
choices  for  A  and  F  are  given  in  terms  of  their  radial  components  as 

A  =  rAr,  (A- 14a) 

and 

F  =  rFr  ,  (A-14b) 


A-3 


where  r  denotes  the  unit  vector  along  the  radial  direction  (see  figure  4-1),  and  the  radial 
components  Ar  and  Fr  are  scalars.  The  substitution  of  equation  (A- 14a)  into  (A-6)  with  J=  0 
yields 

VxVx(r^)- k\r  Ar)  =  -  j(opieV<Pe  .  (A-15) 

Similarly,  the  substitution  of  (A- 14b)  into  (A-l  1)  with  M  =  0  results  in 

V  x  V  x  (jr  Fr)  -  k2(f  Fr )  =  -  jiofxsW d?h  .  (A-16) 


To  eliminate  the  scalar  functions  <Pe  and  <Ph  in  equations  (A-15)  and  (A-16),  respectively, 
gauge  conditions  for  A  and  F  must  be  defined.  To  obtain  the  appropriate  gauge  condition  for 
each  vector  potential,  one  must  first  look  at  the  three  spherical  components  of  each  vector 
potential  equation.  For  example,  in  equation  (A-15),  the  r,  6,  and  cp  components  are 


1 


r2  sin  6  dd 


ndAr 

sm9le, 


+ 


d2A, 


r2  sin2  d  dtp2 


+ 


k2A. 


■■  JO)U£ 


d<PC 


dr 


(A- 17a) 


d2A,  .  „  d®e 


and 


d2Ar 

drdq) 


=  - j  at  pie 


d 0e 
dtp 


(A-l  7b) 

(A- 17c) 


The  6  and  (p  components  of  equation  (A-15),  given  by  expressions  (A- 17b)  and  (A- 17c), 
respectively,  are  satisfied  identically  under  the  following  gauge  condition: 


dA 

=  -jcopi£<Pe 


(A- 18) 


The  substitution  of  the  above  gauge  condition  into  the  radial  component  equation  (A- 17a)  yields 


d2A 


1 


dr2  r 2  sin  6  dd 


„  0Ar 

sm6le, 


1  $  Ar 

H r — ; — — r — y  +  k  A  —  0  . 

r 2  sin2  6  d(p 2 


(A- 19) 


The  first  term  in  the  above  expression  can  be  rewritten  as 


d  \  _  1  d  2  d 
dr2  ~  r  dr  r  dr 


(A-20) 


A-4 


The  substitution  of  (A-20)  into  (A- 19)  followed  by  multiplication  by  Hr  yields 


J _d_ 
r2  dr 


r2_d_K) 

7  dr\r  ) 


+ 


1 


r2  sin  0  dd 


sin6m t 


+ 


i 


r2  sin2  6 


d(p 2' 


(  A,.)  J  A,) 

[-/J  +  k2y-jr- J  =  0 .  (A-21) 


The  above  equation  may  be  written  more  concisely  as 
A, 


(v2  +  k2)  y  =  0  . 


(A-22) 


Therefore,  A  Jr  is  a  solution  of  the  scalar  Helmholtz  equation.  Following  a  similar  procedure  for 
equation  (A- 16)  in  spherical  coordinates  results  in 


(v2  +  k2)  -f  =  0 

provided  that  the  following  gauge  condition  is  satisfied: 


(A-23) 


dF, 


dr 


-  =  -ja)Ll£0h  . 


(A-24) 


Formulas  for  the  EM  field  components  in  terms  of  Ar  and  Fr  can  be  obtained  through 
substitution  of  expressions  (A-14a)  and  (A-14b)  into  (A-12)  and  (A-13).  In  spherical 
coordinates,  the  field  components  are  given  as  follows: 


£  = _ 1 _ 

r  j(OfX£ 


dr2  lr’ 


Ea  = 


_  1  1 


d2A, 
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dF. 
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Hr  =  ^:\T2+k2 
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r  +  ^  1 


b2F. 


ZJ  _  J; _ l_  _ 

6  ju  r  sin  0  dtp  r  jWfJl g  r  drdfl  ’ 


(A-25a) 

(A-25b) 

(A-25c) 

(A-25d) 

(A-25e) 


and 
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(A-25f) 


IT  _ _ x  r  |  i _ 1  r 

nr  dd  jome  r  sin  6  drdcp 

When  Fr  =  0,  then  Hr  =  0,  and  the  fields  are  referred  to  as  transverse  magnetic  (TM)  with  respect 
to  the  radial  direction.  Similarly,  when  Ar  =  0,  then  Er=  0  and  the  fields  are  referred  to  as 
transverse  electric  (TE)  with  respect  to  the  radial  direction. 


A.2  HERTZ  VECTOR  POTENTIALS 

An  alternate  set  of  vector  potentials  that  are  frequently  used  in  EMs  are  the  electric-  and 
magnetic-Hertz  vectors  Ke  and  ith,  respectively.  These  vectors  are  defined  in  terms  of  the 
magnetic-  and  electric -vector  potentials  as  follows  (reference  20) 

A  =  jiou£ne ,  (A-26a) 

and 

F  s  ju)fA't3ih  .  (A-26b) 

Therefore,  in  the  frequency  domain,  the  replacement  of  A  and  F  by  the  Hertz  vectors  amounts  to 
nothing  more  than  the  inclusion  of  a  complex  constant  factor.  The  substitution  of  the  above 
definitions  into  equations  (A- 12)  and  (A- 13)  yields 

£  =  VxVx«c- jmfi  V  x  jih ,  (A-27) 

and 

H  =  j(D£  Vxjre  +  VxVxjrft  .  (A-28) 

The  above  formulas  are  the  EM  fields  in  a  homogeneous  source-free  region  in  terms  of  the  Hertz 
vectors. 

To  represent  the  EM  fields  in  spherical  coordinates  in  terms  of  the  Hertz  vectors,  suitable 
expressions  can  be  obtained  from  only  one  scalar  component  of  both^4  and  F.  Following 
expressions  (A-14a)  and  (A- 14b),  the  following  are  defined: 

3ie  =  fjrer  =  frWc,  (A-29a) 

and 
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(A-29b) 


jth  =  rn;h  =  rrWh  . 

With  the  above  definitions,  the  same  procedure  after  expressions  (A-14a)  and  (A-14b)  is 
followed,  resulting  in 

S/2We  +  k2We  =  0,  (A-30a) 

and 

V2Wh  +  k2Vh  =  o.  (A-30b) 


The  terms  We  and  yh  are  sometimes  referred  to  as  the  Debye  potentials  (reference  23). 


Formulas  for  the  EM  field  components  in  terms  of  *Fe  and  *¥h  can  be  obtained  through  the 
substitution  of  expressions  (A-29a)  and  (A-29b)  into  (A-27)  and  (A-28).  In  spherical 
coordinates,  the  field  components  are  given  as  follows: 


Er  = 


(l^jw 


_  i  d2(rWe)  jeon  dWh 
Eq  ~  ~r  ~drdW~  ~  sm“0  * 

l  d2{rWe )  .  dWh 

v  ~  r  sin  9  drdcp  +  86 


Hr  = 
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dr2 


+  k2  ( rWh ), 


_  j(OE  8We  !  8{rWh ) 
Hd ~  smD  +  7  drdd  ’ 


and 


v  J0}£  dd  r  sin  0  drdq) 


(A-31a) 

(A-31b) 

(A-31c) 

(A-31d) 

(A-31e) 

(A-31f) 


In  the  above  formulas,  ner-r  We  and  Jt^  =  r  Wh .  When  Jif  =  0,  the  fields  are  TM  with  respect 
to  the  radial  direction.  Similarly,  when  Jthr  =  0,  the  fields  are  TE  with  respect  to  the  radial 
direction. 
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APPENDIX  B 

SOLUTION  OF  THE  HELMHOLTZ  EQUATION 
IN  SPHERICAL  COORDINATES 


In  spherical  coordinates,  it  is  shown  in  appendix  A  that  the  electromagnetic  (EM)  fields  can 
be  expressed  in  terms  of  the  two  Debye  potentials,  and  xVh.  In  addition,  at  a  source-free 
point,  it  was  shown  that  these  two  scalar  functions  each  satisfy  the  Helmholtz  equation,  namely 

V2W  +  k2W  =  0.  (B-l) 


In  this  report,  because  the  transverse  magnetic  (TM)  and  transverse  electric  (TE)  fields  are 
derived  for  axisymmetric  sources,  the  Debye  potentials  are  determined  under  axisymmetric 
conditions.  Therefore,  under  axisymmetric  conditions,  the  Helmholtz  equation  is  expressed  in 
spherical  coordinates  as 


1 


r2sin  0  dd 


sin  0 


+  k2W  =  0  , 


0-2) 


where  *F=  «F(r,0). 


To  solve  equation  (B-2),  the  separation  of  variables  method  is  applied.  In  this  method,  let 


'F(r,6)=R(r)T(6). 


0-3) 


The  substitution  of  expression  (B-3)  into  (B-2)  followed  with  the  multiplication  of  each  side  of 
the  equation  by  r  2/RT  yields 


df[r2f)+T£rele[sine%) 


-L  d-[  r2 

R 


+  k2r2  =  0  . 


0-4) 


The  above  equation  can  be  rewritten  as  an  ^--dependent  part  on  one  side  and  a  0-dependent  part 
on  the  other  side,  each  of  which  can  be  set  equal  to  an  arbitrary  constant  v(v+l)  resulting  in 


+  k2r2  =  TSe  le[s'n  e$)m  v(v  +1)  • 


(B-5) 


Expression  0-5)  can  be  decomposed  into  r-dependent  and  0-dependent  differential  equations  as 
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dr{r2^lr)+  k2r2~v(v+V  R  =  0’ 


(B-6a) 


and 


1  d 
sin  6  dd 


sin  6 


dT 

dd 


+  v(v  +1)7  =  0  . 


(B-6b) 


The  solution  of  equation  (B-6a)  is  given  by 

R(r)  =  A h[}\kr)  +  Bh^(kr) ,  (B-7) 

where  h[]\kr)  and  h^\kr)  are  the  spherical  Hankel  functions  of  the  first  and  second  kinds, 
respectively,  of  order  v ,  and^f  and  B  are  arbitrary  constants.  The  spherical  Hankel  functions  are 
expressed  in  terms  of  the  cylindrical  Hankel  functions  as  (reference  24) 

hf{kr)  =  ^ ffW  m(kr) ,  n=  1,2.  (B-8) 

The  solution  of  (B-6b)  is  given  as  follows: 

T(6)  =  CPV(-  cos  6)  +  D Qv(-  cos  6) ,  (B-9) 

where  Pv  and  Qv  are  Legendre  functions  of  the  first  and  second  kinds,  respectively,  of  degree  v 
and  order  zero.  It  should  be  noted  that  Pv  is  singular  when  its  argument  is  -1  (0=  0),  while  Qv  is 

singular  when  its  argument  is  ±1  (6=  0,  it).  The  minus  sign  has  been  included  in  the  argument 
so  that  Pv  is  singular  at  0=  0. 


For  the  three  radiating  sources  under  consideration  in  this  report,  the  EM  fields  are  finite 
everywhere  except  at  the  source  location  (r  =  rs,  6=  0).  Therefore,  Qv  is  not  an  acceptable 

solution  of  equation  (B-6b)  and  D  =  0.  Thus,  the  substitution  of  formulas  (B-7)  and  (B-9)  into 
(B-3)  yields 


Vir, 6)  = 


Ah{}\kr)  +  Bh^\kr) 


Py(-  COS  6) 


(B-10) 


Two  of  the  three  constants  A,  B,  and  v  are  determined  from  the  boundary  conditions  at  the  inner 
(r  =  a)  and  outer  (r  =  a  +  h)  surfaces. 


B-2 


APPENDIX  C 

TRANSVERSE  MAGNETIC  (TM)  AND  TRANSVERSE  ELECTRIC  (TE) 
MODAL  EXPANSIONS  IN  A  SPHERICAL  WAVEGUIDE 


Consider  the  spherical  earth-ionosphere  waveguide  described  in  figure  1-1.  The  main  interest 
of  this  report  is  the  determination  of  the  electromagnetic  (EM)  fields  radiated  in  the  waveguide 
(a  =  r  =  a  +  h)  by  a  dipole  source  located  along  the  north  pole  {6=  0)  of  the  spherical  coordinate 
system.  The  waveguide  is  filled  with  air  and  has  a  permittivity  e0  =  8.854  x  10' 12  F/m,  a 
permeability  [i0  =  471  x  10'7  H/m,  and  an  intrinsic  wave  number  of  kQ  -  (atf\T^~0  ,  where  (0  is  the 
angular  frequency  in  radians  per  second  (rad/s).  The  spherical  waveguide  boundaries  along  the 
earth  (r  =  a)  and  the  ionosphere  (r  =  a  +  h)  are  characterized  by  scalar  surface  impedances. 

These  impedances  are  approximations  of  the  general  tensor  surface  impedances  that  are  more 
representative  of  these  boundaries.  However,  the  scalar  impedances  allow  for  tractable  solutions 
of  the  field  equations. 

In  this  appendix,  the  transverse  magnetic  (TM)  and  transverse  electric  (TE)  modes  that 
radiate  in  a  spherical  earth-ionosphere  waveguide  are  derived  under  the  assumption  that  the 
fields  are  axisymmetric.  The  axisymmetric  TM  modes  are  launched  by  a  vertical  electric  dipole 
(VED)  located  along  the  north  pole  of  the  waveguide.  Similarly,  the  axisymmetric  TE  modes  are 
launched  by  a  vertical  magnetic  dipole  (VMD)  located  in  the  waveguide  along  the  north  pole. 

As  shown  in  appendix  H,  the  HED  fields  involve  the  superposition  of  both  the  TM  and  TE 
modes. 


C.1  TM  MODES 


In  appendix  A,  it  is  shown  that  the  TM  fields  with  respect  to  the  radial  direction  can  be 
obtained  from  the  radial  component  of  the  electric -Hertz  vector  ne.  Under  axisymmetric 
conditions  (?/?<p  =  0),  the  EM  fields  are  obtained  from  equations  (A-31)  and  are  expressed  as 


Er  = 


(  «■»  2 


—  +  k2 
dr*  ° 


( rWe )  =  - 


f)  |  dWe 
r  sin  6  30\S^n  ^  dd  ]  ’ 
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_  j  d2(rWe) 
e  r  drdd 


(C-la) 


(C-lb) 


C-l 


and 


3We 

Hv  =  -JCO£o  -je~  •  (C-Ic) 

Note  that  the  three  remaining  field  components  vanish,  namely  E(p  =  Hr  =  Hq~  0.  In  the  above 
formulas,  the  scalar  function  *Fe  is  a  Debye  potential  and  is  related  to  ner  as 

We=^-,  (C-2) 


where  *FC  is  a  solution  of  the  scalar  Helmholtz  equation  given  by  (A-30a).  The  solution  of  the 
Helmholtz  equation  is  derived  in  appendix  B  and  is  given  by  expression  (B-10).  This  solution 
can  be  rewritten  as  follows: 


xpe  — 


h{y\kj)  +  Bevh^\k0r) 


PI-  cos  6) , 


(C-3) 


where  Bev  and  v  are  constants.  The  constants  Bev  and  v  are  determined  from  the  boundary 

conditions  at  the  inner  and  outer  boundaries  of  the  waveguide  that  are  given  below.  In  expression 
(C-3),  note  that  a  constant  multiplying  h^\k0r)  has  not  been  included  because  it  is  unnecessary 

in  satisfying  the  boundary  conditions. 


Under  the  assumptions  of  a  homogeneous,  isotropic  earth  (inner  sphere)  and  a  homogeneous, 
isotropic  ionosphere  (outer  region)  of  constant  height  h,  the  surface  EM  fields  satisfy  the 
following  boundary  conditions: 
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(C-4a) 


(C-4b) 


where  Zg  and  Ze  denote  the  surface  impedances  along  the  air-earth  and  air-ionosphere  interfaces, 

respectively.  The  substitution  of  the  field  component  expressions  (C-lb)  and  (C-lc)  into  the 
boundary  conditions  (C-4a)  and  (C-4b)  yields 
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(a  +  h) 
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=  -j(oe  Z  -Tgr 
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r  =  a  +  h 


(C-5b) 


Note  that  the  6  derivatives  in  the  above  boundary  conditions  are  continuous  across  each 
interface.  In  addition,  the  surface  impedances  Zg  and  Ze  are  each  constant  along  the  surfaces  r  = 

a  and  r  =  a  +  h,  respectively.  Therefore,  the  boundary  conditions  (C-5a)  and  (C-5b)  reduce  to 


l  d(rWe ) 
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(C-6a) 


(C-6b) 


where  Ag  =  Z^r\Q  and  Ae  =  ZJr\0  denote  the  normalized  surface  impedances  of  the  inner  and 
outer  boundaries,  respectively,  and  r\0  =  p.0/e0  is  the  intrinsic  impedance  of  free  space. 


The  substitution  of  expression  (C-3)  for  We  into  the  boundary  condition  (C-6a)  can  be  solved 
for  Bev  to  yield 


(C-7a) 


where  u  =  k0r  and  ug  =  k0a.  Following  a  similar  procedure,  the  substitution  of  expression  (C-3) 
for  *Fe  into  the  boundary  condition  (C-6b)  can  be  solved  for  Bev  to  yield 
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(C-7b) 
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where  u{  =  kQ(a  +  h).  The  TM  field  components  must  satisfy  the  boundary  conditions  (C-4a)  and 
(C-4b)  simultaneously.  Therefore,  the  constant  Bcv  in  expression  (C-7a)  is  set  equal  to  (C-7b) 
and  yields 
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The  above  expression  denotes  the  characteristic  or  modal  equation  for  the  TM  modes  in  a 
spherical  earth- ionosphere  waveguide.  The  solution  of  expression  (C-8)  yields  a  discrete  set  of 
eigenvalues  v„,  where  n  denotes  an  integer  index  and  corresponds  to  a  particular  mode. 


Given  the  infinite  number  of  discrete  eigenvalues  that  satisfy  the  characteristic  equation  (C- 
8),  *Pe  can  be  expressed  in  terms  of  a  modal  expansion  as 


W( 


=  2  Wev  ( r,9 )  =  2  K  Rev  (u)  Pv  (-  cos  d ) 


«  =  0 


w  =  0 


(C-9) 


where  Re  is  defined  as 

v  n 

K  («)  =  hv\11)  +  K  hv\u)  • 

n  n  n  n 


(C-10) 


In  addition,  Ae  is  referred  to  here  as  the  modal  excitation  coefficient  and  depends  on  the  source 

strength.  The  modal  indices  n  in  expression  (C-9)  correspond  to  the  TM  mode  indices  for  the 
parallel-plate  waveguide  model  as  defined  in  section  2.1.  Expressions  for  the  modal  expansions 
of  the  TM  fields  are  determined  through  the  substitution  of  the  formula  (C-9)  into  equations 
(C-l)  and  are  given  as  follows: 
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The  expression  for  Er  was  obtained  through  use  of  the  relation 


j^  +  ki)  [rft»(V>]  = h «(*/) .  m  =  1, 2  . 


(C-12) 


The  above  formula  can  be  obtained  from  expression  (B-6a)  through  the  following  equivalence: 
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r2fr  hf\kQr)  =  r  jp{rh^\k  0r)\ ,  m=  1,2. 


(C-13) 


C.2  TE  MODES 


In  appendix  A,  it  is  shown  that  the  TE  fields  with  respect  to  the  radial  direction  can  be 
obtained  from  the  radial  component  of  the  magnetic  Hertz  vector  ich.  Under  axisymmetric 
conditions  (?/?  (p  =  0),  the  EM  fields  are  obtained  from  equations  (A-31)  and  are  expressed  as 
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Note  that  the  three  remaining  field  components  vanish,  namely  H,p  =  Er  =  Eq  =  0.  In  the  above 
formulas,  the  scalar  function  *Fe  is  a  Debye  potential  and  is  related  to  as 


xph  _ 


(C-15) 


where  *Fh  is  a  solution  of  the  scalar  Helmholtz  equation  given  by  (A-30b).  The  solution  of  the 
Helmholtz  equation  is  derived  in  appendix  B  and  is  given  by  expression  (B-10).  The  application 
of  the  solution  (B-10)  to  equation  (A-30b)  yields 
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where  Bhv  and  v  are  constants.  The  constants  Bhv  and  v  are  determined  from  the  boundary 
conditions  at  the  inner  and  outer  waveguide  boundaries  that  are  given  below.  As  was  previously 
mentioned  for  the  TM  modes,  a  constant  multiplying  h[}\kor)  has  not  been  included  because  it 
is  unnecessary  in  satisfying  the  boundary  conditions. 


Under  the  assumptions  of  a  homogeneous,  isotropic  earth  (inner  sphere)  and  a  homogeneous, 
isotropic  ionosphere  (outer  region)  of  constant  height  h,  the  TE  surface  EM  fields  satisfy  the 
boundary  conditions 
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where  Zg  and  Zh  denote  the  surface  impedances  along  the  air -earth  and  air- ionosphere  interfaces, 
respectively.  The  surface  impedances  Zg  and  Zh  given  above  are  generally  different  from  those 
corresponding  to  the  TM  modes  given  in  expressions  (C-4a)  and  (C-4b)  because  of  the 
anisotropy  of  the  earth  and  ionosphere.  In  this  report,  the  earth  and  ionosphere  are  each 
represented  as  a  constant  surface  impedance  that  is  independent  of  the  mode  index,  with  the 
surface  impedance  of  the  ionosphere  different  for  TE  and  TM  modes. 


The  substitution  of  the  field  component  expressions  (C-14b)  and  (C-14c)  into  the  boundary 
conditions  (C-17a)  and  (C-17b)  yields 
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Note  that  the  6  derivatives  in  the  above  boundary  conditions  are  continuous  across  each  interface. 
In  addition,  the  surface  impedances  Zg  and  Zh  are  each  constant  along  the  surfaces  r  =  a  and 
r  =  a  +  h,  respectively.  Therefore,  the  boundary  conditions  (C-18a)  and  (C-18b)  reduce  to 
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where  Ag  and  Ah  =  Zh/rj0  denote  the  normalized  surface  impedances  of  the  inner  and  outer 
waveguide  boundaries,  respectively. 


The  substitution  of  expression  (C-16)  for  *Fh  into  the  boundary  condition  (C-19a)  can  be 
solved  for  B1',  to  yield 


(C-20a) 


where  u  and  ug  were  defined  in  the  TM  mode  derivation.  Following  a  similar  procedure,  the 
substitution  of  expression  (C-16)  for  *PA  into  the  boundary  condition  (C-19b)  can  be  solved  for 
B\  to  yield 


(C-20b) 


where  u{  was  defined  in  the  TM  mode  derivation.  Note  that  expression  (C-20a)  for  the  TE  field 
amplitude  coefficient  B*  is  equivalent  to  expression  (C-7a)  for  the  TM  field  amplitude 
coefficient  Bev  if  Ag  is  replaced  by  1  /Ag.  Similarly,  expression  (C-20b)  for  B\,  is  equivalent  to 
expression  (C-7b)  for  Bev  if  Ae  is  replaced  by  1  IAh.  Because  the  TE  field  components  must 
satisfy  the  boundary  conditions  (C-17a)  and  (C-17b)  simultaneously,  the  constant  B*  in 
expression  (C-20a)  is  set  equal  to  expression  (C-20b)  and  yields 
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=  1 .  (C-21) 


The  above  expression  denotes  the  characteristic  or  modal  equation  for  the  TE  modes  in  a 
spherical  earth-ionosphere  waveguide.  Note  that  the  above  TE  characteristic  equation  is 
equivalent  to  the  TM  characteristic  equation  (C-8)  if  Ag  is  replaced  by  1  /Ag  and  Ae  is  replaced  by 
\/Ah.  The  solution  of  equation  (C-21)  yields  a  discrete  set  of  eigenvalues  vm ,  where  m  denotes 
an  integer  index  and  corresponds  to  a  particular  mode. 


Given  the  infinite  number  of  discrete  eigenvalues  that  satisfy  the  characteristic  equation 
(C-21),  *Fh  can  be  expressed  in  terms  of  a  modal  expansion  as 


uu  uu 

Vh=  2  W*  o r,e )  Rhv  (//)  pv  (-cos  d). 


m=  1 


(C-22) 


where  R~,  (u)  is  defined  as 

v  m 


Rh-(u)  =  h(}\u)  +  B*  hf(u) 


(C-23) 


In  formula  (C-22),  Ah-  is  the  TE  mode  excitation  coefficient  and  depends  on  the  source  strength. 
The  modal  indices  m  in  expression  (C-22)  correspond  to  the  TE  mode  indices  for  the  parallel-plate 
waveguide  model  as  defined  in  section  2.2.  Expressions  for  the  modal  expansions  of  the  TE  fields 
are  determined  through  the  substitution  of  formula  (C-22)  into  equations  (C-14)  and  are  given  as 


Ic  00 

Hr  =  -jf  2  K  vm(vm  +  !)  K  00  PV  (-  cos  6) , 

\  m  tu  tti 

He  =  k  |?K.(»)]  4s  P,.(-  9) , 

and 

00  d 

Ecr  =  j(0fAo  2  K  K  00  Jn  py  (-COS  6 )  . 

t  jji  =  j  m  m  V  U  m 


(C-24a) 

(C-24b) 


(C-24c) 


Note  that  the  expression  for  Hr  was  obtained  through  use  of  formula  (C-12). 
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APPENDIX  D 

PROOF  OF  ORTHOGONALITY  OF  RADIAL  FUNCTIONS 


In  appendix  A,  it  is  shown  that  the  transverse  magnetic  (TM)  and  the  transverse  electric  (TE) 
fields  in  a  spherical  waveguide  can  be  expressed  in  terms  of  the  Debye  potentials  We  and 
x¥h,  respectively.  In  appendix  C,  it  is  shown  that  the  radial  dependence  of  each  Debye  potential  is 
given  in  terms  of  a  function  Rv(u),  u  =  kj,  which  satisfies  the  following  differential  equation: 


The  above  equation  may  be  rewritten  as 
m|^(«/?v)  +  [«2_v(v+1)]rv  =  0. 


(P-2) 


To  prove  the  orthogonality  of  the  TM  mode  radial  functions,  consider  the  radial  functions 
Refpi)  and  Rev(u)  corresponding  to  two  different  TM  modes  that  satisfy  the  following  differential 

equations: 

u  iR £)  +  [«2  _  MjU  +i)]  =  o ,  (D-3a) 

and 

u  +  [«2  -  v(v  +1)]  Rev  =  0 .  (D-3b) 


From  formulas  (C-6),  the  boundary  conditions  satisfied  by  each  TM  mode  radial  function  are 
given  as 


I  d(uReJ 

II  g  du 


=  JA.ReM  > 


1 


d(uReJ 

n'  =~j\R‘n(ui), 


ui  du 


U  -  U- 


(D-4a) 


(D-4b) 


where  ug  =  kQa,  ut  =  k0(a + h ),  T)  is  a  complex  eigenvalue  that  satisfies  the  TM  mode 
characteristic  equation  (C-8),  and  Ag  and  Ae  denote  the  normalized  surface  impedances  of  the 
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inner  and  outer  boundaries,  respectively.  If  equation  (D-3a)  is  multiplied  by  Rev(u)  and  equation 
(D-3b)  is  multiplied  by  /?*(/<) ,  the  difference  of  the  resulting  equations  yields 


(««j)  -("R5)  J&K.) = hv+»  - +i> 


ne  V s 
Kv 


(P-5) 


If  the  above  formula  is  integrated  over  the  height  of  the  waveguide  followed  by  integration  by 
parts,  the  following  is  obtained: 


K)  sK)  “ 


d 

du\ 


\uRe\ 

U  -  U{ 

l  W 

It  =  lla 

f 

JU„ 


v( V  + 1 )  -  jw(  n  + 1 )  I  Re{  u)  Rev( u )  dll .  (D-6) 


The  application  of  the  boundary  conditions  (D-4)  to  the  above  result  yields  the  following 
orthogonality  relation: 


Rc;t(")  R%u)du  =  0 


[A  #  V  . 


(P-7) 


Because  the  above  result  is  independent  of  6,  the  radial  functions  may  be  replaced  by  the 
corresponding  Debye  potentials,  i.e., 


W^(u,6)  Wev(u,8)  da  -  0  ,  fi^v  . 


(D-8) 


To  prove  the  orthogonality  of  the  TE  mode  radial  functions,  consider  the  radial  functions 
and  corresponding  to  two  different  TE  modes  that  satisfy  the  differential  equations 

(D-3a)  and  (D-3b),  respectively,  where  the  superscript  e  is  replaced  by  h.  From  formulas  (C-19), 
the  boundary  conditions  satisfied  by  each  TE  mode  radial  function  are  given  as 


j  d(uR») 
“g  ~dii~ 


U  -  U 


and 


i  d(uR^) 

ui  du 


(D-9a) 


(D-9b) 
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where  Ah  denotes  the  normalized  surface  impedance  of  the  outer  boundary  and  77  is  a  complex 
eigenvalue  that  satisfies  the  TE  mode  characteristic  equation  (C-21).  If  the  procedure  below 
expressions  (D-4)  is  followed  for  the  two  TE  radial  functions  RHu)  and  Rhu) ,  the  result 
obtained  is 


K)iK) 


K)iK) 


M  =  «, 
U  =  Ug 


=  [v(v+l)-/r(ju+l)]  P  R%u)  Rhv{u)  du  . 


(D-10) 


The  application  of  the  boundary  conditions  (D-9)  to  the  above  result  yields  the  orthogonality 
relation 


f 

Ju, 


R^(u)  Ry(u)  du  =  0  , 


(D-ll) 


Because  the  above  result  is  independent  of  0,  the  radial  functions  may  be  replaced  by  the 
corresponding  Debye  potentials,  i.e., 


w*(u,d)  W*(u,0)  dll  =  0 


H*v  . 


(D-12) 


The  orthogonality  relations  (D-7)  and  (D- 1 1)  are  valid  provided  that  the  surface  impedances 
of  the  inner  and  outer  boundaries  are  independent  of  the  mode  indices.  If  the  above  procedure  is 
applied  to  two  radial  functions  Re„{u)  and  RHu)  ,  corresponding  to  one  TM  mode  and  one  TE 
mode,  respectively,  it  is  readily  seen  that  these  functions  are  not  orthogonal.  Therefore,  the  TM 
and  TE  mode  Debye  potentials  are  not  orthogonal  with  respect  to  each  other.  (The  evaluation  of 
the  integrals  (D-7)  and  (D-ll)  for  fi  =  vis  presented  in  appendix  I.) 
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APPENDIX  E 

EXCITATION  COEFFICIENTS  FOR  TRANSVERSE  MAGNETIC  (TM) 

AND  TRANSVERSE  ELECTRIC  (TE)  MODES 

In  the  derivation  of  the  modal  expansions  of  the  electromagnetic  (EM)  fields  in  a  spherical 
waveguide  given  in  appendix  C,  each  term  in  the  expansion  includes  an  excitation  coefficient. 
Each  modal  excitation  coefficient  depends  on  the  source  strength  and  provides  a  measure  of  the 
coupling  into  that  mode.  In  this  appendix,  the  excitation  coefficients  corresponding  to  the 
transverse  magnetic  (TM)  and  transverse  electric  (TE)  modes  in  a  spherical  waveguide  are 
derived.  The  excitation  coefficients  for  the  TM  modes  are  applied  in  the  derivation  of  the  vertical 
electric  dipole  (VED)  fields  in  appendix  F.  In  addition,  the  excitation  coefficients  for  the  TE 
modes  are  applied  in  the  derivation  of  the  vertical  magnetic  dipole  (VMD)  fields  in  appendix  G. 

E.1  TM  MODES 

Consider  a  VED  of  moment  p=I  dl  that  is  located  at  the  radial  distance  rs  =  a +zs  along  the 
6=0  axis  of  the  spherical  earth-ionosphere  waveguide  shown  in  figure  1-1.  As  shown  in 
appendix  F,  the  differential  equation  for  the  electric  Debye  potential  *Fe  is  given  by 

(v*  +  kl)v‘  =  -j£j-6(r-r,),  (E-l) 

where  d  denotes  the  Dirac  delta  function.  The  solution  of  the  above  differential  equation  at 
source-free  points  (r  ?  rs)  is  derived  in  appendix  C,  section  C.l,  and  is  given  in  terms  of  a  modal 

expansion  as 

00 

W%r,6)  =  2  Aev  Rev  («)  P  (-  cos  6) ,  (E-2) 

«=o  n  "  n 

where  Rev  has  been  defined  as 

Rev  («)  =  h^\u)  +  Be  h^\u)  .  (E-3) 

n  n  n  n 

In  the  above  formulas,  n  denotes  the  mode  index,  u  =  k0r,  Ae  is  the  excitation  coefficient  of  the 
nth  TM  mode,  Bev  is  a  constant  that  is  given  by  either  expression  (C-7a)  or  (C-7b),  h^\u)  and 

n  vn 

h(y\u)  denote  the  spherical  Hankel  functions  of  the  first  and  second  kinds,  respectively,  of  order 

tl 
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V„,  and  Pv  denotes  the  Legendre  function  of  the  first  kind  of  degree  v„  and  order  zero.  The 

v  n 

complex  constant  v„  is  determined  from  the  solution  of  the  characteristic  equation  (C-8). 


From  appendix  D,  the  orthogonality  condition  for  the  TM  mode  radial  functions  is  given  by 


/" 

Jll, 


K  («)  K  (»)  du  =  R%  2  dmtl , 


m  n 


(E-4) 


where  R$n  denotes  the  L2  norm  of  Rev  as  defined  in  appendix  F,  ug  =  k0a  ,  ui  =  k0(a+h),  and 


8mn  is  the  Kronecker  delta  defined  as 


o  _  I  1 ,  m  =  n 
mn  “10,  nit  n 


(E-5) 


If  each  side  of  expression  (E-2)  is  multiplied  by  Re  and  then  integrated  with  respect  to  u,  the 

vm 

application  of  the  orthogonality  condition  (E-4)  to  the  resulting  equation  yields  an  expression  for 
the  mode  excitation  coefficient  Aev  ,  namely 


Ae  - 
Av~ 


P 

Jun 


We(r,e )  Rev  (//)  du  ,  n  =  0,  1,  2, ... . 


||^J|2PV;1(-COS0) 

It  is  difficult  to  determine  Aev  from  equation  (E-6)  because  *Pe(r,  6)  is  not  known. 


(E-6) 


To  determine  the  excitation  coefficients  Aev  in  formula  (E-6),  let  6  approach  zero,  resulting  in 


At  =  — I 
"  K 


lim  (“ 

»-o  Ju 


lim  We(r,d)  Rev  (//)  du 


lim  P  (-  cos  6) 

6  ~*0  '  n 


(E-7) 


Because  both  the  numerator  and  denominator  are  singular  at  0  =  0,  the  above  expression  is  an 
indeterminate  form.  Then,  from  Magnus  and  Oberhettinger  (reference  25), 


sin  v„Jt  0 
lim  P.,  (-  cos  6)  =  — =7^  In  82 . 
6^0  V  ’  31 


(E-8) 


Before  proceeding  with  the  evaluation  of  the  numerator,  one  needs  to  look  more  closely  at  the 
behavior  of  the  Debye  potential  in  the  vicinity  of  the  source. 
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The  Debye  potential  *Fe  can  be  expressed  as  the  sum  of  a  particular  solution  W*  and  a 
homogeneous  solution  Wf;  of  the  differential  equation  (E-l).  The  particular  solution  satisfies  the 
inhomogeneous  equation  (E-l)  and  the  homogeneous  solution  satisfies  the  Helmholtz  equation 
(v2  +  &2j  W%  =  0.  The  homogenous  solution  is  finite  for  all  points  within  the  spherical 
waveguide  and  W*  is  singular  at  the  source  point  r  =  rs.  Therefore,  in  the  vicinity  of  the  source, 
the  particular  solution  W*  is  the  primary  influence  of  Ye  and,  thus, 


lim  We=  lim  Wf=  lim  We(r,9) . 

r  —  r,  r  —  r,  s  6-0  s 


(E-9) 


For  a  VED  of  moment  p  lying  in  free  space,  the  solution  of  equation  (E-l)  is  given  by 
(reference  23) 


We(r,d)  =  W*(r,e)  =  —£ - T,  , 

j4nwE0rs  K 


(E-10) 


where 


R  =  |  r  -  rs  =  r2  -  2rrs  cos  6  +  r2 


1/2 


(E-11) 


From  the  above  formula,  it  is  easy  to  see  that  We  is  singular  at  r  =  rs. 


Because  the  integrand  in  the  numerator  of  formula  (E-7)  is  finite  everywhere  except  at  the 
source  point  (r  =  rs,  9=  0),  its  only  contribution  to  the  integral  is  from  the  immediate 

neighborhood  of  the  source  point.  Therefore,  expression  (E-7)  can  be  simplified  to 


Mv  (« 


frs(l  +  e) 

x  lim  We(r,6)dr 

)  6-0  Jr(l-e)  sy 


Re  2  lim  (-  cos  9) 
KVn  6  —  0  V 


(E-l  2) 


where  u  has  been  replaced  by  k0r  and  e  «  1.  Note  that  Rev  ( u )  has  been  replaced  by  Rev  (u  ) 

n  y  n 

because  it  is  finite  at  the  source  point.  The  substitution  of  W*(r,6)  in  formula  (E-10)  into  the 
integral  of  equation  (E-l  2)  yields 


fr  (\  +  s) 

lim  yf(r,6)  dr 

6-0Jrs(i-e)  S 


-  lim 

j4nw£0rs  e  —  o 


i 


rsa  +  e ) 
(1-s) 


e~JkoR 

—R~dr- 


(E-13) 
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From  Sommerfeld  (reference  26),  let 


r=r(  1  +i]) 


(E-14) 


where  -  e<ij<  e  ,  e«  l .  The  substitution  of  formula  (E-14)  into  expression  (E- 11)  gives 


R  =  r. 


1  -  2(  1  +  ?7)  cos  0  +  ( 1  +  17)  - 


1/2 


(E-15) 


As  0->  0,  cos  0  =  1-  02/ 2.  Therefore,  R  may  be  approximated  as 

1 1/2 


/?  s  r„ 


=  r . 


l-2(l  +  r7)(l-02/2)  +  (l  +  >7): 


vj2  +  (1  +  i])  62 


1/2 


\J  t] 2  +  0“ 


(E-16) 


where  the  last  approximation  was  made  because  77  «  1 .  Therefore,  with  the  above  approximation 
fori?,  the  numerator  of  the  integrand  in  expression  (E-13)  can  be  approximated  as 


e  -ik°R  «  1 . 


(E-17) 


The  substitution  of  formulas  (E-14),  (E-16),  and  (E-17)  into  (E-13)  yields 

fr  (1  +  e)  p  f  (jji 

I  nre/„  a\  a~  _  ^  I  u'l 

Jrr 


frs(l  +  e)  p 

lim  I  W*(r,d)dr= - 

e-»o7r(i_e)  jAmoe 


-  lim  f  ^ 
A  \/»?2+02 


(E-18) 


From  a  table  of  integrals  (reference  27), 


chf 


\f  t}2  +  6' 


=  In 


77+  \/  7? 2  +  02 


(E-19) 


The  application  of  formula  (E-19)  to  the  integral  in  (E-18)  yields 


ch] 


Le  sfn2+  02 


In 


e  +  \/e2  +  02  -  In  -  £  +  \^£2  +  02 


(E-20) 


For  0«  e,  the  following  approximation  can  be  made: 


\/£2  +  02  =  £  \/ 1  +  (  0  /  £)  2  K  £ 


1  +  ^(0/£): 


=  £  + 
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Therefore,  from  the  above,  the  following  is  obtained: 


c  +  v^c2  +  02  s2e  +  — , 
2s 


and 


-  £  +  £2  +  02  «  —  . 

2e 


(E-21a) 

(E-21b) 


The  substitution  of  the  approximations  (E-21a)  and  (E-21b)  into  formula  (E-20)  yields 


l 


drj 


-e  V»l2+  02 


=  In 


2£  +  f  -in  f 
2c  I  1  2c 


=  In 


2c  +  fl2/2e 
d2/2e 


=  In 


(2  c) : 


02 


=  In  (2c)2 -In  62 


(E-22) 


If  the  limit  of  the  above  result  is  taken  as  6  0  (with  e  finite),  then 


lim  f  *L 


e 2 


=  -In  02  . 


(E-23) 


The  above  limit  is  the  same  as  the  one  given  by  Sommerfeld  (reference  26),  but  differs  in  sign 
from  the  one  given  by  Wait  (reference  1).  The  substitution  of  the  above  result  into  formula 
(E-18)  gives 

frs(l  +  e)  „  P  ", 

lim  I  'Pf(r,0)  dr  - - - - In  62  .  (E-24) 

e^oJrs(i-£)  j4iu»£0rs 


The  modal  excitation  coefficients  Ae  can  now  be  obtained  through  substitution  of  the 

Vn 

limiting  results  (E-8)  and  (E-24)  into  expression  (E-12),  resulting  in 


Ae  =  MoP  1  R''SUJ 

vn  4  r_  sin  v„n  pe  2 
•*  11 


n 


0,  1, 


(E-25) 


where  r\0  =  ~pj  c0  is  the  intrinsic  impedance  of  free  space.  It  should  be  mentioned  that  the 
above  result  differs  in  sign  from  the  one  given  by  Wait  (reference  1). 
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E.2  TE  MODES 


Consider  a  VMD  or,  equivalently,  an  infinitesimally  small  loop  of  electric  current  with  axis 
along  the  radial  direction  that  is  located  at  the  radial  distance  rs  =  a  +  zs  along  the  0  =  0  axis  of 

the  spherical  earth- ionosphere  waveguide  shown  in  figure  1-1.  As  is  shown  in  appendix  G,  the 
differential  equation  for  the  magnetic  Debye  potential  *Fh  is  given  by 

(V2  +  k 2)  Wh  =  -  fs  6(r  -  rs) ,  (E-26) 

where  m  is  the  current  moment  of  the  loop  and  8  denotes  the  Dirac  delta  function.  The  solution 
of  the  above  differential  equation  at  source-free  points  (r  ?  rs)  is  derived  in  appendix  C,  section 

C.2,  and  is  given  in  terms  of  a  modal  expansion  as 

Wh(r,6)  =  f  4  4  (w)  Pv  (-cos  6 )  ,  (E-27) 

nx~  \  m  m  m 

where  R1!  ( u )  is  defined  as 

vm 

4(«)s/i(i)(«)  +  Bj  hf(u)  .  (E-28) 

m  m  tn  m 

In  the  above  formulas,  m  denotes  the  mode  index,  Ah-  is  the  excitation  coefficient  of  the  mth  TE 

V  m 

mode,  and  B!1  is  a  constant  that  is  given  by  either  expression  (C-20a)  or  (C-20b).  The  complex 

i 

constant  vm  is  determined  from  the  solution  of  the  characteristic  equation  (C-21). 

From  appendix  D,  the  orthogonality  condition  for  the  TE  mode  radial  functions  is  given  by 
f'  K  (#>  *J  («)  (III  =  I  R»  I2  d„ ,  (E-29) 

where  ||  R*  j  denotes  the  L2  norm  of  as  defined  in  appendix  G  and  8mn  is  the  Kronecker 
delta  defined  in  expression  (E-5).  If  each  side  of  expression  (E-27)  is  multiplied  by  Ri  and 
then  integrated  with  respect  to  u,  the  application  of  the  orthogonality  condition  (E-29)  to  the 
resulting  equation  yields  an  expression  for  the  mode  excitation  coefficient  4  >  namely 

A! *  =  -p  1 -  P  Wh(r,6 )  R'l  (;/)  du  ,  m  =  1,  2,  3,  ...  .  (E-30) 

4  \\  P9  (-cos  6)  4 
II  V,n  II  1  m 

It  is  difficult  to  determine  4  from  equation  (E-30)  because  *PA(r,  9)  is  unknown. 
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To  determine  the  excitation  coefficients  A1!  in  formula  (E-30),  let  9->  0.  Then, 

V  m 


lim  f Wh(r,6 )  Rl  («)  du 
.  9-*oJu  m 

Ah  =  1 _ £ _ 

ph  2  lim  P„  (-  cos  6) 

Kvm  e-*o  v«  ’ 


(E-31) 


Following  the  same  procedure  as  for  the  TM  mode  excitation  coefficients,  the  magnetic  Debye 
potential  *F/l  satisfies  the  inhomogeneous  differential  equation  (E-26)  and  can  be  expressed  as 
the  sum  of  a  particular  solution  Wj?  and  a  homogeneous  solution  Wft .  The  term  Wjf  is  finite  for 
all  values  of  6  and  r  within  the  spherical  waveguide  and  is  singular  at  the  source  point  r  =  rs. 
The  particular  solution  W*  of  equation  (E-26)  is  the  same  as  the  magnetic  Debye  potential  for  a 
VMD  in  free  space  and  is  given  by  (reference  20) 

Vh(r,d)  =  Wsh(r, 6)  =  ,  (E-32) 

where  R  is  given  by  expression  (E-l  1). 


Because  the  integrand  in  the  numerator  of  formula  (E-31)  is  finite  everywhere  except  at  the 
source  point  (r  =  rs,  9=  0),  its  only  contribution  to  the  integral  is  from  the  immediate  neighbor¬ 
hood  of  the  source  point.  Therefore,  expression  (E-31)  can  be  simplified  to 


A 


h 


frs(  1  +  e) 

.  dA,  x  lim  Ws\r,9)dr 

koK  (Us)  e^°Jrs(  1-8) 


ph  2  lim  P„  (-  cos  9) 

Rvm  e-o  vm 


(E-33) 


where  u  has  been  replaced  by  k0r  and  e  «  1 .  Note  that  («)  has  been  replaced  by  R*l,  (u  ) 

’  m  ’  m  S 


because  it  is  finite  at  the  source  point.  If  the  procedure  described  in  appendix  E,  section  E.  1 ,  is 
followed,  the  integral  in  expression  (E-33)  reduces  to 


frs(  1  +  e)  m  „ 

lim  I  W?(r,0)dr  =  -- — In  92  .  (E-34) 

e^oJrsa-s)  s  4 xrs  v  7 


The  modal  excitation  coefficients  A1!  can  now  be  obtained  through  substitution  of  the 
limiting  results  (E-8)  and  (E-34)  into  expression  (E-33),  resulting  in 


E-7 


(E-35) 


kom  1 
4rs  sin 


R^O 
Ri  \2 


m  =  1,2,3, 


It  should  be  noted  that  the  above  result  differs  in  sign  from  the  one  given  by  Wait  (reference  1). 
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APPENDIX  F 

DERIVATION  OF  VERTICAL  ELECTRIC  DIPOLE  (VED) 
ELECTROMAGNETIC  (EM)  FIELDS 


In  this  appendix,  the  electromagnetic  (EM)  fields  radiated  by  a  vertical  electric  dipole  (VED) 
in  a  spherical  earth-ionosphere  waveguide  are  derived  through  the  use  of  the  spherical  wave 
formulas  derived  in  appendixes  A  and  C.  The  fields  can  be  expressed  in  terms  of  a  transverse 
magnetic  (TM)  mode  expansion. 

F.l  MODAL  EXPANSIONS  OF  EM  FIELDS 

Consider  a  radially -directed  electric  dipole  of  moment  p=Idl  that  is  located  at  the  radial 
distance  rs  =  a +zs  along  the  6  =  0  axis  of  the  spherical  earth-ionosphere  waveguide  shown  in 

figure  1  -1 .  The  electric  dipole  can  be  mathematically  expressed  in  terms  of  the  electric  current 
density  vector  J  as 


J  =  r  Jr  =  r  p  d(r  -  rs) , 


(F-l) 


where  r  is  the  unit  radial  vector,  r  and  rs  denote  the  position  vectors  of  the  observation  and 
source  points,  respectively,  and  S  denotes  the  Dirac  delta  function.  The  function  S(r  -  rs)  has  the 
SI  units  of  m'3  and  can  be  represented  in  spherical  coordinates  as 


d(r-rs) 


djr-r^  m 
2n r2  sin  6 


(F-2) 


The  spherical  coordinate  system  is  illustrated  in  figure  4-1 .  Note  that 

jj  j  d(r  -rs)dv=  1 ,  (F-3) 

where  the  integration  domain  extends  over  all  space. 


To  determine  the  expressions  for  the  EM  fields  radiated  by  the  VED,  one  must  first  obtain  the 
differential  equation  for  the  electric  Debye  potential  within  the  spherical  waveguide  (defined 

by  the  region  a<r<  a+h )  with  an  excitation  current  given  by  expression  (F-l).  From  equation 
(A-6),  appendix  A,  the  differential  equation  for  the  magnetic  vector  potential  A  is  given  by 

V  xV  xA-k20A  =  -  j(op0e0V  <Pe  +  pj,  (F-4) 


F-l 


where  fi0  is  the  permeability  of  free  space,  ea  is  the  permittivity  of  free  space,  k0  =  tosjii0£0  is 
the  wave  number  in  free  space,  and  @e  is  the  electric  scalar  potential.  Following  the  procedure 
shown  in  appendix  A,  section  A.l,  if  one  sets  A  =  r  Ar  and  applies  the  gauge  condition  given  in 
expression  (A- 18),  the  differential  equation  for  A  reduces  to 


(F-5) 


where  V2  is  the  Laplacian  operator.  It  should  be  mentioned  that  with  the  gauge  condition  (A- 18) 
chosen  and  with  only  a  radial  component  of  electric  current,  the  0  and  (p  components  vanish  on 
each  side  of  equation  (F-4).  From  the  definition  (A-26a),  the  magnetic  vector  potential  component 
Ar  can  be  expressed  in  terms  of  the  radial  component  of  the  electric  Hertz  vector  ner  as 

Ar  =  jl°fioEo7ier  •  (F-6) 


From  the  definition  (A-29a),  the  electric  Debye  potential  is  expressed  in  terms  of  ner  as 


(F-7) 


The  substitution  of  formulas  (F-6)  and  (F-7)  into  (F-5)  yields 


(V2  +  *o) 


IJ/e  = 


1  {r 
jW£o  r 


(F-8) 


For  the  VED  excitation  with  current  density  given  by  (F-l),  the  differential  equation  (F-8) 
becomes 


.  (F-9) 


It  should  be  noted  that  on  the  right  side  of  the  above  equation,  r  has  been  replaced  by  rs  because 
the  expression  vanishes  everywhere  except  at  r  —  rs.  The  solution  of  the  above  differential 
equation  at  source-free  points  (r  ?  rs)  is  derived  in  appendix  C,  section  C.  1,  and  is  given  in  terms 

of  a  modal  expansion  as 


We(r,d)  =  2  Aev  Rev  («)  Pv  (-  cos  6) ,  (F-10) 

,i  =  o  "  "  n 


where  n  denotes  the  mode  index,  Re  is  defined  in  formula  (C-10),  u  =  k0r,  Ae  is  the  excitation 
coefficient  of  the  nth  TM  mode,  Bev  is  a  constant  that  is  given  by  either  expression  (C-7a)  or  (C- 


F-2 


7b),  hSy\u)  and  h^\u)  denote  the  spherical  Hankel  functions  of  the  first  and  second  kinds 

n  n 

respectively,  of  order  v„,  and  P  denotes  the  Legendre  function  of  the  first  kind  of  degree  v„ 
and  order  zero.  The  complex  constant  vn  is  determined  from  the  solution  of  the  characteristic 
equation  (C-8).  The  modal  excitation  coefficients  Ae  depend  on  the  source  strength p/coe0rs. 
The  term  Aev  is  derived  in  appendix  E,  section  E.  1 ,  and  is  given  by 


Ae  -MoP 

Av„~  4 r 


1 


S  sin  vn% 


K  («,) 

,  n  =0,  1,2,..., 

K 


(F-ll) 


where  T)0  =  ^fi  /e  is  the  intrinsic  impedance  of  the  free  space.  In  expression  (F-l  1), 
denotes  the  L2  norm  of  Rev  and  is  defined  as 


du  ,  n =  0, 1, 2, ... , 


(F-l 2) 


where  ug  =  k0a  and  ut  =  k0(a+h). 


The  EM  fields  radiated  by  a  VED  are  given  in  formulas  (C-l  1)  with  modal  excitation 
coefficients  given  by  expression  (F-l  1).  Therefore,  the  VED  EM  field  components  in  spherical 
coordinates  are  given  as 


Fve  _  jam0p 

r  "  4  rsu 


y  vw(vB+ !)  KW 
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K(us)  a 


n  =  0  Sin  VMJl  II  Re 
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(F-13a) 


(F-l  3b) 


and 


I  00 

Hve=KP  y  1 

V  4rs  n  =  0  sin 


^n(M)MFv"(-COS  ^ 


In  the  above  formulas,  the  superscript  ve  refers  to  a  VED  source. 


(F-l  3c) 


Following  Wait  (reference  1),  the  effectiveness  of  a  source  can  be  characterized  by  an 
excitation  factor  Aen  and  the  height  variation  of  the  fields  by  the  height-gain  function  Gen(z)  ■ 
These  quantities  are  defined  as  follows: 


F-3 


,n  =0,1,2 


(F-14) 


and 


n  =0, 1,2,..., 


(F-15) 


where  z  =  r  —  a  is  the  radial  distance  measured  from  the  surface  of  the  spherical  earth  (see  figure 
1-1).  Following  Galejs  (reference  2),  a  z-dependent  normalized  impedance  is  defined  as 


A%z)  - 


*1 0 
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mode  n 


uRev  (;/) 


n  =0,1,2,.... 


(F-16) 


The  application  of  the  definitions  (F-14)  to  (F-16)  to  the  VED  formulas  (F-13)  yields 


Ke  =  2  V,,(V,i  +  ~  K  Ge„(zs)  Gen (z)  Pv  (-  cos  0) , 

r  2 hrsu  ,fr0  sin  vHn  "  "  s  "  v* 


Tve  _  VoP 


2hrs  »  =  o  sin  v„ji 


1  A* 4,fc)  G'(zs)  Gen{z )  ~  Pv(- cos  0) , 


and 


H*  =  tA-  2  -v-1 —  Ae„  GeJzs)  Gf.(z)  4a  pv  (-  cos  0) 
V  2 rsh  ,“0  sin  v,(ir  "  "  5  "  dd  v"  7 


(F-17a) 

(F-17b) 

(F-17c) 


The  above  expressions  are  similar  to  those  given  by  Wait  (reference  1)  and  Galejs  (reference  2). 
The  radial  variations  in  the  formulas  for  each  of  the  above  field  components  may  be 
approximated  in  order  to  make  them  suitable  for  numerical  computation. 


If  the  VED  and  the  field  point  are  each  located  on  the  inner  boundary  of  the  waveguide 
(earth’s  surface),  then zs  =  z  =  0  or,  equivalently,  rs  =  r  =  a.  Thus, 

Ge„(zs)  =  Ge„(z)  =  G% 0)  =  1  ,  n  =  0,  1, 2, ... .  (F-18) 

In  addition,  from  equation  (C-6a),  the  z-dependent  normalized  impedance  AeJt{z)  defined  in 
expression  (F-16)  reduces  to 


F-4 


Ae„(z)  =  Aen{0)  =  -Ag,  n  =0,1,2 


(F-19) 


where  Ag  =  Z^r\Q  denotes  the  normalized  surface  impedance  of  the  inner  boundary.  The 
substitution  of  the  above  simplifications  into  the  VED  field  expressions  (F-17)  gives 


E>^Zv-^AlP^e)- 

(F-20a) 

E0  --  2 ha  1  sin  v„Jt  Afl  36  PvS~  cos  6)  ’ 

(F-20b) 

and 

Hv*=p  2  1  K  za  Pv  (-cos  6)  . 

v  2 ha  tfr*0  sin  v„ji  n  36  v«  ’ 

(F-20c) 

F.2  QUASI-TEM  FIELDS 

The  tables  of  section  4.1  show  that  in  the  ELF  band,  all  modes  are  nonpropagating  except  the 
n  —  0  mode.  This  mode  is  often  referred  to  as  the  quasi-transverse  EM,  or  quasi-TEM,  mode.  If 
only  the  quasi-TEM  mode  is  considered,  the  VED  fields  for  a  source  and  receiver  each  located  at 
the  surface  of  the  inner  sphere  (earth)  are  given  as 

Ev,  Vv.+  1) 

r  2 haug  sin  v0n  °  V  7 

(F-21a) 

£ve  _  60PAg  1  9  p  (  COo  Q) 

6  2 ha  sin  v0n  °  36  v<r  '  ’ 

(F-21b) 

and 

HVJ=P  .  1  Ae  ^  P  (-  cos  6)  . 
v  2 ha  sin  v  re  0  36  V  ' 

(F-21c) 

The  above  field  expressions  are  of  primary  interest  in  this  report. 
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APPENDIX  G 

DERIVATION  OF  VERTICAL  MAGNETIC  DIPOLE  (VMD) 
ELECTROMAGNETIC  (EM)  FIELDS 


In  this  appendix,  the  electromagnetic  (EM)  fields  radiated  by  a  vertical  magnetic  dipole 
(VMD)  in  a  spherical  earth-ionosphere  waveguide  are  derived  through  the  use  of  the  formulas 
derived  in  appendixes  A  and  C.  The  fields  can  be  expressed  in  terms  of  a  TE  mode  expansion. 

Consider  a  radially -directed  magnetic  dipole  of  moment  K  dl  that  is  located  at  the  radial 
distance  rs  =  a  +  zs  along  the  6  =  0  axis  of  the  spherical  earth-ionosphere  waveguide  shown  in 

figure  1  -1 .  The  magnetic  dipole  can  be  mathematically  expressed  in  terms  of  the  magnetic 
current  density  vector  M  as 

M  =  r  Mr  =  r  Kdl  d(r - rs)  ,  (G-l) 

where  r  is  the  unit  radial  vector,  r  and  rs  denote  the  position  vectors  of  the  observation  and 
source  points,  respectively,  K  denotes  the  magnetic  current,  and  8  denotes  the  Dirac  delta 
function.  Because  an  infinitesimal  dipole  of  magnetic  current  is  equivalent  to  an  infinitesimal 
loop  of  electric  current  (where  the  loop  and  dipole  axes  are  parallel)  through  the  relation  Kdl- 
j(OHo  Ida  (references  19  and  20),  where  I  is  the  loop  current  and  da  is  the  differential  area 

enclosed  by  the  loop,  equation  (G-l)  may  also  be  written  as 

M  =  r  joifxjn  8(r  -  rs) ,  (G-2) 

where  m  =  I  da  is  the  current  moment  of  the  loop. 

To  determine  the  expressions  for  the  EM  field  radiated  by  the  VMD,  one  first  needs  to  obtain 
the  differential  equation  for  the  magnetic  Debye  potential  *Fh  within  the  spherical  waveguide 
(defined  by  the  region  a<r<  a+h)  with  an  excitation  current  given  by  expression  (G-2).  From 
equation  (A-l  1),  the  differential  equation  for  the  electric  vector  potential  F  is  given  by 

VxVxf-  k2QF  =  -  ja>ii0£01V<Ph  +  e0M ,  (G-3) 


where  <Ph  is  the  magnetic  scalar  potential.  Following  the  procedure  shown  in  appendix  A, 
section  A.l,  by  setting  F  =  rFr  and  applying  the  gauge  condition  (A-24),  the  differential 
equation  for  F  reduces  to 


G-l 


(G-4) 


It  should  be  mentioned  that  with  the  gauge  condition  (A-24)  chosen  and  with  only  a  radial 

component  of  magnetic  current,  the  0  and  (p  components  vanish  on  each  side  of  equation  (G-3). 
From  the  definition  (A-26b),  the  electric  vector  potential  component  Fr  can  be  expressed  in 

terms  of  the  radial  component  of  the  magnetic  Hertz  vector  nhr  as 

Fr  =  j(0t*0£oJtr  •  (G-5) 

From  the  definition  (A-29b),  the  magnetic  Debye  potential  'f/h  is  expressed  in  terms  of  nhr  as 

Wh  =  ^ .  (G-6) 


The  substitution  of  formulas  (G-5)  and  (G-6)  into  (G-4)  yields 


V2  +  £2)  wh 


1  Mr 
V 


(G-7) 


For  the  VMD  excitation  with  current  density  given  by  expression  (G-2),  the  differential 
equation  (G-7)  becomes 

(V2  +  *2)  tyh  =-y-  <5 (r  -  rs)  .  (G-8) 


It  should  be  noted  that  on  the  right  side  of  equation  (G-8),  r  has  been  replaced  by  rs  because  the 
expression  vanishes  everywhere  except  at  r  —  rs.  The  solution  of  the  above  differential  equation 
at  source-free  points  (r  ?  rs)  is  derived  in  appendix  C,  section  C.2,  and  is  given  in  terms  of  a 

modal  expansion  as 


Vh(r,8)=  2  A% 

m  =  1 


RU")Py(-  cos  6), 

*  m  m 


(G-9) 


where  m  denotes  the  mode  index,  R1!  is  defined  in  formula  (C-23),  Ah-  is  the  excitation 

'  m  Vm 

coefficient  of  the  mth  TE  mode,  B1!  is  a  constant  that  is  given  by  either  expression  (C-20a)  or 

*  m 

(C-20b),  and  the  complex  constant  vm  is  determined  from  the  solution  of  the  characteristic 
equation  (C -21).  The  modal  excitation  coefficients  Ah-  depend  on  the  source  strength  m/rs.  The 
term  A*  is  derived  in  appendix  E,  section  E.2,  and  is  given  by 
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In  the  above  expression,  i?*  denotes  the  L2  norm  of  R%  and  is  defined  as 


Rh  = 


f 


du  , 


(G-ll) 


where  u„  =  k0a  and  «,•  —  k 0(a+h). 


The  EM  fields  radiated  by  a  VMD  are  given  in  formulas  (C-24)  with  modal  excitation 
coefficients  given  by  expression  (G-10).  Therefore,  the  VMD  EM  field  components  in  spherical 
coordinates  are  given  as 
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(G-12a) 


(G-12b) 


(G-12c) 


In  the  formulas  (G-12),  the  superscript  vm  refers  to  a  VMD  source.  Following  Wait  (reference  1) 
and  Galejs  (reference  2),  the  excitation  factor  A*  ,  height-gain  function  G*(z) ,  andz-dependent 
normalized  impedance  A*  (z)  for  TE  modes  are  defined  as 
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(G-15) 


The  application  of  the  definitions  (G-13)  to  (G-15)  to  the  VMD  formulas  (G-12)  yields 
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(G-16a) 

(G-16b) 

(G-16c) 


The  above  expressions  are  similar  to  those  given  by  Wait  (reference  1)  and  Galejs  (reference  2). 
The  radial  variations  in  the  formulas  for  each  of  the  above  field  components  may  be 
approximated  in  order  to  make  them  suitable  for  numerical  computation. 

If  the  VMD  and  the  field  point  are  each  located  on  the  inner  boundary  of  the  waveguide 
(earth),  then zs  =  z  =  0  or,  equivalently,  rs-  r  =  a.  Thus, 


Gi(zs)  =  Ghm{z)  =  Ghm( 0)  =  1  ,  m  =  1,  2, 3, ... , 
and  the  z-dependent  normalized  impedance  A^n (c)  reduces  to 
At(z)  =  A>l,(0)=-Ap,  m  =1,2,3,.... 


(G-17) 


(G-18) 


The  substitution  of  the  above  simplifications  into  the  VMD  field  expressions  (G-16)  gives 
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The  above  field  expressions  are  of  primary  interest  in  this  report. 
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APPENDIX  H 

DERIVATION  OF  HORIZONTAL  ELECTRIC  DIPOLE  (HED) 
ELECTROMAGNETIC  (EM)  FIELDS 

In  this  appendix,  the  electromagnetic  (EM)  fields  radiated  by  a  horizontal  electric  dipole 
(HED)  in  a  spherical  earth-ionosphere  waveguide  are  derived  from  the  vertical  electric  dipole 
(VED)  and  vertical  magnetic  dipole  (VMD)  fields  through  the  use  of  the  reciprocity  theorem. 

The  expressions  for  the  nonradial  field  components  involve  both  TM  and  TE  mode  expansions. 

H.l  MODAL  EXPANSIONS  OF  EM  FIELDS 

The  Lorentz  reciprocity  theorem  (references  19  and  20)  relates  a  response  at  one  source  due 
to  a  second  source  to  the  response  at  the  second  source  due  to  the  first  source.  Consider  two  sets 
of  impressed  line-current  sources  (/ h  K{)  and  (/2,  K2)  excited  at  the  same  frequency  f  where  I 
and  K  refer  to  electric  and  magnetic  line  currents,  respectively.  The  EM  fields  produced  by  (7j, 
Kx)  are  (Ei}  H{)  and  the  fields  produced  by  (/2,  K2)  are  (E2,  H2).  If  the  line  currents  lie  in  a 
linear,  reciprocal  medium,  the  reciprocity  theorem  applies  and  is  given  as 

/  (e,  •  l2dl2  -  H,  •  Ejd/j)  =  f  (e2  •  /,<«;  -  «2  •  *,<«,) ,  (H-l) 

where  dlx  and  dl2  denote  the  differential  length  vectors  along  line  current  sets  one  and  two, 
respectively. 

To  derive  the  fields  produced  by  a  HED,  consider  the  two  coordinate  systems  defined  in 
figure  H- 1 .  In  this  illustration,  the  coordinates  (x,  y,  z)  refer  to  those  of  the  HED,  and  the  primed 
coordinates  (x\y\ z')  refer  to  those  of  the  VED  (case  1)  and  VMD  (case  2).  The  origin  of  each 
coordinate  system  lies  at  a  point  on  the  earth’s  surface,  with  the  x-y  and  x-j/ planes  each  lying  on 
the  tangent  plane  defined  at  their  respective  origins.  The  distance  a#  is  the  great-circle  path  length 
between  the  origins  of  the  two  coordinate  systems,  where  6  is  the  polar  angle  defined  in  figure 
1-1.  The  reciprocity  theorem  is  applied  for  two  different  cases.  In  case  1,  anx-directed  HED  with 
electric  dipole  moment  Ix  dlx  is  located  at  coordinates  (x,  y,  z)  =  (0,  0,  zr)  and  a  z -directed  VED 
with  electric  dipole  moment  I2  dl2  is  located  at  coordinates  (x',y\z')  =  (0,  0,  zs).  In  case  2,  an 
x-directed  HED  with  electric  dipole  moment  7j  dlx  is  located  at  coordinates  (x,  y,  z)  =  (0,  0,  zr ) 
and  a  z -directed  VMD  with  magnetic  dipole  moment  K2  dl2  is  located  at  coordinates  (x',y',z')  = 


H-l 


(0,  0,  zs).  The  radial  component  of  electric  field  E*e  produced  by  the  HED  is  determined  from 
reciprocity  in  case  1  and  the  radial  component  of  magnetic  field  Hhrc  produced  by  the  HED  is 
obtained  from  reciprocity  in  case  2,  where  the  superscript  he  refers  to  a  HED  source. 
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Figure  H-l.  Coordinate  Systems  Defined  in  the  Application  of  the  Reciprocity 
Theorem  for  Determination  of  the  HED  Fields 


To  determine  Ehre ,  consider  case  1  in  figure  H-l .  The  application  of  the  reciprocity  theorem 
(H-l)  to  the  HED  and  VED  sources  yields 


E'ltix’  =  0,  /  =  0,  z'  =  zs)  1 2  dl2  =  E\cx{x  =  0,y  =  0,z  =  zr)  f  df  , 


(H-2) 


where  E*e.,(x'  =  0,  y'  =  0,  z'  =  zs)  is  the  vertical  (radial)  component  of  electric  field  from  the  HED 
evaluated  at  the  VED  location  and  E2x(x  =  0,  y  =  0,  z  =  Zr)  is  the  x  component  of  electric  field 
from  the  VED  evaluated  at  the  HED  location.  Setting  I1  dl\  —I2  dl2,  expression  (H-2)  reduces  to 


El l(x'  =  0,  y  =  0,  z!  -  Zs)  =  E™(x  =  0,  y  =  0,z  =  Zr) 


(H-3) 


H-2 


The  left  side  of  equation  (H-3)  can  be  expressed  in  spherical  coordinates  as  £^(r  =  rs,  6,  (p) , 
where  the  spherical  coordinates  of  the  source  point  are  ( r ,  6)  =  (rr,  0)  and  the  spherical 
coordinates  associated  with  the  field  point  are  (rs,  0,  (p).  Similarly,  the  right  side  of  equation 
(H-3)  can  be  expressed  in  spherical  coordinates  as  -  Ev2&{r  =  rr,  &)  cos  cp ,  where  the  spherical 
coordinates  of  the  source  point  are  (r,  O')  =  ( rs ,  0)  and  the  spherical  coordinates  associated  with 
the  field  point  are  ( rs ,  6',  (pA),  where  O'-  6 and  6’  =  - 6.  The  projection  of  Ev2e.  onto  the  x  axis 
is  illustrated  in  figure  H-2.  Note  that  the  azimuthal  dependence  is  missing  as  an  argument  of 
£20'  because  the  VED  fields  are  axisymmetric. 
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Figure  H-2.  Electric  Field  Component  E2&  Radiated  by  a  VED  ( Source  2) 
Projected  onto  the  x-Axis  at  the  HED  ( Source  1 )  Location 


If  the  above  spherical  coordinate  substitutions  are  applied  to  equation  (H-3),  then 
E?(rs,  6,  <p)  =  -  Evee,(rr,  6 ')  cos  cp ,  (H-4) 


where  the  source  numbers  have  been  dropped.  With  the  radial  coordinate  of  the  VED  given  by 
r-rs  and  the  coordinates  of  the  field  point  given  by  ( rr,6 from  formula  (F-15),  EVq,  is 
expressed  as 
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cos  &) ,  (H-5) 


where  ur  =  k0rr  The  substitution  of  expression  (H-5)  into  (H-4)  yields 
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The  above  formula  applies  for  an  ^-directed  HED  located  at  ( r ,  0)  =  (rr,  0)  and  the  field  point  at 
(rs,  0,  (p).  If  the  HED  is  located  at  (r,  0)  =  (rs,  0)  and  the  field  point  coordinates  are  (r,  0,  <p),  then 
rs  is  replaced  by  r  and  rr  is  replaced  by  rs  in  expression  (H-6),  resulting  in 


Ehr\r,  0,  cp)  =  J--^y  cos  cp 
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where  0'-  0and  dld0'=  -  d/d0 . 


To  determine  Hze ,  consider  case  2  in  figure  H-l.  The  application  of  the  reciprocity  theorem 
(H-l)  to  the  HED  and  VMD  sources  yields 


H^(x’  =  0,  y’  =  0,  z'  =  zs)=-  E™(x  =  0,y  =  0,z  =  zr) 


Elm' 


K2  dl2 


(H-8) 


where  H*ez,(x'  =  0 ,y'  =  0,  z'  =  zs)  is  the  vertical  (radial)  component  of  magnetic  field  from  the 
HED  evaluated  at  the  VMD  location  and  E™(x  =  0,  y  =  0,  z  =  Zr)  is  the  x  component  of  electric 
field  from  the  VMD  evaluated  at  the  HED  location.  The  left  side  of  equation  (H-8)  can  be 
expressed  in  spherical  coordinates  as  H*er(r  =  rs,  6,  (p) ,  where  the  spherical  coordinates  of  the 
source  point  are  (r,  0)  =  (rr,  0)  and  the  spherical  coordinates  associated  with  the  field  point  are  (rs, 
0,  <p).  Similarly,  the  right  side  of  equation  (H-3)  can  be  expressed  in  spherical  coordinates  as 
E™(r  =  rr,  6')  sin  cp ,  where  the  spherical  coordinates  of  the  source  point  are  (r,  0')  -  (rs,  0)  and 
the  spherical  coordinates  associated  with  the  field  point  are  ( rs ,  0'  ip'),  where  0'-  (9  and  q>’  =  -  q> . 
The  projection  of  E2“,  onto  the  x  axis  is  illustrated  in  figure  H-3.  Note  that  the  azimuthal 
dependence  is  missing  as  an  argument  of  E2“,  because  the  VMD  fields  are  axisymmetric. 
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Figure  H-3.  Electric  Field  Component  E2™  Radiated  by  a 
VMD  (Source  2)  Projected  onto  the  x-Axis  at  the 
HED  (Source  1)  Location 


If  the  above  spherical  coordinate  substitutions  are  applied  to  equation  (H-8),  then 


(H-9) 


where  the  magnetic  dipole  has  been  replaced  by  the  equivalent  infinitesimal  electric  current  loop 
of  moment  l2  da2 ,  where  K2  dl2  =  jo)ju0 12  da2.  With  the  radial  coordinate  of  the  VMD  given  by 
r=rs  and  the  coordinates  of  the  field  point  given  by  (rr,  O'),  from  formula  (G-12c),  ZT"".  is 
expressed  as 
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The  substitution  of  expression  (H-10)  into  (H-9)  yields 
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where  p1  =  Ix  dlv  The  above  formula  applies  for  an  ^-directed  HED  located  at  (r,  6)  =  (rr,  0) 
and  the  field  point  at  ( rs ,  6,  <p).  If  the  HED  is  located  at  (r,  6)  -  (rs,  0)  and  the  field  point 
coordinates  are  (r,  6 \  (p),  then  rs  is  replaced  by  r  and  rr  is  replaced  by  rs  in  expression  (H-l  1), 
resulting  in 
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where  6'=  dznd  did  6'=  -  dldd,  and  the  subscript  1  has  been  removed. 


The  remaining  field  components  for  a  HED  source  can  be  determined  from  the  Debye 
potentials  as  given  by  equations  (A-31).  Because  the  Debye  potentials  for  a  HED  have  not  yet 
been  determined,  they  must  be  derived  from  the  radial  field  components  given  by  expressions 
(H-7)  and  (H-12).  From  equations  (A-31a)  and  (A-31d),  the  radial  field  components  are 
expressed  in  terms  of  the  Debye  potentials  as 
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(H-l  3a) 


(H-l  3b) 


where  We  and  *Fh  denote  the  electric  and  magnetic  Debye  potentials,  respectively.  Given  the 
expressions  for  Ehre  and  H*e  in  formulas  (H-7)  and  (H-12),  respectively,  suitable  forms  for  We 
and  5P*  are  as  follows: 
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where  Aev  and  A*  denote  the  TM  and  TE  mode  expansion  coefficients,  respectively. 


From  the  relation  (C-12), 
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The  substitution  of  expressions  (H-15a)  and  (H-15b)  into  (H-14a)  and  (H-14b),  respectively, 
yields 
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A  comparison  of  formulas  (H-7)  and  (H-16a)  shows  that 
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Similarly,  a  comparison  of  formulas  (H-12)  and  (H-16b)  yields 
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The  substitution  of  expressions  (H-17a)  and  (H-17b)  into  (H-14a)  and  (H-14b),  respectively, 
yields 
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From  formulas  (A-31),  the  remaining  HED  field  components  are  given  in  terms  of  the  Debye 
potentials  as 
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The  substitution  of  the  Debye  potential  expressions  (H-18)  into  the  HED  field  formulas  given 
above  yields 
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Formulas  (H-20)  show  that,  unlike  the  VED  and  VMD  fields,  four  of  the  six  HED  field 
components  involve  the  superposition  of  both  TM  and  TE  modes. 

Following  Wait  (reference  1)  and  Galejs  (reference  2),  the  HED  formulas  can  be  expressed  in 
terms  of  the  excitation  factor,  height-gain  function,  and  normalized  impedance  defined  in 
expressions  (F-14)  to  (F-16)  and  (G-13)  to  (G-15).  The  application  of  these  functions  to  the 
HED  formulas  (H-7),  (H-12),  and  (H-20)  yields 
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The  above  expressions  are  similar  to  those  given  by  Wait  (reference  1)  and  Galejs  (reference  2). 
The  radial  variations  in  the  formulas  for  each  of  the  above  field  components  may  be 
approximated  in  order  to  make  them  suitable  for  numerical  computation. 

If  the  HED  and  the  field  point  are  each  located  on  the  earth’s  surface,  then  zs 
Therefore,  the  height-gain  functions  reduce  to 
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In  addition,  the  z-dependent  normalized  impedances  A*(z)  and  Ajjt(z)  defined  in  expressions 
(F-16)  and  (G-15),  respectively,  reduce  to 
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The  substitution  of  the  above  simplifications  into  the  HED  formulas  (H-21)  gives  the  following: 
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H.2  QUASI-TEM  FIELDS 


The  tables  of  section  4.1  showed  that  in  the  ELF  band,  all  modes  are  nonpropagating  except 
the  n  =  0  TM  mode.  This  mode  is  often  referred  to  as  the  quasi-transverse  EM,  or  quasi- TEM 
mode.  If  only  the  quasi- TEM  mode  is  considered,  the  HED  fields  for  a  source  and  receiver  each 
located  at  the  surface  of  the  earth  are  given  as 
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?  2/j  v0(v0  +  1)  sin  v0Jt  sin  6  dd  '<> 


„he_jko\PK  1 
"»  "  2A 


1  1  0 

TTfl  Pv  (“  cos  e)  sin  <P  . 


vo(vo+  1)  sin  v03i  sin  6  08  v« 


and 


JKAgPK  i  i 

^  2/i  v0(v0  +  1)  sin  v0ji  36 


2  Fv  (-  cos  0)  COS  (p  . 


(H-25a) 


(H-25b) 


(H-25c) 


(H-25d) 


(H-25e) 


Note  that  the  radial  magnetic  field  component  H^e  is  not  listed  above  because  it  is  expressed 
only  in  terms  of  TE  modes. 
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APPENDIX  I 

THIN-SHELL  APPROXIMATION  TO  THE 
RADIAL  DEPENDENCE  OF  THE  FIELDS 


The  spherical  waveguide  formulas  derived  in  appendices  F,  G,  and  H  involve  the  Legendre 
function  of  the  first  kind  for  the  6  dependence  and  the  spherical  Bessel  functions  for  the  r  (radial) 
dependence,  where  each  function  is  of  complex  order.  Exact  and  approximate  representations 
for  the  Legendre  function  are  presented  in  appendices  J  and  K.  Galejs  (reference  2)  has 
presented  an  approximate  formula  for  the  radial  variation  of  the  fields  under  the  condition  hlr  « 

1.  This  approximation,  referred  to  as  the  thin-shell  approximation,  is  suitable  in  the  extremely 
low  frequency  (ELF)  band  and  provides  results  of  qualitative  accuracy  in  the  very  low  frequency 
(VLF)  band.  In  this  appendix,  a  derivation  of  the  thin-shell  approximation  is  presented  and 
applied  to  obtain  approximate  solutions  of  the  transverse  magnetic  (TM)  and  transverse  electric 
(TE)  mode  characteristic  equations.  Thin-shell  approximations  for  the  mode  excitation  factors, 
the  height-gain  functions,  and  the  z-dependent  normalized  impedances  are  also  presented. 

1.1  APPROXIMATE  SOLUTION  OF  THE  RADIAL  DIFFERENTIAL  EQUATION 


In  appendix  B,  it  is  shown  that  the  radial  dependence  of  each  Debye  potential  is  given  in 
terms  of  a  function  Rv(u),  u  =  k0r,  which  satisfies  the  following  differential  equation: 


d_ 

dr 


dR\ 

r  i 

KJ+ 

(k0r)2-v(v+ 1) 

Rv  =  0 


(1-1) 


As  shown  in  appendix  B,  two  linearly  independent  solutions  of  the  above  differential  equation 
are  given  by  the  spherical  Hankel  functions  of  the  first  and  second  kinds,  h(y\u)  and  h^\u) , 

respectively.  Consider  the  following  function: 

Wv{r)  =  r  Rv(k0r)  .  (1-2) 


The  substitution  of  the  above  function  into  equation  (1-1)  yields  the  following  differential 
equation  for  Wy(r)  : 


Wy(r)  + 


,.2  V(  V  +1) 
Ko  r2 


Wv(r)  =  0  . 


(1-3) 
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For  r  ->  8  or  for  fields  in  a  thin  shell  of  h/r«  1 ,  the  radial  variable  r  that  appears  in  brackets 
in  the  above  equation  may  be  approximated  by  its  average  value  r  in  the  shell  that  is  given  by 

r  s  r  =  a +  e,  (1-4) 


where  e  =  h  la  «  1 .  The  application  of  the  above  approximation  into  the  differential  equation  (I- 
3)  yields 


W';(r)  +  k2n,Wv(r)e  0, 
where 


k2  _  /,2  v(v±l) 

Krv~Ko  $2 


a-5) 

a-6) 


If  |  V  |  »  1,  then  from  formulas  (K-26)  and  (K-27), 
v(v +1)  s  ( ka )2  , 


a-7) 


where  A:  is  the  wave  number  in  the  6  direction.  From  the  approximations  (1-4)  and  (1-7),  krv  in 
formula  (1-6)  is  approximated  as 

k2 


b  ^  c*  b  2  _ 


—  k2  -  k 2|  1  -  ^ j  a  k2  -  k2  ,  e/a  «  1  . 


d-8) 


From  the  final  expression  in  the  above  approximation,  krv  can  be  interpreted  as  the  wave  number 
in  the  radial  direction.  Also  note  that  the  above  approximation  slightly  underestimates  the  radial 
wave  number  for  all  values  of  r. 


The  solution  of  equation  (1-5)  is  given  by 

Wv(r)  =  A  ejkrS  +  B  e~Jk"r ,  (1-9) 

where  A  and  B  are  arbitrary  constants.  From  definition  (1-2),  Rv  is  approximated  as 


WJr)  A  eJKs  .  „  e-JkrS 


Rv  m  =  A 

v  r  r 


+  B 


a-io) 


A  comparison  of  the  exact  solution  (B-7)  for  Rv  with  the  approximation  given  above  yields 

eJKrS 


h[ }\u)  s 


(I- Ha) 


and 
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,  ™  e~jK”r 
h(^(u)  m 


a-nb) 


The  above  approximations  are  used  to  solve  the  TM  and  TE  mode  characteristic  equations  that 
are  derived  in  appendix  C. 

1.2  APPLICATION  OF  THIN-SHELL  APPROXIMATION  TO  TM  MODES 


Consider  the  characteristic  equation  for  the  TM  modes  given  in  expression  (C-8).  The 

application  of  the  thin-shell  approximation  (1-1  la)  to  the  derivative  terms  in  formula  (C-8)  that 
involve  h^\u)  yields  the  following: 


and 


=  h^\u)  +  u  A  J>'(«)  =  (l  +  jknr) 


(M2) 


(M3) 


Following  a  similar  procedure  for  hS^\u) ,  from  approximation  (1-1  lb),  gives 


(M4) 


The  substitutions  of  the  thin-shell  approximations  (1-1 1),  (1-13),  and  (1-14)  into  the  quotients  in 
equation  (C-8)  yields 


e-j2krva 


and 


AJ k\ 


krv+  "g-'o 
krv-Agkoj 


,-jlka 


(I- 15a) 
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1  <\uh^Xu) 
ui  du 


+  jAeh^Xut) 


a  -  u. 


1  d[«/i(v2)(w)] 

ui  dii 


+  jAJi^X„i) 


«HJ, J  +  #,v+W 


1 


(a  +  h ) 


-j(krv-Aek0) 


J2kja+h) 


^  ^  rv  +  A<M0\  j  2k  rv(a+h) 
\krv-Acko) 


(I-15b) 


The  final  forms  of  approximations  (I- 15a)  and  (I -15b)  are  suitable  because  \ia  «  1.  The 
substitution  of  the  above  quotients  into  the  TM  mode  characteristic  equation  (C-8)  yields 


K%+v»'' 

lkrvn  +  AJ<o\ 

-v.j 

J2krv  h 
£  n 


1  =  ej2,,JT 


n  =  0,  1,2,.... 


(1-16) 


In  the  above  equation,  the  index  n  corresponds  to  the  order  of  the  TM  mode.  Note  that  there  are 
an  infinite  number  of  solutions  to  the  above  equation. 


The  TM  mode  characteristic  equation  (1-16)  may  be  rewritten  as 


1  +  Kf 


1  +  -"'  =e~j WrVk-»x)  ,  „  =  o,  1,2,..., 


1C. 


(1-17) 


where  Kn  and  are  defined  as 


JK 

HI 

>  &o 

3  O 

(1-18) 

f  _Aeko 
b 

(W9) 

In  the  above  expressions,  note  that  a  subscript  has  been  added  to  the  radial  wave  number  to 
indicate  the  mode  index.  The  natural  logarithm  of  equation  (I- 1 7)  yields 


In  the  ELF  band,  the  spherical  waveguide  boundaries  appear  as  nearly  perfect  electric 
conductors.  Therefore,  «  1  and  \Ae\  «  1  and  the  following  conditions  are  valid: 


K„  «  1.  «  1  • 


a-21) 


From  reference  27, 


wH)=2hf+f+ 


^  J ,  IjtI  «  1 


(1-22) 


The  application  of  the  approximations  (1-21)  and  (1-22)  into  the  TM  mode  characteristic  equation 
(1-20)  results  in 


Ko  j*3 

Kn  +  ~f  +  £n  +  ^  s  “  j[krvnh  ~  njT)  >  »  =  0,  1,  2,  ... . 


(1-23) 


The  substitution  of  the  variable  definitions  (1-18)  and  (1-19)  into  the  above  equation  yields  the 
following  quartic  equation  for  krv  : 


1 4  _  tin  j.3 

rVn  h  rVn 


j(A  +Ae)k0 


h 


k2 

rv„ 


j(A3g  +  A3e)k30 


3  h 


=  0  ,  n  =  0,  1,2, 


(1-24) 


Approximate  solutions  of  the  above  equation  for  the  quasi-TEM  mode  (n  =  0)  will  be  determined 
separately  from  the  higher-order  TM  modes. 


From  the  parallel-plate  waveguide  results  given  in  section  2,  it  is  known  that  the  quasi-TEM 

mode  (n  =  0)  is  the  only  propagating  mode  over  most  of  the  ELF  band.  If  n  =  0,  the  quartic 
equation  for  k  given  above  reduces  to  the  following  quadratic  equation  for  k2  : 

'  1  n  TV  n 


k4  - 

rv„ 


j(Ag  +  Ae)k0 

lr2 

'j(A3g  +  A3e)k30 

h 

krv0 

3  h 

=  0 


a-  25) 


A  solution  of  the  above  equation  is 

,  mji  wf+4?) 

A  '-k  .  .  r 


2  _Mg  +  Ae>ko 


kz  = 

Krv 


2 h 


1  + 


3 


1/2  i 


(1-26) 
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Because  the  earth  is  more  highly  conducting  than  the  ionosphere,  the  following  simplification  is 
allowed: 


M3  +  43)  3AAe 

— - - —  =  A  4 A _ -  «  a  -  'M 

( As  +  Ae )2  As  +  A*  Ag  +  Ar  -A‘  ~V 


g  e 


« A 


(1-27) 


The  substitution  of  the  above  approximation  to  the  solution  (1-26)  yields 


k2  — 

n' 


j(Ag  +  Ae)k0 

2 h 


1  + 


j4k  h 

\-J-^-(Ae-2Ag) 


j(Ag  +  Ae)k0 
h 


Aek0h  «1 


(1-28) 


The  above  result  corresponds  with  expression  (87)  on  page  89  of  Galejs  (reference  2).  Note  that 
as  A  -*  0  and  A  -*  0 ,  corresponding  to  perfectly  conducting  boundaries,  then  k  ,  -*■  0  and 

S  r'o 

k&kQ,  resulting  in  a  purely  TEM  mode  as  described  in  section  2. 


For  the  higher-order  TM  modes  (//  *  0 ),  if  only  the  first-order  term  in  the  series  expansion  (I- 

22)  is  applied  to  the  TM  mode  characteristic  equation  (1-20),  the  following  quadratic  equation  for 
krv  is  obtained: 


k2  _  ML  k  - 

rvn  h  rvn 


j(A  +  Ae)k0 


h 


=  0  ,  n=  1,2,3, 


(1-29) 


Note  that  a  factor  of  k 2  has  been  removed  from  the  original  quartic  equation  because  this 

r *  n 

yields  the  trivial  solution  k  =  0.  A  solution  of  equation  (1-29)  is  given  by 


k  -M. 

rvn-  2 h 


1  + 


1  + 


J4 (VW 

(tm) 2 


1/2 ' 


^  mr_  , 

-  h  m i 


Aek0h  «!;//  = 


1  2  3 

i,  ,  ...  . 


(1-30) 


The  above  result  corresponds  with  expression  (88)  on  page  90  of  Galejs  (reference  2).  Note  that 
the  first  term  in  the  above  expression  corresponds  to  the  cutoff  wave  number  for  the  TM  modes 
in  the  parallel-plate  waveguide  with  perfectly  conducting  walls  as  given  in  formula  (2-15). 
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To  determine  the  thin-shell  approximations  for  the  excitation  factor  Aen ,  the  height-gain 
function  Gen(z) ,  and  the  z-dependent  normalized  impedance  Aen{z)  for  the  TM  modes  as  defined 
in  expressions  (F-14),  (F-15),  and  (F-16),  respectively,  one  must  first  obtain  the  thin-shell 
approximation  for  the  constant  Bev  that  appears  in  the  radial  function  Re  that  is  defined  in 
formula  (C-10).  The  constant  Be  is  given  by  either  formula  (C-7a)  or  (C-7b).  The  substitution 

*  n 

of  the  thin-shell  approximations  (1-1 1),  (1-13),  and  (1-14)  into  formula  (C-7a)  yields 


n  =  0, 1, 2, ... 


a-3i> 


Therefore,  the  thin-shell  approximation  to  the  radial  function  Re  is 


K  («)  * 


1+7 

(v-v.* 

e-jkrv(r-2a) 

1  -7*| 

r 

n  =  0, 1,2,.... 


a-32) 


The  thin-shell  approximation  for  the  height-gain  function  Gen{z)  is  obtained  through  substitution 
of  the  approximation  (1-32)  into  (F-15)  and  is  given  by 


Gen(z)  m  cos  krv  z.  - 


i _ 


krv,a 


sin  k  z  ,  n  =  0, 1,2, 


(1-33) 


Note  that  in  the  above  expression,  r  -  a  has  been  replaced  by  z,  the  vertical  distance  above  the 
earth.  At  ELF,  the  second  term  in  the  above  expression  is  much  less  than  the  first  term,  resulting 
in  the  following: 


Gen(z)  =  cos  kr 


■  Agko  ■  , 

z  +  J  —jr  sm  k 


rv 


Z  ,  n  =0,  1,2, 


a-34) 


The  above  result  corresponds  with  expression  (93)  on  page  93  of  Galejs  (reference  2). 


II  II 0 

Rev  ,  is  obtained 

through  substitution  of  the  approximation  (1-32)  into  the  definition  (F-12)  and  is  given  by 
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ra  +  h 


Re  z  -  k 

v  o 


Jkrv / 


Ja 


1 +j[Kv„«-\K«) 


l-j[krva  +  Akoa 


e-jkrv(r-2a) 


dr 


lJL  J2kn  “ 

\  ~  2k 


rv„ 


(ej2kr'»h  -l)-2h 


1  + 


j[krva-Ask0o 


~{k  n-,“ 


+  Agkoa 


+ 


2k. 


1  +j[krva-Ak0a 


l-j[krva  +  Ak0a 


(e-J2k"  J’  -  lH  ,  «  =  0,1,2, 


(1-35) 


If  the  square  of  approximation  (1-32)  is  evaluated  at  r  =  a,  then 


w 


j2kr a 

2  e  r'» 


r2 


1  -2 


1  +j[krva-AJc0a 


\-j[krva  +  Agk0a* 


+ 


+  j(kry,f-Askoa) 


a  +  Agkao 


n  =  0, 1,  2, .... 


(1-36) 


The  thin-shell  approximation  for  the  TM  mode  excitation  factor  Aen  is  found  through  the 
substitution  of  approximations  (1-35)  and  (1-36)  into  the  definition  (F-14),  resulting  in 


sin  2 k  h  /sin  2 k  h 

Aes*Jl  + - — 2_  +  - — —  i 

«  t  2k...  h  2k, „  h 


'Vii2 

v  I 


(knA 

V  n  ) 
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P-\K 


k2rva 


+  [  1  -  cos  2 krv  h 


+ 


k2  h 


,  w  =  0,  1,2, 


a-37) 
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Note  that  in  the  ELF  band,  krv  h  ->  0 ,  |  Ag  | «  1,  and  |  Ae\ «  1,  and,  therefore,  the  above  result 
reduces  to  Aen  =  0.5. 


The  z-dependent  normalized  impedance  for  TM  modes,  Aen(z) ,  has  been  defined  in  formula 
(F-16).  The  substitution  of  definition  (C-10)  and  the  thin-shell  approximations  (1-13),  (1-14),  and 
(1-31)  into  the  numerator  of  formula  (F-16),  in  appendix  F,  yields 


jL 

du 


uRev(u)  *(l  +jkrvr)j 


}-jkrj) 


e-jkrVn(r-2a) 

f 


a-38) 


The  substitution  of  approximations  (1-32)  and  (1-38)  into  expression  (F-16)  gives 


—  krvn  C Z+a)Agk0a  +  jz 

cos  k  z-j 

n 

1  + 

^vn(^+«))i 

sin  krvz 

vn 

k0(z+a) 

krv  a  cos  krvz-[ 

\-jAgk0a 

)  sin  krvz 

n  =  0,1,2,....  (1-39) 

Note  that  for  points  on  the  surface  of  the  earth,  z  =  0,  the  normalized  impedance  reduces  to 

^(0)  =-4. 

1.3  APPLICATION  OF  THIN-SHELL  APPROXIMATION  TO  TE  MODES 


Consider  the  characteristic  equation  for  the  TE  modes  as  given  in  expression  (C-21).  As 

mentioned  in  appendix  C,  the  TE  characteristic  equation  is  equivalent  to  the  TM  characteristic 
equation  (C-8)  if  Ag  is  replaced  by  1  /Ag  and  Ae  is  replaced  by  1  IAh.  Therefore,  with  the 

appropriate  substitutions  applied  to  the  approximate  TM  mode  characteristic  equation  (1-16),  the 
TE  mode  characteristic  equation  is  given  by 


k r9 m^h  +  ko\ 
krVaAh-ko) 


J2krv  h  „  1  ^jltnjT 

e  m  =  1  =  e 


m~  1,2,3, .... 


(MO) 


In  the  above  equation,  the  index  m  corresponds  to  the  order  of  the  TE  mode.  Note  that  there  are 
an  infinite  number  of  solutions  to  the  above  equation. 
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The  TE  mode  characteristic  equation  (1-40)  may  be  rewritten  as 


=  «=  1,2,3,..., 

\l~Km  \  1-tfl 


where  if  and  £  are  defined  as 


if,,,  3 


k  -  A 

rvm  8 


m~  b  > 
Ko 


and 


^  rv  A  h 


(1-41) 


(1-42) 


(1-43) 


In  the  above  expressions,  note  that  a  subscript  has  been  added  to  the  radial  wave  number  to 
indicate  the  mode  index.  The  natural  logarithm  of  equation  (1-41)  yields  the  following: 


In 


1  +  if 


m 


K. 


+  In 


1  +  f 


1  -  em 


m  - 


=  -  j 2  k  h  -  Tim  ,  m  =  1 , 2, 3, ... . 


(1-44) 


In  the  ELF  band,  because  the  spherical  waveguide  boundaries  appear  as  nearly  perfect  electric 
conductors,  then  \Ag\  «  1  and  \Ah\  «  1,  resulting  in  the  following  conditions: 


«  1, 


fl-45) 


The  application  of  the  approximations  (1-22)  and  (1-45)  into  the  TE  mode  characteristic  equation 
(1-44)  results  in 
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m 


if3 

+  — +  t 
3 


C3 

^>m 


+  “f1  s  “  j[krv  h  ~  mjr)  »  m  =  2,  3,  ...  . 


(1-46) 


The  substitution  of  the  variable  definitions  (1-42)  and  (I-  43)  into  the  approximate  TE  mode 
characteristic  equation  (1-46)  yields  the  following  cubic  equation  for  k r-  : 


'(^  +  *2)' 

k3  4- 

3  kl 

*  r9  + 

rn 

(4  +  ^/,) 


+  jh 


krv  ~  imK  ’  m~  2’3’ 


(1-47) 


The  application  of  the  conditions  (1-45)  to  the  above  equation  results  in  a  negligible  cubic  term  in 
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k  resulting  in  the  following  approximate  solution: 

v  m 


mix 

h 


1- 


j  (Ag  +  Ah) 


kji 


-l 


,  m  =  1,2,3, 


a-48) 


The  above  result  corresponds  with  expression  (90)  on  page  90  of  Galejs  (reference  2).  Note  that 
as  Ag  -*  0  and  Ah—>0,  then  kri>  -*  irnilli ,  corresponding  to  the  cutoff  wave  number  for  the  TE 
modes  in  the  parallel-plate  waveguide  with  perfectly  conducting  walls  as  given  in  formula  (2-32). 
Also  note  that  if  m  =  0,  then  jfc  =0  and  the  TE  fields  are  not  excited. 


To  determine  the  thin-shell  approximations  for  the  excitation  factor  Ahm ,  the  height-gain 
function  G*(z) ,  and  the  z-dependent  normalized  impedance  A*(z)  for  the  TE  modes  as  defined 
in  expressions  (G-13),  (G-14),  and  (G-15),  respectively,  one  must  first  obtain  the  thin-shell 
approximation  for  the  constant  B\  that  appears  in  the  radial  function  R]l  that  is  defined  in 
formula  (C-23).  is  given  by  either  formula  (C-20a)  or  (C-20b).  The  substitution  of  the  thin- 

Vm 

shell  approximations  (1-11),  (1-13),  and  (1-14)  into  formula  (C-20a)  yields  the  following: 


bj  =- 


\  +  }{kr*.aAl-k<p) 

\-j(kr*.a\  +  Ka) 


J2kr*a 


m=  1,  2,  3, 


a- 49) 


Therefore,  the  thin-shell  approximation  to  the  radial  function  R.r  is 


Jkr 


K  («) 


Ag  +  j{kr9aAg-k0a 


AR-i\kryJaA^koa 


~jkrv  (r  ~2a) 


,  m  =  1,2,3, 


a-50) 


The  thin-shell  approximation  for  the  height-gain  function  G*(z)  follows  directly  from  Gen(z) 
with  Ag  replaced  by  1  /Ag.  Therefore,  from  formula  (1-34), 


Gfn(z)  ■  cos  krV  z  +j 


Askryn 


sin  A: 


rv„ 


m  =  1,2,3, 


a-51) 


The  above  result  is  applicable  at  ELF  and  corresponds  with  expression  (99)  on  page  94  of  Galejs 
(reference  2). 
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The  thin -shell  approximation  for  the  TE  mode  excitation  factor  Ahm  follows  directly  from 
Aen  with  Ag  is  replaced  by  VAg.  Therefore,  from  formula  (1-37), 


A"  «■ 

m 


A  h 

zl Z^JE 
ko  J 


{l V 

\koht 


1  + 


V 


1  - 


k0aj 


sin  2k  h 

_ 1  tn 

~krvJl 


1- 


koa 


-  1 


koa 


(Ak  -\2( 
8  rv* 

ko  J 


1  + 


l 


sin  2k  „  h 

' 1  tn 

2Kr9h 


1  -  cos  2k  rX,  h 


,  m=  1,2,3,.... 


(1-52) 


Note  that  in  the  ELF  band,  k,  h-*  0  and  with  the  conditions  (1-45),  then  A"  ->  0.  Therefore, 

*  tn 

the  TE  modes  are  weakly  excited  in  the  ELF  band.  The  2-dependent  normalized  impedance  for 
the  TE  modes,  A%t(z) ,  is  obtained  through  inversion  of  the  formula  (1-39)  for  Aen(z)  and  with  Ag 
replaced  by  1  /Ag  and  is  given  as  follows: 


-  k0(z+a) 

A  k  a  cask  z- 

o  '  '  m  m 

K  -ps, 

1  sin  z 

'  '  m 

> 

3 

(z+a)koa  +  jzAg 

cos  k  z  +  j 

m 

A  +  A  i 

s 

[k^z+a)) 

{Kr°)-Jko° 

sin  krTl  z 

'  m 

m=  1,2,  3, ... . 


(1-53) 


Note  that  for  points  on  the  surface  of  the  earth,  z  =  0,  the  normalized  impedance  reduces  to 

4*(0)—. 4r 
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APPENDIX  J 

LEGENDRE  FUNCTION  Pv(-cos  0)  AND  ITS  DERIVATIVES 


The  range  dependencies  of  the  spherical  wave  propagation  formulas  derived  in  this  report 
involve  the  Legendre  function  of  the  first  kind  Pv  and  its  first  two  derivatives.  In  this  appendix, 
both  exact  and  approximate  formulas  for  Pv  are  presented.  The  exact  formula  is  expressed  in 
terms  of  an  infinite  series.  The  approximate  formulas  for  Pv  are  restricted  to  particular  ranges  of 
the  argument.  Comparisons  of  the  exact  and  approximate  formulas  for  Pv  are  given  for  several 
values  of  vthat  correspond  to  particular  propagation  conditions. 

J.l  ASYMPTOTIC  APPROXIMATION 

A  traveling  wave  representation  of  the  spherical  wave  formulas  derived  in  this  report  is 
obtained  from  the  first  term  in  the  asymptotic  series  for  Pv(-  cos  0).  From  Erdelyi  (reference  28), 
if  |  V  |  »  1 ,  Im{  v}  >  0,  and  0  is  not  near  0  or  n,  a  suitable  asymptotic  approximation  for  Pv(-  cos  0) 
is 


Pv(-  cos  0) 


3i  (v  +  1/2)  sin  0 


1/2 


COS 


(v  +  l/2)(3i  -  0)  -  ji/4 


(J-l) 


If  the  cosine  is  expressed  as  the  sum  of  two  complex  exponentials,  the  above  approximation 
becomes 


PS-  cos  0) 


1 


11/2 


2ji  (v  +  1/2)  sin  0 
+  exp  -j(v+  l/2)(jt  —  0)  +  731/4  \  . 


exp  j(v  +  l/2)(3t  -  0)  -  /jc/4| 


(J-2) 


The  asymptotic  formula  (J-2)  may  be  rewritten  as  the  superposition  of  two  waves  traveling  in 
opposite  directions  as 


PS-  cos  0) 


1 


j2nka  sin  0 


1/2  ejk<m  \e-jkp  +  je-jkPi  1  , 


(J-3) 


where 


k  =  p-ja  =  vV/2 


(J-4) 


J-l 


In  the  above  formulas,  k  is  the  wave  number  in  the  earth-ionosphere  waveguide,  /J  =  2n!X  and  a 
are  the  phase  and  attenuation  constants,  respectively,  and  X  is  the  wavelength.  In  addition,  p  = 
a6  and  p,  =  2k a  -  p  are  the  direct  and  indirect  great-circle  path  distances  from  the  source  to  the 

receiver,  respectively,  as  shown  in  figure  4.  Note  that  there  is  a  90°  phase  advance  acquired  by 
the  indirect  great -circle  path  wave  as  it  passes  through  the  antipode  (6-n). 

As  shown  in  appendices  F,  G,  and  H,  the  range  dependencies  of  many  of  the  field  expressions 
involve  either  the  first  or  the  second  derivative  of  Pv(-  cos  0)  with  respect  to  0.  From  equation 
(J-3),  the  asymptotic  formulas  for  the  first  and  second  derivatives  of  Pv(-  cos  0)  with  respect  to 
0are 


Jq  Pv(-  cos  6)~  - 


jka 

1/2 

jkcm 

2 Jt  sin  6 

C 

e-Pp_je~M 


and 


d^_ 

dd2 


P(-  cos  0)  ~  - 


1 


2jn  sin  6 


1/2 


(ka)3'2  ejkan 


e~PP  +  je-Jkp, 


(J-5) 


(J-6) 


The  above  formulas  are  based  on  taking  only  the  derivative  of  the  exponential  sum  in  expression 
(J-3)  because  this  is  the  dominant  spatially  varying  part  of  the  asymptotic  formula  for  Pv(-  cos  6). 

Note  that  there  is  a  90°  phase  advance  in  the  indirect  wave  in  the  first  derivative  formula  and  a 
90°  phase  decrease  in  the  indirect  wave  in  the  second  derivative  formula. 


The  asymptotic  formulas  (J-3),  (J-5),  and  (J-6)  are  useful  in  the  determination  of  the  relative 
phase  differences  between  the  direct  and  indirect  great-circle  path  fields.  For  ELF  propagation 
formulas  based  on  the  earth-flattening  approximation,  each  indirect  great-circle  path  field 
component  is  obtained  from  the  corresponding  direct  great-circle  path  component  through 
replacement  of  p  by  p,  and  with  an  appropriate  phase  adjustment  as  obtained  from  the  asymptotic 

formulas  (J-3),  (J-5),  and  (J-6).  From  these  asymptotic  formulas,  the  phase  adjustments  for  the 
indirect  great-circle  path  fields  are  given  in  the  table  J-l.  Note  that  the  spherical  components  r 
and  0  of  the  fields  are  replaced  by  the  equivalent  cylindrical  components  2  and  p,  respectively. 
This  replacement  is  made  because  table  J-l  is  applied  to  the  formulas  based  on  the  earth¬ 
flattening  approximation  that  are  presented  in  cylindrical  coordinates  in  section  3. 


J-2 


Table  J-l.  Phase  Adjustment  Terms  for  the  Indirect  Great-Circle  Path  Fields 


Source 

Field  Component(s) 

Primary  Range 
Dependence 

Phase  Adjustment 
Term 

VED 

4 

pv(-  COS  9) 

j 

Ep,Hv 

■JfjPyi-  COS  9) 

~j 

HED 

4  4p ,  Hp 

JqPv(~  cos  9) 

-j 

Ep ,H(p 

&**-"*•> 

j 

J.2  APPROXIMATE  SERIES  FORMULA  SUITABLE  IN  THE  VICINITY  OF 
THE  ANTIPODE 


In  the  previous  section,  the  asymptotic  formulas  for  Pv(-  cos  ff)  and  its  first  two  derivatives 
are  not  valid  for  9  near  0  or  n  (i.e.,  for  field  points  located  near  the  source  or  the  antipode). 
Therefore,  these  formulas  are  not  suitable  for  observation  points  lying  in  the  vicinity  of  the 
antipode.  MacDonald  (reference  29)  has  derived  an  approximate  formula  for  Pv  that  is  valid 

when  its  argument  is  at  or  close  to  one  and  is,  therefore,  valid  in  the  vicinity  of  the  antipode. 
MacDonald’s  formula  for  Pv(-  cos  9)  is  given  as 


Pv(-  cos  9)  =  J0(r])  +  sin: 


jt  -  0 


O 

+ 

f . . 

sinf-f) 

-  - 

,  (J-7) 


where 

r]  =  2(v  +  1/2)  sin(^^)  =  2 ka  sin[^y^)  .  (J-8) 

In  expression  (J-7),  Jn  denotes  the  Bessel  function  of  the  first  kind  and  order  n.  Note  that 

formula  (J-4)  is  applied  in  expression  (J-8).  Expression  (J-7)  can  be  simplified  through  the  use 
of  the  following  recurrence  relation  (reference  24) 

Wz)  =  x-/»(z)-',»-i<z).  P-9> 


where  n  is  a  positive  integer  or  zero,  and  z  is  an  arbitrary  complex  constant.  The  application  of 
the  above  relation  to  the  first  bracketed  term  in  formula  (J-7)  yields 


J-3 


2^  J3(i])  =  |  J0(i])  -  (l  +  ?72)  JjC?;) 


6*j 


(J-10) 


For  observation  points  in  the  vicinity  of  the  antipode,  {n-  6)1  2  «  1  and 

^  n  —  8 1  s  jt  -  6 


sin 


(J-ll) 


As  a  result,  r\  can  be  approximated  as 

nskPa> 


(J-12) 


where  pa  =  a{n  -  0)  is  the  great-circle  path  distance  from  the  antipode  to  the  field  point  (figure 

4).  The  substitution  of  expression  (J-10)  and  approximations  (J-ll)  and  (J-12)  to  the  two  leading 
terms  in  the  series  (J-7)  yields 


P v(~  cos  0)  —  J0(kpa) 


1  + 


1  Pa 


12  v  a 


pja_ 
24 ka 


l+(kpa)2}j,(kpa)  ,  pa«l  .  (J-13) 


For  field  points  in  the  vicinity  of  the  antipode,  the  first  term  in  the  series  (J-13)  is  dominant. 
Therefore,  for  small  pa,  consider  only  the  leading  term  of  the  above  series,  i.e., 


P v(-  cos  d)ej0(kpa) 


1  +  _L  \£a 

12V  a 


(J-14) 


In  the  above  approximation,  the  second  term  within  the  brackets  accounts  for  the  curvature  of  the 
earth  as  pala  increases.  From  Wait  (reference  1),  the  bracketed  term  in  formula  (J-14)  can  be 

approximated  as 


\  +  J_(^) 

12V  a  ) 


Pala 


sin  pjo 


1/2 


(J-15) 


where  the  error  in  the  in  the  above  approximation  is  O  [( pala  )4].  The  substitution  of  the  above 
approximation  into  formula  (J-14)  yields 


Pv(-  cos  e )  —  j0(kpa) 


pja 


sin  palo 


1/2 


(J-16) 


J-4 


In  the  above  formula,  J0(kpa)  corresponds  to  the  range  dependence  for  a  flat  earth  and  the 
square-root  term  accounts  for  the  curvature.  Therefore,  formula  (J-16)  is  an  earth-flattening 
approximation  to  Pv(-  cos  9)  that  accounts  for  curvature  and  is  valid  for  field  points  located  in 
the  vicinity  of  the  antipode. 


The  approximation  (J-16)  has  been  applied  in  the  development  of  some  antipode-centered 
ELF  propagation  formulas  that  were  previously  derived  by  the  author  (reference  7).  From 
formula  (J-16),  the  first  and  second  derivatives  of  Pv(-  cos  9)  with  respect  to  6  are 


jg  Pv(-  cos  9)  s  ka  J^kpJ 


Pja 


sin  pja 


1/2 


and 


d2 


■J^2  Pv(~  COS  B)-~  (kaY 


Pa/a 


sin  pla 


1/2 


(J-17) 


(J-18) 


The  approximate  formulas  for  Pv(-  cos  9)  and  its  first  two  derivatives  as  given  in  expressions  (J- 
16),  (J-17),  and  (J-18)  are  required  in  order  to  reduce  the  vertical  electric  dipole  (VED)  and 
horizontal  electric  dipole  (HED)  spherical  waveguide  formulas  to  the  corresponding  antipode- 
centered  ELF  propagation  formulas  presented  in  section  4.3. 


J.3  INFINITE  SERIES  FORMULA 


The  Legendre  function  of  the  first  kind  and  degree  v  and  order  zero,  Pv,  is  related  to  the 
hypergeometric  series  Fas  (reference  24) 


Pv(z)  =  F  -v,  v+1;  1; 


1  -z 


(J-19) 


where  F  is  given  as 


F(n  h.  r.  _  V  («)».(*)•  zn  _  F(C)  V  r(a+n)  r(b+n )  zn 
’  ,,J  (c)„  nl  ~  r(a)  r(b)  „40  r(c+n)  n\  ‘ 


(J-20) 


In  the  above  series,  F  denotes  the  gamma  function  and  ( •  )n  denotes  Pochhammer’s  symbol 
defined  as 


J-5 


(z),= 


z(z  +  l)(z  +  2) 


1 ,  n  -  0 

(z  +  n-  1) ,  »=  1,2,3,... 


(J-21) 


With  the  recurrence  formula  for  the  gamma  function  given  by  T(z+1)  -z  T\z),  Pochhammer’s 
symbol  can  be  expressed  in  terms  of  the  gamma  function  as 


(z)„ 


r(z  +  n ) 

m 


n  =  0,  1,  2, ... . 


(J-22) 


According  to  Abramowitz  and  Stegun  (reference  24),  the  hypergeometric  series  (J-20)  converges 
for  \z  |  <  1.  From  expressions  (J-19)  and  (J-20),  Pv(-  cos  6)  is  given  by 


f 


PA-  cos  d)  =  F 


— v,  v  +1 ;  1; 


1  +  cos  6  |  _  ^  (— v) n  (v  +l)n  (1  +  cos  Q)n 


n  =  0 


2"  («!)2 


(J-23) 


Therefore,  in  the  evaluation  of  Pv(-  cos  6),  the  hypergeometric  series  (J-23)  converges  for 
0<  9=  n. 


The  derivative  of  Pv(-  cos  6)  with  respect  to  6  is  given  as 

jq  Pv(-  cos  6 )  =  P'v(-  cos  6)  —  jq  -  =  Pv(-  cos  0)  sin  6 ,  (J-24) 

where  the  prime  denotes  the  derivative  with  respect  to  the  argument  -  cos  6.  From  equation 
(J-24),  the  second  derivative  of  Pv(-cos  0)  with  respect  to  6  is  given  as 


jLpv(- cos  e)  =  ± 


PA-  cos  6)  sin  6 


=  Pv(-  cos  6)  cos  6  +  sin  6  P'v(-  cos  6) 


=  Pv(-  cos  6)  cos  6  +  Py(—  cos  6)  sin  2 6  . 


(J-25) 


From  the  series  formula  (J-20),  the  first  and  second  derivatives  of  the  hypergeometric  series  are 

4;F{a,b\c\z)  =  g£-F{a+\,b+\\c+\\z)  ,  (J-26) 

and 


£_ 

dz2 


F(a,  b\  c\  z)  = 


(a) 2(fr)2 
(c)2 


F(o  +2,  b  +2;  c  +2;  z)  . 


(J-27) 


J-6 


From  equations  (J-19),  (J-26),  and  (J-27),  the  first  and  second  derivatives  of  Pv(-cos  0)  with 
respect  to  the  argument  -  cos  0  are  given  as 


P’v(-  cos  0 )  =  V(^+1)  f|-v  +1,  v  +2;  2;  *  +f-sd 


and 


P'v(-  cos  0)  =  (  -V---2  ^V'hl)2  f|-v  +2,  V  +3 ;  3 ;  ^-|os  . 


(J-28) 


(J-29) 


The  substitution  of  equations  (J-20)  and  (J-28)  into  (J-24)  yields 


^  Pv(-  cos  0)  =  sin  6  2 

^  n  -  0 


dd 


(-v+l)H  (v  +2)n  (1  +  cos  0)” 
2"  «!  («  +1)! 


(J-30) 


Similarly,  the  substitution  of  equations  (J-20),  (J-28),  and  (J-29)  into  (J-25)  yields 


d2  n  r  a\  V(v+1) 
~2PV(-  cos  9)  =  ^— 


cos  6  2 

n  =  0 


(-v)2(v+l)2  .  2 
+ - - -  sinz 


(-v+l)y>(v+2)n(l+cos  QY_ 
2”  n\  ( n  +1)! 

fl2  ^_v+2^^v+3)«(1  +cos 


»  =  o 


2"  «!  (n  +2)! 


(J-31) 


In  the  evaluation  of  the  spherical  wave  propagation  formulas  derived  in  this  report,  the  Legendre 
function  and  its  first  two  derivatives  are  computed  from  formulas  (J-23),  (J-30),  and  (J-31). 


Before  computing  the  Legendre  function  and  its  first  two  derivatives  from  the  series  formulas 
given  above,  a  discussion  of  numerical  convergence  is  in  order.  As  an  example,  consider  the 
propagation  parameters  for  the  quasi-TEM  mode  at  76  Hz  that  are  given  in  table  5-1.  Under 
daytime  conditions,  v  =  1 1 .7  +  j  0.880.  Table  J-2  provides  a  listing  of  the  number  of  iterations 
required  for  the  series  formulas  to  convergence  to  15  decimal  places  for  various  values  of  6.  The 
table  shows  that  the  number  of  iterations  required  for  each  function  to  converge  increases  with 
distance  from  the  antipode  (0=  7t).  As  0->  0  (source  point),  the  number  of  required  iterations 
increases  rapidly.  This  behavior  is  expected  because  formulas  (J-23),  (J-30),  and  (J-31)  are 
ascending  series  in  the  argument  (1  +  cos  0).  In  summary,  for  0not  too  close  to  zero,  not  much 
CPU  time  is  required  in  order  for  these  series  formulas  to  converge. 


J-7 


Table  J-2.  Number  of  Iterations  Required  for  the  Infinite  Series  Formulas  for  Pv(-  cos  6)  and 
its  First  Two  Derivatives  to  Converge  to  15  Decimal  Places  with  v  =  11.7+ j  0.880 


0  (deg) 

Pv(-  cos  6) 

35pv<-CO! s»> 

rfe2Pv(-cosff) 

10 

3414 

3827 

4240 

30 

423 

454 

486 

50 

160 

168 

176 

70 

84 

87 

89 

90 

51 

52 

53 

110 

34 

34 

35 

130 

24 

24 

23 

150 

17 

16 

16 

170 

11 

10 

10 

The  series  formulas  for  the  Legendre  function  and  its  first  two  derivatives  can  be  used  to 
check  the  accuracies  of  the  approximate  formulas  (J-16),  (J-17),  and  (J-18).  Figures  J-l  show 

plots  of  the  real  and  imaginary  parts  of  the  Legendre  function  and  its  first  two  derivatives  as 
computed  by  the  series  and  approximate  formulas  given  as  functions  of  the  polar  angle  Qa 

measured  from  the  antipode,  where  6a  =  n  -  6. 

The  ranges  of  validity  of  the  approximate  formulas  for  the  Legendre  function  and  its  first  two 
6  derivatives  as  given  in  expressions  (J-16),  (J-17),  and  (J-18)  can  be  determined  through 
comparison  with  the  infinite  series  (exact)  formulas.  Figures  J-la  and  J-lb  show  comparisons  of 
the  magnitude  and  phase,  respectively,  of  the  exact  and  approximate  formulas  for  Pv(-cos  6)  as  a 
function  of  6a  for  v  =  4.75  -  j  0.440.  This  value  of  v  corresponds  to  typical  daytime  propagation 
conditions  at  30  Hz  as  given  in  table  5-1.  Figures  J-2  and  J-3  show  similar  comparisons  for 
dPv(-  cos  9)1  dd  and  d2Pv(-  cos  6)/dd2 ,  respectively.  The  plots  show  that  the  exact  and 
approximate  formulas  are  in  excellent  agreement  out  to  6a  =  45°.  Beyond  this  angle,  the 
approximate  formula  for  the  second  derivative  starts  to  slowly  depart  from  the  exact  formula. 

For  example,  at  Qa  =  45.1°,  the  approximate  formula  for  the  second  derivative  differs  by  0.067 
dB  in  magnitude  and  0.95°  in  phase  with  the  exact  formula.  In  comparison,  at  6a  =  90°,  the 
approximate  formula  for  the  second  derivative  differs  by  0.46  dB  in  magnitude  and  8.32°  in 
phase  with  the  exact  formula.  At  6a  =  170.1°,  the  approximate  formula  for  the  second  derivative 
differs  by  2. 1 8  dB  in  magnitude  and  69.6°  in  phase  with  the  exact  formula.  The  approximate 
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formula  for  dPv(-  cos  6)/dd  shows  better  agreement  with  the  exact  formula  than  the  second 
derivative  and  the  approximate  formula  for  P^-cos  6)  shows  an  even  better  agreement.  This 
observation  is  attributed  to  the  fact  that  each  successive  derivative  of  the  approximate  formula 
for  P^-cos  0)  introduces  additional  error  into  the  result.  The  approximate  formula  for  P^- cos  6) 
does  not  start  to  noticeably  depart  from  the  exact  formula  until  6a  =  150°,  where  the  formulas 
differ  by  0.14  dB  in  magnitude  and  2.26°  in  phase. 

Comparison  plots  of  the  exact  and  approximate  formulas  for  Py(- cos  0),  dPv(-  cos  8)1  dd , 
and  d2Pv(-  cos  8)1  dd2  as  a  function  of  Qa  are  given  in  figures  J-4,  J-5,  and  J-6,  respectively,  for 
v=  15.9  -  j  1.32.  This  value  of  v  corresponds  to  daytime  propagation  conditions  at  100  Hz  as 
given  in  table  5-2.  In  these  plots,  the  oscillations  near  the  antipode  are  considerably  damped 
because  of  the  larger  attenuation  at  this  higher  frequency.  Each  of  these  plots  show  a  noticeable 
improvement  in  agreement  as  compared  to  the  previous  value  of  v.  For  example,  at  6a  -  90°,  the 
approximate  formula  for  the  second  derivative  differs  by  0.012  dB  in  magnitude  and  1.88°  in 
phase  with  the  exact  formula;  and  at  0a  =  150°,  the  approximate  formula  for  the  second 
derivative  differs  by  0.065  dB  in  magnitude  and  6.47°  in  phase  with  the  exact  formula.  As  in  the 
previous  set  of  plots,  the  approximate  formula  for  dPv(-  cos  8)1  dd  shows  better  agreement  with 
the  exact  formula  than  the  second  derivative,  and  the  approximate  formula  for  Pv(-  cos  6)  shows 
an  even  better  agreement.  Further  examination  of  the  approximate  formulas  at  other  ELF 
frequencies  indicates  that  the  accuracies  of  the  approximate  formulas  for  Pv(—  cos  6)  and  its 
derivatives  improve  with  increasing  frequency. 
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Figure  J-l.  Comparison  of  the  Infinite  Series  (Exact)  and  Approximate  Formulas 
for  Pv(-  cos  6)  for  v=  4.75  -  j  0.440;  (a)  Magnitude;  (b)  Phase 


Figure  J-2.  Comparison  of  the  Infinite  Series  (Exact)  and  Approximate  Formulas 
for  dPv(-  cos  6)/dd  for  v=  4.75  -  j  0.440;  (a)  Magnitude;  (b)  Phase 


9a  (deg) 

Figure  J-3.  Comparison  of  the  Infinite  Series  (Exact)  and  Approximate  Formulas 
for  d2Py(-  cos  6)/dd2  for  v=  4.75  - j  0.440;  (a)  Magnitude;  (b)  Phase 


Figure  J-4.  Comparison  of  the  Infinite  Series  (Exact)  and  Approximate  Formulas  for 
PA-  cos  6)  for  v- 15.9  -j  1.32;  (a)  Magnitude;  (b)  Phase 


Figure  J-5.  Comparison  of  the  Infinite  Series  (Exact)  and  Approximate  Formulas  for 
dP  {-  cos  6)/dd  for  v=  15.9 -j  1.32;  (a)  Magnitude;  (b)  Phase 


Phase  { d2P v(-  cos  0 )  /  d0  2 }  (deg ) 


Figure  J-6.  Comparison  of  the  Infinite  Series  (Exact)  and  Approximate  Formulas  for 
d2Pv(-  cos  9)/dd2  for  v- 15.9  -j  1.32;  (a)  Magnitude;  (b)  Phase 


I 
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APPENDIX  K 

EARTH-FLATTENING  APPROXIMATION 


As  shown  in  appendix  A,  the  electromagnetic  (EM)  fields  in  a  spherical  earth-ionosphere 
waveguide  can  be  expressed  in  terms  of  two  Debye  potentials.  At  a  source-free  point,  each  of 
these  potentials  satisfies  the  scalar  Helmholtz  equation.  As  shown  in  appendix  B,  the  Helmholtz 
equation  can  be  solved  in  spherical  coordinates  by  the  separation  of  variables  method  in  which 
the  range  or  6  dependence  is  expressed  in  terms  of  a  function  T(6),  which  satisfies  the  following 
differential  equation: 

0  +  c°te^  +  v(v+i)r  =  °.  (K-i) 

In  equation  (K-l),  v(v+l)  is  the  separation  constant.  As  shown  in  appendix  B,  for  a  source 
located  at  9=  0,  the  solution  of  equation  (K-l)  that  is  finite  at  the  antipode  (0=n)  is  given  by 

T(6)  =  CPV{- cos  0)  ,  (K-2) 

where  Pv  is  the  Legendre  function  of  the  first  degree  v  and  order  zero,  and  C  is  an  arbitrary 
constant.  Both  exact  and  approximate  formulas  for  Pv(-cos  0)  are  given  in  appendix  J. 


In  the  earth-flattening  approximation,  the  polar  angle  0  in  equation  (K-l)  is  replaced  by  the 
range  variable  p,  where  p  =  a6  and  a  is  the  radius  of  the  earth.  Note  that  p  is  the  direct  great- 
circle  path  distance  along  the  spherical  earth  that  connects  the  source  and  the  field  points  (see 
figure  3-2).  The  replacement  of  0by  p/a  in  equation  (K-l)  yields 


d2T 

dp2 


+  dcot 


^  +  s2T  =  0 
dp 
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(K-3) 


where  the  constant  s  is  defined  as 
v(v+l) 

O  —  ^ 


(K-4) 


In  his  estimation  of  the  error  involved  in  the  earth- flattening  approximation,  Pekeris  (reference 
12)  assumed  that  T{6)  can  be  expanded  in  the  following  asymptotic  series: 

oo 

7X0)  =  T(pla)  =  2  (as)~2n  T2n(p)  =  TQ{p)  +(as)~2T2(p )  +  (as)~4  T4(p )  +  ... .  (K-5) 

w  =  0 


K-l 


The  series  expansion  of  cot(p/a)  is  given  by  (reference  27) 


(K-6) 


The  substitution  of  the  series  expansions  (K-5)  and  (K-6)  into  the  differential  equation  (K-3) 
followed  by  the  multiplication  of  s  ~2  yields 


d2T/ 


dx2 


O  ,  ,-2<l2T2  ,  x -4  d2T4 

r  +  (as)  +  (as)  4  — rri1  +  ... 


dx2 


1  1  4- 

X  3 


(as)2  45  (us)~ 


£0 

dx 


+  ( as)' 


dT2 

dx 


+  (as)' 


dTj. 

dx 


+ 


T0+(as)~2T2  +(asr4T4  + 


0  , 


(K-7) 


where  x  =  sp.  Following  Pekeris  (reference  12),  the  coefficients  of  like  powers  of  (as)~2  are 
equated  to  obtain  a  system  of  simultaneous  differential  equations.  The  first  three  equations  are 
given  as 


L(Ta)  =  0  , 

(K-8a) 

L(j)-x  dTo 

3  dx  ’ 

(K-8b) 

and 

y3  dT n  v  dT 2 

^->  =  45  *  +f  d  ’ 

(K-8c) 

where  the  differential  operator  L  is  given  as 

L=  d2  +  1  d  +  1 
dx2  x  dx  + 

(K-9) 

The  solution  of  equation  (K-8a)  is 

T0(p)  =  H&(sp)  .  (K-10) 

The  above  formula  gives  the  range  dependence  for  a  planar  earth  (with  a  source  located  at  p  =  0) 
and  represents  the  earth-flattening  approximation.  To  obtain  the  solution  of  equation  (K-8b),  let 
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us  consider  the  following  differential  equation  that  is  valid  for  any  cylindrical  Bessel  function 
Z„(x)  of  arbitrary  order  n: 

(K-ll) 

The  above  result  was  obtained  from  Koo  and  Katzin  (reference  13).  Note  that  Zn(x)  refers  to  the 
Bessel  function  of  the  first  kind  Jn(x),  the  Bessel  function  of  the  second  kind  Yn(x),  the  Hankel 
function  of  the  first  kind  H^Xx) ,  or  the  Hankel  function  of  the  second  kind  H\^\x)  .  With  the 
assumed  time-harmonic  dependence  of  eim  and  an  outward  propagating  wave  from  the  source, 
Zn(x)  =  H^\x)  in  the  present  application. 


The  substitution  of  Zn(x)  =  H^\x)  and  n  =  2  into  formula  (K-l  1)  yields 

The  derivative  of  x2  H^\x)  is  given  by 

[*2 H$\x)\  =  x2H^'(x)  +  IxHfXx) , 


(K-l  2) 


(K-13) 


where  the  prime  denotes  the  derivative  with  respect  to  the  argument.  The  substitution  of  the 
recurrence  relation  (reference  24) 

HF'(x)  =  H?Xx)-ZHf?Xx) , 

into  equation  (K-13)  yields 


^  [x2  H?Xx)]  =  x2  Hf\x)  =  -x2  H^'(x)  . 

The  substitution  of  the  above  result  into  expression  (K-l  2)  gives 
L  [x2  H(2Kx)}  =  -4x  H^'(x)  . 

If  H[f](x)  is  replaced  by  T0(p),  formula  (K-l 5)  may  be  rewritten  as 


-  jo  x2H{iKx)\^  = 


x*To 
3  dx 


(K-l  4) 


(K-l  5) 


(K-l  6) 


A  comparison  of  expressions  (K-8b)  and  (K-l  6)  shows  that 
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(K-17) 


T2(p)=-±(sp)2H?\sp), 


where  x  has  been  replaced  by  sp. 


The  substitution  of  formulas  (K-10)  and  (K-17)  into  the  asymptotic  series  (K-5)  yields 
T(pla)  *  Hf\sp)  -  -jL  (§)2H>2>(sp)  +  o(§)4  .  (K-18) 

In  the  above  equation,  the  first  term  corresponds  to  the  model  of  a  flat  earth  and  the  remaining 
terms  are  corrections  for  curvature.  Therefore,  the  earth-flattening  approximation  is  given  by 

T(p/a)  *  H^(sp)  .  (K-19) 

The  earth-flattening  approximation  is  valid  when  the  great-circle  path  distance  p  is  small 
compared  with  the  earth’s  radius  a.  In  particular,  Pekeris  (reference  12)  found  that  the  earth¬ 
flattening  approximation  is  correct  to  within  2%  of  the  exact  value  for  ranges  out  to  about  half 
the  radius  of  the  earth.  It  should  be  noted  that  this  result  is  independent  of  the  frequency. 


Next  one  needs  to  determine  the  constant  C  in  formula  (K-2)  so  that  T(6)  can  be  related  to 
Pv(-  cos  6) .  To  obtain  C,  let  0  -» 0  on  each  side  of  expression  (K-2),  giving 

lim  T(8)  =  C  lim  Pv(- cos  8)  .  (K-20) 

e  —  o  e-»o  v 

From  equation  (K-18), 

Hmo  T(6)=  ^lim()  T(pla)  -  H^\sp)  .  (K-21) 

From  Abramowitz  and  Stegun  (reference  24), 

limoH«Xsp)  =  -ji  ln(jp)--yf  lnp  .  (K-22) 

From  appendix  E,  formula  (E-8), 

lim  PI-  cos  9)  =  aiua  |„  e2  =  2mva  |n(g)  ^  Iskva  |n  p  (K.23) 

The  substitution  of  the  limits  (K-21),  (K-22),  and  (K-23)  into  expression  (K-20)  gives 
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(K-24) 


j  sinvn 

From  expressions  (K-2),  (K-18),  and  (K-24), 

=  r<e>  =  r(§)  *  -  £  (§)  •  0«5) 

where  the  error  in  the  above  approximation  is  0(p  la  )4.  From  formula  (J-4),  the  wave  number  k 
in  the  earth-ionosphere  waveguide  is  given  by 

k  =  p-ja  =  ^. 1/2  (K-26) 

If  |  v| »  1,  then 

(v  +  1/2) 2  *  v(v  +  1)  .  (K-27) 

The  above  result  is  generally  valid  in  the  ELF  band.  The  application  of  the  above  approximation 
to  definitions  (K-4)  and  (K-26)  results  in  the  following  approximation: 

5  S  k ,  Iv  I  »  1  .  (K-28) 

Therefore,  expression  (K-25)  can  be  approximated  out  to  order  (p/a)2  as 

fjsi  ‘  HoKkp)  -  n  (§)2«22><*P) .  M  »  1  •  (K-29) 

The  above  formula  is  similar  to  the  one  applied  by  Wait  (reference  1)  in  his  development  of  ELF 
propagation  formulas.  The  formulas  differ  by  a  factor  of  j. 

For  |  kp  |  »  1,  the  Hankel  functions  in  expression  (K-29)  can  be  replaced  by  the  first  term  of 
their  respective  asymptotic  expansions: 

Hfikp)  ~  ,  (K-30a) 

and 

H2)(kp)  ~  e~m>  5nl4)  •  (K_30b) 
The  substitution  of  the  above  approximations  into  formula  (K-29)  yields 


K-5 


Py(-  COS  6) 
j  sin  vji 


,  Ivl  »  1,  Iftpl  »  1  . 


(K-31) 


e~MP  -  rt/4) 


2 


From  Wait  (reference  1), 


-JL.-1  +  §1 

sinO  ~  12  ’ 


(K-32) 


where  the  error  in  the  above  approximation  is  0(0 4).  The  application  of  the  above  approximation 
to  expression  (K-31)  results  in 


Pv(-  cos  0) 

j  sin  vk 


p/a 

sin  pi  a  ’ 


Ivl  »  1,  \kp\  »  1 , 


(K-33) 


where  the  asymptotic  approximation  in  expression  (K-30a)  has  been  replaced  by  H^\kp) .  In 
the  above  approximation  (K-33)  for  Pv(-  cos  0) ,  H^\kp)  corresponds  to  the  range  dependence 
for  a  flat  earth  and  the  square-root  term  is  the  correction  for  curvature.  As  was  noted  following 
approximation  (K-29),  expression  (K-33)  also  differs  with  the  one  derived  by  Wait  (reference  1) 
by  a  factor  of  j.  Formula  (K-33)  is  referred  to  as  the  earth-flattening  approximation  with 
curvature  correction  for  Pv(-  cos  0) . 
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